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APPLICATION OF THE ILKOVIC EQUATION TO 
QUANTITATIVE POLAROGRAPHY 


By Floyd Buckley and John Keenan Taylor 


ABSTRACT 


Conditions of applicability of the Ilkovié equation to polarographic analysis 
were investigated for re presentative ions singly and in mixtures. Gelatin was 
used to suppress maxima found in current-voltage curves. 

Results obtained with a capillary of usual characteristics showed that a sup- 
presser is necessary to secure agreement with the Ilkovié equation over any con- 

iderable drop-time range, and that the need increases with the drop-rate of the 
e electrode and the dilution of reducible ion. A criterion of diffusion control satis- 
factory for classification of polarograms by visua] examination was found. The 
application of the Ilkovié equation to quantitative analysis was shown to be 

nadviss able unless it is definitely established that the experimental conditions are 
well within the range for insuring diffusion control. A logarithmic relation 
etween critical drop-time and suppresser concentration was found. 
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I. INTRODUCTION 


The polarographic method of chemical analysis is based upon the 
interpretation of the current-voltage curves obtained when solutions 
containing electroreducible or electrooxidizable substances are elec- 
trolyzed with electrodes, one of which is very small (usually mercury 
dropping from a fine-bore capillary) so that it is completely polarized, 
while the other is practically unpolarizable. From the position and 
magnitude of certain steps observed in these curves, both qualitative 
and quantitative information concerning the electrolyzed material 
can be obtained. 

The polarograph, an instrument that automatically obtains and 
records these current-voltage curves, has numerous applications in 

chemical analysis. Most of the published procedures make use of the 
instrument as an interpolating device for comparing known and un- 
known solutions. If the standard and unknown solutions are suffi- 
ciently alike, and the same dropping electrode is employed for each 

ander identical conditions, the analysis resolves into a comparison of 
ie heights of the waves or steps on the current-voltage curves. 
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The phenomenon of a limiting current which gives rise to the char. 
acteristic polarographic waves 1s determined, under ideal conditions 
by the rate of diffusion of the electrolyzable ion from an inexhaustib|, 
supply in the surrounding solution to the surface of the electrode 4 
which the concentration is zero because the ions are removed as rapidly 
as they arrive. The expression relating this diffusion current to the 
constants of the system under investigation was derived by Ilkovié 
[1]! in the form 

ig=knD'"?Cm?*¢'", 


where 74 is the diffusion current in microamperes, n is the number of 
electron equivalents per mole of electrode reaction, D is the diffusion 
coefficient of the reducible (or oxidizable) substance, C is its concen- 
tration in millimoles per liter, m is the mass flow of mercury from the 
dropping electrode in milligrams per second, ¢ is the drop-time ir 
seconds, and k is a proportionality constant. In this equation th 
quantity knD'” is a constant for the particular ion under investigation 
whereas the other factor depends upon the concentration of the solv- 
tion and the physical constants of the apparatus employed. Lingan 
[2] has accordingly proposed that the quantity knD'”, which is exper: 
mentally determined by i,/Cm?“t'”, be designated as the ‘“‘diffusion- 
current constant.’”’? The use of this constant, together with th 
physical characteristics of the capillary and of the electrical systen 
would eliminate the need for calibration experiments and put polar- 
ography on an absolute basis. 

The validity of the above procedure depends on strict agreement 
with the Ilkovié equation under the experimental conditions of th 
analysis. However, it is found that the polarographic waves are not 
determined by the diffusion current alone, but are often complicat 
by the presence of more or less pronounced maxima, caused by adsorp- 
tion of the reducible ion on the mercury electrode, stirring of t! 
electrolyte near the electrode, or other factors. Certain surface- 
active materials added to the electrolyte have been found to suppres 
these maxima. Although the Ilkovié equation has been tested by 
several investigators [3] and definite statements have been made 2s 
to the conditions for its validity, none of these researches has included 
a systematic study of the effect of suppressers of maxima on the diffu- 
sion current, but there are general statements in the literature whic! 
describe variations of wave height caused by the addition of a larg 
excess of suppressing agents [4]. Indeed, no satisfactory criteria for 
determining when maxima have been suppressed and the electro 
current is controlled by a simple diffusion process (hereafter referred 
to as diffusion control) have been proposed. Observations in this 
laboratory on the occurrence and nature of maxima have suggeste 
the need of an exploratory investigation to determine the influence: 
of the drop-rate of the capillary and the concentration of the sup- 
presser on the validity of the fundamental Ilkovié equation. This 
information is necessary before Lingane’s proposal can be generally 
applied to quantitative polarography. 

! Figures in brackets indicate the literature references at the end of this paper. 


§ This diffusion-current constant differs from that previously employed by Kolthoff and Lingane 
Rev. 24, 1 (1939) by the factor m2/%t\/* 
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II. APPARATUS AND METHOD 


The polarograph used was of the photographic recording type and 
has been described, together with its auxiliary equipment, in a previous 
communication from this laboratory [5]. The characteristics of the 
electrical measuring circuit were determined by the Electrical Division 
of this Bureau. In addition, the calibration of the galvanometer was 
checked by substituting a standard resistor for the dropping-electrode 
cell. ab acele sensitivity of the galvanometer was 2.2210~° 
ampere per millimeter at 1 meter. Nominal sensitivity ratios differed 
by no more than 1 percent from observed values. 

‘The dropping electrode was constructed from a drawn-out capillary 
so that a maximum drop-time of about 6 seconds was obtained in the 
solutions under investigation. This capillary was sealed to a long 
burette of large cross section so that the driving head was not altered 
significantly during a 10-minute dropping interval. The burette was 
fitted with an enlarged funnel top and a stopcock sidearm near the 
bottom to facilitate adjustment of the driving head to obtain any 
desired drop-time in the range of 1 to 6 seconds. 

The flow of mercury in the capillary was determined as a function 
of the driving head. The apparatus described by Lingane [6] was 
particularly suitable for this measurement. The flow is practically 
independent of the immersion medium (air, water, solution) but the 
drop-time is, of course, quite sensitive to the composition of the sup- 
porting electrolyte. A scale attached to the burette furnished a con- 
venient means of ascertaining the driving head and hence the flow of 
mercury for each polarogram recorded. The drop-time of the capil- 
lary was taken by means of a stopwatch at the half-wave potential of 
the particular ion under investigation. 

Galvanometer sensitivities were selected to produce wave heights 
of such size that errors in their measurement were reduced to about 1 
percent, except for the solutions of extreme dilution, in which case 
the errors were no more than two or three times this amount. 

Characteristic types of polarograms encountered in this investiga- 
tion and the usual methods of measuring wave heights are shown in 
figure l. The segment preceding the wave is extended, as shown in the 
figure, to compensate for rise in residual current. The method for 
estimating the current at the completion of a wave depends upon the 
shape of the wave and the slope of the following segment. In the 
absence of prominent maxima, three types of waves are observed. 
The slope of the following segment may be equal to, as in A, less than, 
as in B, or greater than, as in C, that of the initial segment. In cases 
A and B, a line is drawn parallel to the initial portion of the wave and 
through the maximum of the final portion. Case C permits the back- 
ward extension of the following segment. In each case the wave 
height is taken as the vertical distance at the half-wave potential. 
When sharp maxima are present as in D and E, the peak is ignored 
and estimation of wave height is made as in A and B. 

lhe above-described methods of measuring wave heights are fre- 
quently used and are satisfactory for comparative purposes, but it is 
later shown that they do not provide true measures of diffusion cur- 
rents except when the form of the polarographic wave is like that 
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shown in figure 1, A, as this is the only type of wave obtained when th 
electrode current is determined by asimple diffusion process. This fae; 
forms the basis of a criterion for the minimum concentration of Sup. 
presser necessary to insure diffusion control of the electrode current. 

Nickel, cadmium, bismuth, and thallium were selected as the 
cations for study. This group provided a bivalent ion (Cd) having 
a half-wave potential very near that of the electrocapillary maximun 
one bivalent (Ni) and one univalent ion (Tl) reduced at potentials oy 
either side of the electrocapillary maximum, and one trivalent joy 
(Bi). The selection is a reasonable cross section of the type of cation 
encountered in elementary analysis. The supporting electrolyte 
was 0.1 N in potassium chloride for nickel, cadmium, and thalliuy 


and 1 WN in nitric acid for bismuth. 
i é 








CURRENT 
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VOLTAGE 


Figure 1.— Methods of measurement for characteristic types of polarographic wave 


Solutions of nickel, cadmium, and bismuth were prepared by quat- 
titative dilution of standards made from spectrochemically pur 
metals. The thallium standard was prepared from pure TICl. 

Polarograms for each solution were taken over a drop-time rang’ 
of 1- to 6-second intervals of approximately 0.5 second for nickel ané 
cadmium and approximately 1 second for thallium and bismuth. Th 
gelatin concentration was increased in increments of 0.002 percent 12 
all cases except for cadmium. The range of gelatin concentratio 
was 0 to 0.008 percent for cadmium, thallium, and bismuth and 0 t 
0.020 percent for nickel. 

Oxygen was removed by bubbling purified nitrogen through th 
solution in the cell for 15 minutes preceding the recording of the 
polarogram, and its reentrance was prevented by a continuous flow © 
nitrogen over the solution during the measurement. 

All measurements were moe in a constant-temperature rool 
maintained at 25.5° +0.3° C. 
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III. RESULTS 


Values of the diffusion-current constant [p—i,/Cm?t'” were cal- 
culated from the measured wave heights and the characteristics of 
the apparatus and electrode. The results for various solutions of 
nickel as a function of drop-time of the electrode and gelatin concen- 
tration are given in figure 2. The failure to comply with the law 
expressed by the equation for Jp is strikingly evident in the absence 
of a suppresser in that no constant value of Jp over an appreciable 
range of drop-times is observed. With the small gelatin content of 
0.002 percent, however, constant values of Jp are obtained for all 
drop-times greater than 3 seconds. Ata concentration of 0.004 percent 
of gelatin the critical drop-time is decreased to 2.1 seconds and at 
0.006 percent to 1.8 seconds. Solutions containing between 0.008 
and 0.020 percent of gelatin give, for fixed gelatin concentrations, 
constant values of Jp over the entire drop-time range of 1 to 6 seconds.* 

A representative selection of the polarograms obtained for a solution 
(.142< 107° M in nickel is shown in figure 3. For all the polarograms 
in horizontal alinement, the gelatin content was that indicated at the 
right side of the figure, and for all those in vertical alinement the 
drop-time did not differ by more than 0.1 second from the nominal 
value listed at the top of the figure. Values calculated for the 
diffusion-current constant are inserted at the base of each polarogram. 
At first inspection, all the curves with the exception of the last two 
on the right at 0 percent and the last one on the right at 0.002 percent 
of gelatin appear regular and ‘“‘well defined.’”’ A study of the figure 
reveals the danger of gross error incurred from the assumption of 
diffusion control as determined by the loose criterion of “good 
definition.’”’ A positive criterion that has been found to be sufficient 
in this and all other cases of this study is the near identity of slope 
at the beginning and end of a wave. 

The effect of gelatin, at least for concentrations less than 0.02 
percent, in diminishing the diffusion-current constant, was found to 
be significant only at high dilution of the reducible ion. The family 
of curves shown in figure 4 illustrates the magnitude of this effect 
in solutions of nickel. These curves show that diffusion control 
exists over the entire drop-time range of 1 to 6 seconds for a concentra- 
tion of gelatin in the neighborhood of 0.006 percent. An increase in 
suppresser concentration to 0.014 percent (0.020 percent for solutions 
0.568 X 10-* M and 0.0284 10-* M) gives no further change in Jp for 
the two higher concentrations of nickel but causes a small continuous 
decrease at the higher dilutions. The effect is most pronounced at 
the lowest concentration of nickel. 

The results depicted in figure 4 offer an explanation of conflicting 
statements which have been made regarding the strict proportionality 
of concentration of the reducible ion and its diffusion current. The 
curves show that this proportionality depends both upon drop-time 
' For a given gelatin content, a slight increase in the values for Jp at the higher drop-times is noted, especi- 
aly for the lowest concentrations of the reducible ion. The magnitude of this increase is small and of the 
order of the experimental errors involved in its determination. It may be due to some systematic error 
introduced by the method of measurement of the wave heights. However, the errors incurred on the 
assumption of a constant value for Jp, when diffusion control is effective, should not exceed +2 percent. 
‘The conclusions of this paper are based upon the interpretation of approximately 800 current-voltage 


curves, which obviously can not be published in entirety. The polarograms shown have been selected 
decause they illustrate characteristic behavior. 
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Ficure 2.—-Effect of gelatin and drop-time on the diffusion-current constant of nickel 
in a supporting electrolyte of 0.1 N KCI, 
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and gelatin content. Thus constant values of Jp are obtained over 
the entire drop-time range of 1 to 6 seconds only if the gelatin con- 
centration is very near 0.006 percent. At lower gelatin content the 
proportionality is valid only in a very restricted drop-time range. 
At high gelatin concentrations, proportionality exists only if the ion 
concentration is in the vicinity of 0.1107 M or higher. 

The behavior of cadmium is of special interest as this element is 
reduced at a potential very close to the electrocapillary maximum. 
Both the adsorption theory of Heyrovsky [7] and Ilkovié [8], and the 
electrostatic streaming theory of Antweiler [9] demand that, at this 
potential, maxima be either absent or of small magnitude. Hence, 
deviations from the Ilkovié equation might be expected to be very 
small for this element. Indeed, positive assertions that cadmium 
shows no maximum in a 0.1-N solution of potassium chloride have 
been made [10]. 
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Figure 3.—Characteristic polarographic waves for nickel (0.142 X 10-* M) ina 
supporting electrolyte of 0.1 N KCI. 


All polarograms were taken at a sensitivity of 1/20, except those at a nominal drop-time of 1.1 seconds, 
which were taken at 1/40. 


Figures 5 and 6 show the results fora solution 0.0148107* M 
in cadmium. The general behavior is similar to that of nickel. 
Gelatin is required to secure agreement with the Ilkovié equation 
and deviations at high drop-rates in the presence of insufficient sup- 
presser are showao by the appearance of a decreased slope at the end 
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Ficure 4.—Effect of gelatin on the diffusion-current constant of nickel for specific 
drop-times in a supporting electrolyte of 0.1 N KCl. 











The symbol @ denotes the superposition of points corresponding to the different drop-times. 
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of the wave (figure 6). Sharp maxima are absent and all waves are 
“well defined.”” The writers’ conclusion that the appearance of well- 
defined waves is not sufficient evidence for diffusion control is con- 
frmned. The positive criterion found to be valid for nickel, that is, 
the approximate equality of the slopes preceding and succ eeding the 
wave, was also applicable to all cadmium solutions. 
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Effect of gelatin and drop-time on the diffusion-current constant for 
cadmium in a supporting electrolyte of 0.1 N KCl. 
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Figure 6.—Characteristic polarographic waves of cadmium (0.0148X 10-* M) ina 
supporting electrolyte of 0.1 N KCl. 


The fractions under the curves{denote sensitivities. 
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The deviations shown by cadmium when electrolyzed from solutions 
of potassium chloride of various concentrations and in the absence of 
gelatin are shown in figure 7. All polarograms are regular, ‘well de- 
fined,”’ and show the characteristic appearance of a decreased slope at 
the end of the wave for drop-times less than 3 seconds. Constant 
values of Jp are obtained only for drop-times greater than 3 seconds 
The value of Jp found with concentrations of supporting electrolyte 
less than 1 N is 3.54, which is in excellent agreement with Lingan 
value [2] of 3.51 obtained in a solution 0.1 N in potassium chloride 
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[t is interesting to observe from figure 6 that the average valu 
Ip is 3.35. This low value is in agreement with the results found f 
nickel where low values of Jp are observed for low concentrations 
reducible ion. 

Some results for univalent thallium are given in figures 8 and 
The general behavior of this ion parallels that of nickel in ever 
respect. Both exhibit sharp maxima at high concentrations of re- 
ducible ion (polarograms for nickel under these conditions not show! 
and deviate from ideal behavior at low concentration in the manne! 
already described. The appearance of the two types of abnormalit 
is clearly shown in the polarograms given in figure 9. It is observ: 
that there is a sharp (peak) maximum superimposed on the regul: 
type of polarogram that is obtained at low concentrations of reducib! 
ion. With successive additions of gelatin, these sharp maxima grac- 
ually decrease and eventually disappear, but additional gelatin 1s re- 
quired to destroy the abnormal negative slope and insure diffusio! 
control. The value of Jp for a 1X10-*M solution in the presence ©! 
0.008 percent of gelatin is 2.71, in good agreement with the value 0 
2.69 reported by Lingane and Kolthoff [3]. The decrease of Jp wit! 
extreme dilution of reducible ion previously observed for nickel an¢ 
cadmium is confirmed in the present instance. 
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Effect of gelatin and drop-time on the diffusion-current constant of 
thallium in a supporting electrolyte of 0.1 N KCI. 
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9.—Characteristic polarographic waves of thallium (2x 10-°M) in a sup- 
porting electrolyte of 0.1 N KCI. 
larograms were taken at a sensitivity of 1/150, except those at a nominal drop-time of 1.3 seconds for 


0.002 percent of gelatin and for 2.3 seconds at 0.000 percent of gelatin, which were taken at a sensi- 
1/200. 
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Figures 10 and 11 summarize some results obtained for bismut) 
in 1 N nitric acid. The behavior of this ion differs from that of th, 
preceding in the following ways: (1) Deviations from the Ilkoyjé 
equation in the absence of a suppresser are considerably greater thay 
hitherto noted, (2) prominent rounded maxima appear only at the 
higher drop rates, (3) at concentrations sufficiently large to give very 
pronounced maxima in the range 1 to 3 seconds, no departure from 
regular and well-defined waves is noted for drop-times greater thay 
4 seconds (see fig. 11). In some aspects the behavior of bismuth iy 
an environment of nitric acid is intermediate between that of cadmium 
which shows no prominent maxima, and thallium, which exhibits 
the characteristic peak maximum. The value of Jp for a 1X10~ Vf 
solution is 4.70, in good agreement with the value of 4.64 report 
by Lingane [2]. 

The applicability of the Ikovié equation to the analysis of mixtures 
of reducible ions was studied in the case of Cd—Ni (0.297 and 0.568 
millimolar, respectively) and ZJl-€d (1.00 and 0.897 millimolar 
respectively) in a 0.1-N solution of potassium chloride. The hall- 
wave potentials differ by 0.5 volt for the first pair and only 0.1 voli 
for the second. Hence, in accordance with expectation, the analysis 
of the Cd—Ni mixture was very satisfactory and that of T/-Cd quit: 
uncertain. Two series of polarograms over a drop-time range of | | 
6 seconds at intervals of approximately 1 second were taken for each 
mixture. One series was without gelatin and the other with 0.0! 
percent of the suppresser. Typical polarograms are shown in figur 


) 


In the series of polarograms for the mixture of cadmium and nickel 


in the absence of gelatin, the behavior of cadmium is similar to th 
of the single ion. Nickel differs from its single ion behavior in that th 
small sharp maxima previously observed at an identical concentrati 
are now absent. At high drop-rates the characteristic decreas 
slope at the end of the wave remains. Cadmium thus appears | 
play the role of a suppresser. This is confirmed by the results give: 
in table 1. In the presence of 0.010 percent of gelatin the values of / 
for both ions agree very well with those obtained with single ions. 


TABLE 1.—Diffusion-current constants obtained from electrolysis of a solution 0.5 
millimolar in Ni and 0.297 millimolar in Cd 


[Supporting electrolyte was 0.1 N KC]] 


Cadmium Nickel 
Mercury j|__. _ ie =) 
flow 
| 0% gelatin | 0.01% gelatin 0% gelatin | 0.01% gelatin 


} 
mg/sec | 
1, 35 3. 47 | 3. 49 
1.79 3. 57 | 3. 46 
2. 40 3. 58 3. 46 
3. 40 4.09 3. 40 
5. 80 6. 24 3. 35 


The mixture of cadmium and thallium is of particular interes! 
because of the information it yields regarding the general resolving 
power of the method of analysis. In the polarograms for this mixtur 
without gelatin, the thallium wave is regular and well defined over thi 
entire drop-time range. The deposition of cadmium immediately 
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\J polarograms were taken at a sensitivity of 1/20, except those at a nominal drop-time of 1.1 seconds for 
“” percent and 0.002 percent of gelatin, which were taken at a sensitivity of 1/40. 
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after that of thallium completely prevents the appearance of the 
characteristic decreased slope at the end of the wave. In fact, this 
intermediate portion of the polarogram has an abnormally large 
positive slope. The shortness of this connecting segment makes it 
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Figure 12.—Polarograms for a cadmium-nickel and a cadmium-thallium miztur 
taken in the absence and presence of gelatin as a suppresser. | 
Polarograms are ordered from left to right with increasing flow of mercury. Those for the cadmiut ‘ 
<el mixture were taken at a galvanometer sensitivity, S=1/70, and those for the cadmium-thalliu 
ire were taken at S= 1/100, except the last one in each series for which S=1/200. The vaaes for Jp cal i 
from these curves are given in tables 1 and 2. 


very difficult to determine the real wave heights except at the highest 

drop-rate, where the segment is sensibly parallel to the residual 
current. The magnitude of the error for thallium as calculated 
from the Ilkovié equation (which assumes diffusion control) can be 
seen from table 2. : 
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TaBLE 2.—Diffusion-current constants oblained from electrolysis of a solution 
1.0 by 10-* M in thallium and 0.897 by 10-3 M in cadmium 


[Supporting electrolyte was 0.1 N KC1.] 


Cadmium Thallium 
Mercury 
flow 


0% gelatin | 0.01% gelatin!) 0% gelatin | 0.01% gelatin 


mg/sec 
1.35 
1.75 
2.40 
3. 40 
5. 80 


In this mixture the cadmium wave has the form observed when 
the ion is alone. However, the series is very abnormal in that the 
difusion current for this ion at high drop-rates is much too small. 
The reason for this is that the total current at the end of the cadmium 
wave is chiefly limited by the total normal diffusion current of both 
ions, and therefore thallium, which is characterized by a prominent 
maximum and is reduced first, assumes a disproportionately large 
share of the total current. The polarograms clearly show that the 
fractional reduction of the total current in the presence of 0.010 
percent of gelatin is very much less than that for thallium alone. 
In the presence of a suppresser the cadmium wave becomes normal, 
and the form of the thallium wave remains unchanged. The diffi- 
culty of determining the end of the thallium wave and the beginning 
of the cadmium wave can be seen from the fact that the values of 
I, for thallium are too low and those for cadmium too high. This 
reflects the difficulty of determining wave-height when the segment 
connecting two waves is very short and of abnormal slope. 

The three types of maxima encountered in this study are illustrated 
in figures 6, 9, and 11. The sharp maxima characteristic of univalent 
thallium and nickel (polarograms not shown) and the particularly 
simple departure from normal behavior exhibited by cadmium (fig. 6) 
can be taken as extreme types. The rounded maxima appearing at 
intermediate and high drop-rates in the case of bismuth (fig. 11) can 
be regarded as a transition type. The change caused in these polaro- 
grams by the addition of successive increments of gelatin is instruc- 
tive. Whereas relatively large concentrations of gelatin are required 
to suppress maxima of the cadmium type, only small amounts of 
suppresser are necessary to eliminate sharp or rounded peaks. The 
resulting polarographic waves, however, show departure from ideal 
behavior like those of cadmium and exhibit the same response to 
further additions of suppresser. It is also noted that for the lowest 
drop-time and at intermediate concentrations of gelatin the waves for 
thallium appear to be the resultant of a small peak maximum super- 
imposed upon a curve of the cadmium type. At low concentrations 
of reducible ion (see, for example, fig. 3) only the cadmium type of 
maximum is found. These observations suggest that sharp maxima 
may be primarily caused by adsorption and those of the more general 
and persistent cadmium type chiefly by streaming phenomena. 

The appearance of a reduced slope for that segment following the 
polarographic wave is not caused by a decrease in the value of the 

627528—45——_2 
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factor m?*t'” with increase in applied potential. The results of this 
study show this quantity is independent of gelatin concentrations as 
high as 0.012 percent and decreases by only 2.0 percent in the range 
of voltage from cadmium to nickel. This confirms the earlier work of 
Lingane and Kolthoff [3] who found a 2.2-percent decrease from —0.5 
to —1.1 volts. It is therefore incorrect to attribute the large decrease 
in slope so frequently observed to the operation of this factor. Greater 
deviations from parallelism of initial and final portions of a wave must 
be attributed to lack of diffusion control of the electrode process. 

Although there is conclusive evidence of a real deviation from the 
Iikovié equation for conditions under which it is generally assumed to 
be applicable, there might be the objection that the bulk of the dis. 
crepancies observed in this investigation are in reality caused by a 
peculiarity of the dropping electrode. However, if the mass of the 
drop delivered by a capillary (practically constant over a wide range 
of flow, gelatin content, and voltage) is considered a satisfactory 
capillary characteristic, a drop of 7.7 mg mass delivered by the electrode 
used in this study is well within the range 6 to 9 mg of several electrodes 
reported in the literature [11]. Moreover, experiments with a special 
capillary delivering drops of 40.7 mg show the same type of deviations 
It was observed, however, that conditions of drop-time, concentration 
of reducible ion, and gelatin content under which departure from 
ideality occur were noticeably different. 

An examination of all families of curves similar to figure 2 reveals 
the singular fact that for a fixed gelatin content the critical drop-time 
at which diffusion control ceases is (within experimental error) inde- 
pendent of the nature of the reducible ion, its concentration, and the 
supporting electrolyte. This fact suggests that for the simple metal 
cations diffusion control is primarily determined by characteristics of 
the dropping electrode and the nature and concentration of th 
suppressor. Table 3 summarizes the data obtained in this investiga- 
tion, and figure 13 is the corresponding logarithmic plot for th: 
individual ions. 

The results shown in figure 13 are in accord with the simple adsorp- 
tion theory for the suppression of maxima.’ From a consideration 


cible ions in the presence of gelatir 


TABLE 3.—Critical drop-times for various redu 


Nickel Cadmium | Thallium Bismuth 


Concentre- Critical) Con- Critical} Con- Critical; Con- | 
tion Gelatin! drop- | centra- Gelatin) drop- (| centra-|Gelatin) drop- | centra- 'Gelatin 
time tion | time tion | | time tion 


| Per- Milli- | Per- Milli- | Per- | | Milli- | Per- 

Millimolar | cent |Seconds\ molar cent \Seconds| molar | cent |Seconds| molar cent 
0. 002 2.90 |0. 00593 0.05 | 0.002 | 2.90 | 0.0311 | 0.002 
0.0284 .004 | 1.80 0.004 | 2.20 .004) 2.20] --| .004 
. 006 | os . 006 te | . 006 


002 3.70 , 002 | . .164 | .002 
2. 10 004; 23 . 004 2. | .004 
. 006 } . | . 006 


. 90 , 002; 3 1.04 . 002 
00 . 004 y . 004 
80 | eeee . 006 | . 006 
. . 002 | 
SO | ,004 
. 006 


| 


4 It is to be emphasized that the data in figure 13 do not provide information concerning the mechani 
for the formation of maxima. 
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of this theory and the kinetics of the adsorptive process, the concen- 
tration of suppresser required to give a minimum coverage of the 
drop should, for a given size of drop, be an exponential function of 
the drop-time. The dependence of this function on the size of drops 
has not been fully investigated. However, preliminary experiments 
with a capillary producing drops of approximately 41 mg indicate a 
significant dependence. his behavior is now under study. 
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FicureE 13.—Effect of gelatin on critical drop-times for several reducible ions. 


Each section of the figure is a logarithmic plot of the corresponding data of table 3. A given symbol in a 
t represents a fixed concentration. There is no necessary correspondence between the same symbol used 
0 different plots. 


IV. CONCLUSIONS 


One of the important conclusions to be drawn from this research is 
that gelatin (or some other suppresser) is necessary to obtain agree- 
ment with the Ilkovié equation over a considerable range of drop- 
times for capillaries giving drops in the range 6 to 9 mg. The need of 
a suppresser increases with the dilution of reducible ion and is in all 
cases most acute for drop-times less than 3 seconds. 

The concentration of gelatin required to effect diffusion control over 
a given drop-time interval is, for a given capillary, practically inde- 
pendent of the nature and concentration of the reducible ion. The 
effect of an excess of suppresser over that required to insure a con- 
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stant value of Jp is significant only at high dilutions of the reducible 
ion. 

A linear relation between the logarithm of the critical drop-time 
and logarithm of the gelatin concentration has been found. This re. 
lation has been shown to be independent of the nature of the reducib}; 
ion, its concentration, and the supporting electrolyte. 

The occurrence of “well-defined” waves so frequently mentioned 
in the literature has been found unreliable as evidence of diffusion 
control, unless a critical examination of the polarograms shows that 
the slopes of the segments preceding and following the wave do not 
differ by more than a few percent. Waves that conform to this cri- 
terion have always shown agreement with the Ilkovié equation 
whereas large deviations from this form have resulted in abnormal] 
diffusion currents. 

Evidence has been obtained that suggests that at least two factors 
are operative in causing maxima. One factor gives rise to sharp o: 
rounded peaks, the other to deviation from normal slope of that part 
of the polarogram following the polarographic wave. 

The interpretation of a polarographic wave is complicated by th 
presence of a maximum, the formation of which is a complex function 
of the characteristics of the capillary electrode, the concentration of 
the reducible ion, the presence of simultaneously reducible ions, and 
concentration of suppressive agents. Consequently, it is not advis- 
able to apply the diffusion-current constant and the Ilkovié equation 
to quantitative analyses unless it has been established that the ex- 
perimental conditions are within its range of validity. 

Since this paper was written an article by Lingane and Loveridg 
[13] has appeared, in which the dependence of the diffusion-current 
constant as a function of m?t'/ was studied for several capillaries 
Their results, which were for lead in a 1.0—N solution of potassiun 
chloride and in the presence of 0.01 percent of gelatin, when plotted 
like those in the present paper, give values for critical drop-times in 
good agreement with those shown in figure 13. The slight increase of 
diffusion-current constant with decrease in the value of m?“t'” (and 
hence with increasing drop-time) noted by these authors has also been 
observed and is commented upon in footnote 3 of the present paper. 
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COMPARATIVE LIQUID-JUNCTION POTENTIALS OF SOME 
pH BUFFER STANDARDS AND THE CALIBRATION OF 
pH METERS 


By George G. Manov, Nicholas J. DeLollis, and S. F. Acree 


ABSTRACT 


$y application of the equation pH=(H— E,,;— E;)/k to solutions whose pH 
values were known accurately, the sum of the potentials of the reference electrode 
and of the liquid-junction potential, F,.,+ E;, was obtained at 25° C by the method 
of cells with liquid junction for seven solutions suitable for standards of pH. The 
solutions used were 0.1018—m hydrochloric acid, pH 1.081; 0.01—m hydrochloric 
acid+0.09-m sodium chloride, pH 2.101; 0.05-m potassium acid phthalate 
+0.02—m potassium chloride, pH 3.989; 0.02—m potassium dihydrogen phos- 
phate +0.02—m disodium hydrogen phosphate +0.02—m sodium chloride, pH 6.863; 
0.02—m potassium phenolsulfonate+0.02—m potassium sodium phenolatesulfonate 
+-0.02—m sodium chloride, pH 8.795; 0.02—m boric acid+-0.02—m sodium borate 
added together as sodium tetraborate)+0.02—m sodium chloride, pH 9.155; 
and 0.01727—m calcium hydroxide +0.01819—m sodium chloride, pH 12.38. Silver- 
silver-chloride electrodes immersed in saturated potassium chloride solution were 
used rather than the calomel electrodes customarily employed. 

As E,,7 remains constant when the buffer is changed, values for the differences 
in the liquid-junction potentials of various buffers in contact with saturated 
potassium chloride solution were obtained from the data. These differences 
were then used to calibrate seven Type 015 and three “low-alkali error’ glass 
electrodes of commercial manufacture. The average agreement between the 
true pH of the buffer-chloride solution (determined from cells without liquid 

nections) and that read on various commercial pH meters when corrected for 

e difference in the liquid-junction potentials and the alkali error of the electrode 
was +0.01 pH unit. The data also furnish a critical test of the consistency of 
the pH values assigned to the various buffer solutions recommended by this 
Bureau for the calibration of the pH scale and for checking pH meters. 

Recommendations are made for checking pH meters. 


CONTENTS 


Introduction 
Experimental technic 
1. Method 
2. Electrodes 
3. Chemicals 
4. Equipment 
Measurements of electromotive force 
. Comparative liquid-junction potentials 
Comparison of available buffer standards 
Recommendations for the calibration of pH meters 
teferences 








116 Journal of Research of the National Bureau of Standards 
I. INTRODUCTION 


Recent publications from this Bureau have given provisional pH 
values for certain standard buffer solutions [1] ' and detailed informa- 
tion concerning phthalate [2], phosphate [3], phenolsulfonate [4], and 
borax buffers [5]. The precision of each standard, established by the 
method of cells without liquid junctions [6], is considered to be 0.002 
pH unit. 

In the establishment of the pH scale by electrometric methods, it is 
assumed that any solution in which the emf of a hydrogen electrode 
differs from that in a second solution by 59.14 mv at 25°C (2.3026RT7/F 
at other temperatures) has a pH value 1 unit higher or lower than 
that of the first. The establishment of the pH scale therefore involves 
the use of one or more standard buffer solutions of known pH to obtain 
reference points and the subdivision of the scale into desired fractions 
(e.g. 0.1) of the theoretical pH unit. 

It is evident that as many independent sets of calibrations of the 
scale of a pH meter can be made as there are available buffer stand- 
ards. As the number of reference materials increases, it becomes 
important to determine how well the corresponding calibrations agree 
Because the potential at the junction between the buffer solution and 
the electrolyte of the reference electrode (the so-called liquid-junction 
potential) enters into practically all measurements of pH, data on the 
comparative values for these liquid-junction potentials are needed 
One of the purposes of this paper is to present data on the comparative 
liquid-junction potentials at 25° C for seven standard buffers that 
range in pH from 1.081 to 12.38. The method is independent of the 
potential and of the type of reference electrode. 

As glass and hydrogen electrodes agree closely over the range of pH 
values from 1 to 9 [7], it is possible to use glass electrodes in pH meters, 
together with the values for the comparative liquid-junction poten- 
tials, to determine the accuracy with which the scale of the pH meter 
can be calibrated by the use of different buffers; that is, how closely the 
theoretical relation of 59.14 mv (at 25° C) for 1 pH unit is followed 
Conversely, if it is assumed that the standards are spaced correctly 
with relation to one another along the pH scale, it is possible to deter- 
mine accurately the performance of various types of glass electrodes 
immersed in different buffer solutions. 


II. EXPERIMENTAL TECHNIC 
1. METHOD 


The pH of a buffer or unknown solution can be obtained by meas- 
urement of the electromotive force of the cell 


—H,(g)|buffer or unknown solution|saturated KCljreference 
Ey E; Ever 
electrode + (1) 
and the equation 


pH=(E—E,.,— E,)/k, (2) 


in which £ is the measured emf, (E,+£,+E£,,.s) of the cell between 
the hydrogen and the reference electrodes, E, and £; in figure 1: 


! Figures in brackets indicate literature references at the end of this paper. 
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E..; is the potential of the reference electrode toward the hydrogen 
electrode in a solution of unit hydrogen ion activity when £; is zero. 
E, is the potential at the liquid junction between the buffer solution 
and the saturated potassium chloride, and k is the value for 
2 3026RT/F at the temperature of the measurements [8]. Conversely, 
if solutions are available whose pH values are known accurately, a 
measurement of & permits a simple calculation of E,,;+-E, for each 
buffer. Since the value for E,., does not change with the buffer, it 
is possible to obtain accurate values for the differences in E, for 
various buffers by a method that is independent of the nature of the 
reference electrode. 

If the data pertaining to two buffers of known pH are designated 
by the subscripts I and II, the difference in the liquid-junction 
potentials is given by the expression 


FE, — Ey, = (44— En) —k(pHi— pH), (3) 


directly in terms of the measured quantities. 

Independently, values for £,.;+E; can be obtained from the emf 
of two silver—silver-chloride electrodes, A, and E;, immersed, re- 
spectively, in solutions of the buffer (I) containing known concen- 
trations of chloride ions and in saturated potassium chloride (III). 
The emf of such a cell is given by the equation 


(Eyes t+ Ey— E,)=k log (@exy/Acr,,) + Ej. (4) 


If the ratios of the activities of the chloride ions in solutions I and 
[Il were known, eq 4 would yield directly values for £; for each buffer. 
It is not possible to obtain this ratio, however, and the data permit 
only the calculation of E,,,+ E;. 

The activity of the chloride ion in the two solutions can be expressed 
by the following equations: 


k log Qc, =k pH— (E,+ E,)+ E° (5) 
and 
k log ae1,,=E°—Enn (6) 


where E° is the normal potential of the silver-silver-chloride electrode, 
).22238 v at 25° C [9]. Equations 4, 5, and 6 can be combined to 
give 


Eyes + Ey=(E,+ Ey) + (Exes t+ E,— E,)—k pH. (7) 


Comparison of the sum of (£,+£,) and (£,.;+£,—£,) with the 
values for (#,+4,+£E,,.;) obtained by direct measurement indicate 
a consistency of 0.07 mv (0.001 pH unit) in the two sets of data for 
each buffer. This comparison is important because the emf between 
E, and E; involves electrodes reversible to the same ion, whereas the 
emf between EZ, and FE, and between FE, and &, involve electrodes 
reversible to different ions. 

[t should be emphasized that it is actually the difference in E, for 
various buffers in contact with the saturated potassium chloride 
which is used in the calibration of pH meters rather than the absolute 
values for EZ, for each buffer. The following example will make this 
pont clear. When a pH meter is calibrated with buffer I and the 
dial set to read the corresponding pH value, the setting of the “zero 
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control” knob, E,.,., compensates for the liquid-junction potentig! 
of the first buffer in contact with the reference electrolyte and fo; 
the departure,” AE yes, of the emf of the reference electrode from its 
normal value. The emf impressed upon the potentiometer system 
is then given correctly by the equation . 


Eon = Enep+ E+ Ee, + Ecerot+ SE nop 


gy |! 2 


where £,, is the potential of the glass (or hydrogen) electrode iy 
contact with buffer I. When buffer I is replaced by buffer II ang 
the setting F,,,, left unchanged, the observed emf is given by 


| Eves T Ein t | t ) t AE ves. 


The true difference in the pH of the two buffers is given by 


PH puey— PH truey= (Ens — Ee) [ke 


whereas the observed difference (as read on the meter) is 


PH 0001; — PH ove, = (E over: — Ene) /k 


(Ey — Ent Enn— Ex) /k. 


Since pHl,»,, and pH,,,., are equal, the correction to be applied to th 
I el | 


reading of the meter is 
PH trners = PH 0001, + (Ly, — Ey) /k. 


Because most meters use saturated potassium chloride solution as 
the reference electrolyte, the measurements reported here were co 
fined to this solution. 

2. ELECTRODES 

The hydrogen and the silver-silver-chloride electrodes were pre- 
pared by methods described previously [5]. Palladinized electrodes 
were used in solutions of potassium acid phthalate [2]. The silver- 
silver-chloride electrodes were of greater utiliy than the calome 
electrodes generally used, because they were less affected by air and 
showed considerably smaller thermal hysteresis. When the potassiun 
chloride solution was saturated with silver chloride, successi 
batches of electrodes agreed in 24 hours to within 0.02 my (0.00! 
pH unit) and remained stable within this limit for several months 
If the potassium chloride was not saturated with silver chloride, th 
emf drifted and became increasingly electropositive with tim: 
the silver chloride dissolved from the surface of the electrode. Drifts 
as high as 17 mv were observed. It is thus possible to introdu 
serious errors when microelectrodes are immersed directly in con- 
paratively large volumes of the solution. 

? This departure can be caused by exposure of the calomel to air [10] or by use of potassium « 


taining traces of bromides [11]. The zero control knob also compensates for the asymmetry | 
the glass electrode, but this is assumed to remain constant when the buffer is changed. 
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3. CHEMICALS 


Potassium chloride of reagent-grade quality stated to contain 
0.01 percent bromides was recrystallized twice from conductivity 
water. The recrystallized material and the corresponding mother 
liquor were tested by the following method at each stage of the 
purification process. Two electrodes that previously agreed to within 
),02 mv in a saturated solution of the purest potassium chloride were 
immersed, respectively, in (a) a solution of the crop of crystals and 
b) the corresponding mother liquor, both saturated at 25°C. Oxygen 
was excluded from the cell system, and the solutions were saturated 
with silver chloride. The compartments containing the two solutions 
were connected by a rubber tube with a pinchcock. Measurements 
were made of the emf of such a cell, the potential of which is due 
solely to the presence of different quantities of bromide ion in the 
two solutions. When tested as above, a solution saturated with 
potassium chloride containing 0.05-mole percent of potassium bromide 
and 99.95-mole percent of the purest potassium chloride was found 
to give an emf 1.46 mv higher than the purest potassium chloride 
alone. From measurements of this type it was concluded that each 
recrystallization removed approximately 90 percent of the bromides. 
The final product contained less than 0.0002-mole percent of bromides, 
which would produce a difference in potential of less than 0.01 mv. 
The hydrochloric acid used for the preparation of the silver—silver- 
chloride electrodes was also freed from bromides [5]. Although the 
values for £j,—;,, are independent of the emf and of the nature of 
the reference electrode, the above precautions were employed in 
order that the values for /,,;+ £; for each buffer would not include a 
systematic error because of the presence of potassium bromide or air 
in the reference electrolyte. 

The solutions were prepared in quantities of 10 liters from weighed 
amounts of conductivity water and specially purified chemicals. The 
conductivity water was prepared by C. G. Malmberg and had a specific 
conductance of 0.2107* mho/em® at 25° C. It was stored in an 
atmosphere of hydrogen until needed. The phosphates and phenol- 
sulfonate were purified by R. G. Bates [4]. Particular care was taken 
to prevent contamination with atmospheric carbon dioxide. The 
solutions were transferred when necessary by a vacuum-hydrogen 
technic [5]. The calcium hydroxide was prepared by carefully wash- 
ng pure calcium carbonate (low-alkali grade) with distilled water, 
after which the material was ignited at 1,000° C. The resulting oxide 
was allowed to react with an excess of water in a stoppered, paraffin- 
lined bottle and the caleitum hydroxide washed several times with 
conductivity water. The final solution was diluted slightly with water 
to avoid precipitation at higher temperatures, and the concentration 
determined by titration with standard hydrochloric acid. This pro- 
cedure for the preparation of solutions of pure calcium hydroxide was 
recommended by Flint and Wells [12]. 


4. EQUIPMENT 
The apparatus is shown in figure 1. Reservoir A was filled with a 


solution of potassium chloride and silver chloride saturated at a tem- 
perature slightly higher than 25° C and kept free from air by bubbling 
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hydrogen through the solution. Chamber B contained two silver 
silver-chloride electrodes, E, and E;, and was provided with a liquid. 
junction compartment, C, of the type used previously [13, 14] and g 
stopcock, D, leading to a waste jar, E. A second chamber, F, served 
to prevent accidental contamination of the main body of the buffe; 
solution in @ by possible diffusion of the saturated potassium chloric 
solution from the liquid-junction compartment C. Bubbler H was 
used to saturate the hydrogen gas with water vapor and was filled 
with the same solution used in F and G. The buffer solution under 
study was transferred to F and G without exposure to the air. Cel! 
G@ contained two silver-silver-chloride electrodes, E, and E,, and a 
double hydrogen electrode, E,, [14] immersed in the buffer solution 
The apparatus was then immersed in a water thermostat, the temper- 
ature of which was controlled at 25° C to within 0.01 degree. 


MYDROGE 
Pr rayages* 


Cai dCs 


ee 





Figure 1.—Apparatus used for the determination of E,.s+ E; for buffer-chloride 
solutions by the method of cells with liquid-junctions. 


A, Reservoir for the saturated potassium-chloride solution; B, chamber containing the two silver-silv« 
chloride electrodes E, and Es; C, the liquidjunction compartment with stopcock D ieading to waste )8! 
E; F, chamber to prevent contamination of the buffer chloride solution in G by acciaental mixing of t! 
saturated potassium chloride solution from B; H, the bubbler which saturates the hydrogen with water 
vapor before the gas enters the double hydrogen electrode chamber E;; E;and £;, silver-silver-chlorict 
electrodes immersed in the buffer-chloride solution. 


III. MEASUREMENTS OF ELECTROMOTIVE FORCE 


The stopcock between F and G@ was closed, and hydrogen gas was 
bubbled over the hydrogen electrodes in G. Measurements were 
made of the emf (£,+ £,) and of the agreement between the pairs of 
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electrodes, £3 and £, &; and £,, and of the double hydrogen electrode. 
The hydrogen electrodes agreed to within 0.01 mv in a number of 

preliminary trials, and for the balance of the measurements the two 
hydrogen electrodes were connected by a pool of mercury with a 
common terminal /, as shown in figure 1. Equilibrium was attained 
na few hours, after which the emf did not vary by more than 0.05 mv 
for several days. All values were corrected to l-atmosphere pressure 
of hydrogen gas. The values for (#,+,) reported here for hydro- 
chlorie acid, “hydroc hloric acid-sodium chloride, phthalate, phenol- 
sulfonate, and borax agreed within 0.05 mv with data published 
previously for these buffers; the emf for the phosphate buffer was 0.11 
my (0.0018 pH unit) higher than that subsequently found by Bates 
in an extension of his work [3], but is within the precision claimed for 
these measurements. The pH value for the phenolsulfonate buffer 
was recalculated by a somewhat more detailed method and is 0.006 
init higher then that previously reported [4a]. The pH for the hydro- 
chloric acid-sodium chloride solution was calculated from the data 
of Harned and Ehlers [9] and differs slightly from that used by Hitch- 
cock and Taylor [15]. 

After measurements of (4,+,) had been obtained for a buffer 
solution, cock D was opened, and the liquid junction between the 
potassium chloride and the buffer solution was made in C. The emf 
GC +E;+Ey7) and (E,.;+H,—H,) were measured at frequent 

itervals for various combined rates of flow and for various ratios of 
potassium chloride and buffer solution. The liquid junction for each 
buffer was quite sharp even at the highest rates of flow and the steadi- 
ness of the emf justified measurements to the nearest 0.01 mv. In 
table 1 is given a portion of the observed data and the calculated 
value for E,.¢+#; for a typical buffer, potassium acid phthalate. 


laBLE 1.—Electromotive force and other data at 25°C for the cell: 


Pt, H, | 0.05-m KH Phthalate, 0.02-m KCl || saturated KC] | AgCl, Ag+ 
Ag, AgCl (pH 3.989) 


f flow, drops 
per minute Time in min- 
| utes aftera | 
change in ee Sa ——— 
flow rate 


Electromotive force in international volts 


(Brest Ei— | (Ei +s) (Ey+ Es) (Es+ Ey) 


43) 


| 
—0. 13023 | 0. 56443 0. 00002 0. 00001 
. 13021 | 
13021 
13021 
. 13021 
. 13019 
13023 56442 00000 | —. 00003 
. 13021 
. 13024 
. 13026 
13026 56445 .00002 | +. 00001 
13026 


rage values Values of E,s+ E; 

Flow of buffer > flow of KC] | 0. 43432 | —0. 13023 | From E;:+E;+ Eve 0. 19839 
Flow of buffer ~ flow of KCl 43432 . 130: From Eyves+ Ej— Ex . 19827 
Flow of buffer < flow of KCl 43432 | , 


rand average r 0. 43432 , — 0. 19833 
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For comparison with other data in the literature involving calomel 
electrodes, 0.0455 v shouldibe added to E,.,+ E; for each of the buffers 
given in table 2. 


TaBLE 2.—Electromotive force and other data at 25° C for the cell: 


buffer, NaCl or KCl saturated KCilAgCl, Ag + 


—Pt, H,| 
Ag, AgCl 


Electromotive force in international volts Values of Eyes+ E; frox 


Buffer * 
(Ext+Ej+ | (EnectEi— | (E+ Ej+ 
Eves I 2 Eves 


\ FE + Ey 


b 0. 26292 b—(). 08850 
Hydrochloric acid ] 0. 35149 ©, 26206 ©, 08849 0. 19909 
1, 26318 1 O8829 


A verage 0. 26302 —(0). 08843 


b0. 32147 b—0O. 09065 
Hydrochloric acid+chlor- 2 0. 41209 e 3 5, © 09068 
ide. 1 32155 4 99061 


A verage 3214¢ 0. 09065 


013023 
Phthalate+chloride- : 0. 56443 |: 3432 13022 


Average 
Phosphate+chlorice 5 3 0. 73499 °. 6037 © 13118 197S¢ 
Average 


Phenolsulfonate+chlor- | 795 0. 84811 
ide | ‘ 726 é 13081 


Average 71725 13081 


12056 
Borax +chloride 9. 155 0. 86919 7 ©, 12054 
} 1 12954 


Average 73976 0. 12954 


92900 b—(), 13445 
Calcium hydroxide+ 2.3 5 92900 e. 13447 
chloride | 1, 92892 1. 13450 


A verage 92897 0. 13447 


* Composition of buffer solutions 

0.1018 m HC! 

0.01 m HC1+0.09 m NaCl 

0.05 m KHCsH,yO¢+0.02 m KCI 

0.02 m KH2P0,+0.02 m Na;3HPO,+0.02 m NaCl. 

0.02 m KHPs+0.02 m KNaPs+0.02 m NaCl; (Ps=paraphenolatesulfonate ion) 

0.02 m H;BO;+0.02 m NaBO;+0.02 m NaC! 

0.01727 m Ca(OH)2+0.01819 m NaCl 
> Flow of buffer > flow of KCl. 
e Flow of buffer ~ flow of KC] 
i Flow of buffer < flow of KCl. 


The data reported here for the emf and the pH of calcium hydrox- 
ide—sodium chloride solutions are believed to be the first obtained bj 
the method of cells without liquid junctions. The values 4.5A fo 
a, and 0.1 for 8 for this solution were estimated from data for other 
buffer-chloride mixtures. Reasonable departures therefrom (+0.° 
in a, and +0.05 in 8) are unlikely to give rise to uncertainties large? 
than 0.02 pH unit, a quantity which is negligible in the calibration 0! 


the majority of pH meters. A comparison with the available pH 
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data of other workers is given in table 3. The measurements of 

Wells [12a] and of Flint and Wells [12b, c] made at 30° C compare 

favorably with those of Lea and Bessey at 25° C [16] when correction 

is made for the change in pH with temperature for highly alkaline 

solutions. Measurements of the emf and the pH at 30° C by the 

method of cells without liquid junctions indicate that dpH/dT= 
0.035 pH unit per degree in this region. 


TaBLE 3.—pH values interpolated for 0.01727-m calcium hydroaide solution at 25°C 


Method pH at 25° C Investigator 


Conductivity 2.45 | Noyes and Eastman [17] 
Electromotive force, Hycalomel 2.36 | Wells [12a] ». 
do ‘ Flint and Wells [12b] ». 
do 36 | Flint and Wells [12c] . 


do 2. 4: Lea and Bessey [16]. 
do 2.45 | Fritsch [19] >. 
do : Kalousek, Jumper, and Tregoning [18]. 
Electromotive force, celis without 2. ; This investigation. 
liquid junctions 


» Corrected to 25° C by means of d pH/d T=—0.035 pH/degree at pH 12.4. 
mperature varied from 24° to 27° during the investigation. 


IV. COMPARATIVE LIQUID-JUNCTION POTENTIALS 


Values of E;,—;,, can be obtained for various pairs of buffers by 
~ gegen of eq 3 and 7 to the data in table 2. The entries in 
table 4 represent the Po values of (Z;,—£F;,,)/k in terms of pH 
units. In each case, Z,, designates the buffer listed at the top of 


columns 2 to 8 and £,,, is the corresponding buffer in column 1. The 
consistency of the data presented in tables 1 and 2 indicates that the 
values for (E;,—£;,,)/k in table 4 have the same degree of precision 
as claimed for the ‘pH measurements themselves, +0.002 unit. In 
general, the acid and the alkaline buffers have liquid-junction poten- 
tials which are, respectively, positive and negative to the saturated 
potassium chloride solution. 


Experimental values for (Ej,— Ej,,)/k at 25° C for various pairs of buffers 


in contact with saturated potassium chlorid solution 


- Hydro- > _— ~—— Calcium 
ph chloric ——_ ; we .. — Borax+ | hydrox- 
ehloric | ooiq4 ate phate ulfonate | chioride | ide+ 
acid chloride | CDloride | chloride |+chloride chloride 


Hydrochloric acid 032 013 | —0.020 | —0,033 014 | 03 


chloric acid+chloride | . 032 . 000 019 + —.001 018 | 01 
halate+chloride +. O13 —. 019 +. 000 - —. 020 .001 | .02 
phate+chloride . 0% —.012 . 007 +. —. 013 . 006 .02 
Phenolsulfonate+chloride . 03% +. 001 . 020 +. O1 +. 000 -, 019 .01 
Borax+chloride —. 018 001) — 


-.019 . 000 .02 
Calcium hydroxide+chloride . 0% +.01 +. 02 +. 01 02 +. 00 


It is beyond the scope of this paper to discuss the actual value of 
E,.; to be used in pH measurements performed with the hydrogen 
or the glass electrode. It is recommended instead that the meter be 
calibrated with two or more buffers at the start of the work. In 
most cases these extra measurements form but a small part of the total 
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effort expended. More important, however, such a calibration 
eliminates the errors otherwise present which may be caused by 
insufficient purification of the potassium chloride, improperly pre- 
pared reference electrodes including the effect of air on their 
potentials, hysteresis, and the possibility of undersaturation of the 
potassium chloride solution. For purposes of comparison, however, 
it is of interest to observe the values for F,,;+£;, which have been 
assigned to the saturated calomel electrode by various workers, 
These data are given in table 5. It should be noted that the valu 
243.4 +0.4 mv (+0.007 pH unit) for the use with average buffers 
(pH 4 to 9) is identical with that estimated by Hamer [22], 243.4, 
for buffers of pH 4 to 5. The variations in F,,;+£, given in table 
5 are not surprising in view of the possible sources of error. For 
example, 0.05-mole percent of bromides in the potassium chloride 
of the reference electrolyte increases the potential by 1.5 my; an 
additional increase of 1.5 mv is found if the electrolyte is saturated 
with air. 


TaBLe 5.—Comparison of the potential of the saturated calomel electrode at 25° ( 
obtained by various investigators 


Eves+ EF in millivolts 
Investigator ies ” 
Strong Average 
acids | buffers 


1925 Seatchard [20] : 245. 245. 4 | 
1928 Clark [21] 245 245.8 | 
1937 .| Hamer [22 : ; 245.7 | 243. 4 
tine Hitchcock and Taylor [15] 245. 244. 1 
1938 ....| MaclInnes, Belcher, and Shedlovsky [23]--.. ial 244. 6 
1044 .-| This research ®- al P 244. 6 b 243.4 








* Corrected for the difference of 45.5 mv in Extech, Evel pe-agcr’ 


> +0.4 mv. 


V. COMPARISON OF AVAILABLE BUFFER STANDARDS 


The buffer solutions in table 2 and the values of (E,,—£E,,,)/k in 
table 4 were used to calibrate a total of seven Type 015 and thre 
‘low alkali-error” glass electrodes. The apparent pH values at 25° C 
were measured in a constant-temperature room by means of the 
Beckman model G, the Coleman model 3C, the Hellige model 7040, 
and the Leeds & Northrup model 7661—Al commercial vacuum-tube 
pH meters. 

The results showed that the sum of the dial readings on the meter 
and the corresponding value of (£;,—£;,,)/k, together with the alkali- 
error correction of the glass electrode whenever this was significant, 
equalled the true pH of the buffer solution to within en average of 
+0.01 unit, regardless of the order in which the buffers were selected 
for use in the initial calibration of the meter. The averages of the 
data for the Type 015 and for the “low alkali-error’”’ electrodes in 
which a 0.05-m solution of potassium acid phthalate (NBS Standard 
Sample 84 b) was used for the initial calibration are given in table 6. 
The correction for the alkali error in 0.01727-m calcium hydroxide 
+-0.01819-m sodium chloride solution, 0.21 pH unit, was obtained by 
subtracting the corrected pH reading from the true pH. This 
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procedure is considered preferable to the use of the value 0.26 obtained 
by Fritsch [19], because the temperature was not carefully controlled 
during his experiments. 


VI. RECOMMENDATIONS FOR THE CALIBRATION OF 
pH METERS 


The buffer standards listed in table 6 can be used for calibrating 
pH meters at 25° C when appropriate corrections are made for 
differences in the liquid-junction potentials and the alkali errors of 
the glass electrode. The maximum difference found in the calibration 
of four commerical meters by the use of these buffers was 0.01 unit 
and is the limit of readability of most of these instruments. In the 
analysis of unknown solutions, it is recommended that a standardizing 
buffer be selected that most nearly approximates the unknown in pH 


and in composition. 


TaBLE 6.—Comparison of the calibrations of type 015 and “‘low-alkali error” glass 
electrodes in various buffer solutions at 25° C 


(0.05-m potassium acid phthalate (pH 4.008) was used for the initial calibration of the pH meters] 


Corrected 
} | pH read- | | Correc- value for the 
: | Known ing ob- . + yy ,| tion for | observed pH | pHaneen— 
Buffer « | pH value! served on (Eq~ Ei, ik alkali (sum of col- [PH ode (eorr.) 
meter | error umns 3, 4, 
| and 5) 


AVERAGE TYPE 015 GLASS ELECTRODE 


ig 
— 
oS 
= 
~ 


Phthalate . 008 | 
Hydrochloric acid 081 | 
Hydrochloric acid +chloride 2. 101 | 
Phthalate+chloride - 3. 989 | 
Phosphate-+chloride }. 863 | 
Phenolsulfonate+chloride 795 
Borax+chloride _. . 155 
Calcium hydroxide+chloride - - - 2. 38 


0. 00 | 0.00 | 4.01 (0. 00) 
—.01 .00 | 1.08 00 
+. 02 .00 | 2.10 00 
00 .00 | 3.99 00 

6 

s 

4 


Sze 


+.01 00 87 —.01 


+.02 .00 | 8.79 val +.01 


.00 02 | 16 .00 


NOSRON - 


+. 02 4.21 | 12.38 } 00 


I 
ena 


~ 


Average dif- | 0.00 
ferenee. | 


V-ALKALI ERROR” GLASS ELECTRODE 


AVERAGE “LOV 


| | 
00s | »(4.01)| 0.00 | 4.01 
O81 | 1.08 | ‘ 1.07 
101 | ‘ . 02 | .00 | 2.11 
9s . .00 | 3.99 
863 | 5 ° 6. 87 
795 | | .00 | 8.80 
155 5 9. 16 
38 35 | +. 02 | : 12. 37 


hthalate 
‘hlorie acid 
jrochloric acid+chloride 
thalate+chloride 
hate+chloride 
Phenolsulfonate+chloride 
Borax +chloride 
Calcium hydroxide+chloride 


roe rors 


° 
— 


| Average dif- 
| | ference. 
| 


* Composition of the buffer solutions 
0.05 m KHCsH,O,. 
1018 m HCl. 
1 m HC1+-0.09 m NaCl. 
05 m KHCsH,O,¢+0.02 m KCl. 
02m KH3P0,4+0.02 m Na3H PO;«+0. 02 m NaCl. 
(2m KH Ps+0.02 m KNaPs+0.02 m NaCl. 
02 m HyBO;+0.02 m NaBO;+0.02 m NaCl. 
01727 m Ca(OH)3+0.01819 m NaCl. 
* The dial of the pH meter was set to read this value, and the instrument was balanced with the “zero 
‘wntrol” knob. 
* Provisional value, +0.02 pH unit. 
* Obtained by difference. 
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The concentrations of the solutions that can be used for the calj. 
bration of pH meters are not limited to those given in table 6 but may 
include others that are either more dilute or more concentrated 
(references 1 to 5). Since the experimental values for (£,,—£,, 


in table 4 were determined for buffer-chloride solutions at 25° ¢ 
it is of interest to ascertain how much these quantities are change 
when chloride-free buffers and other temperatures are employed 
From the considerations given below, it can be said that the ch; ang 

are considerably less than 0.01 pH unit. 

Because (£,+/,—E,,7) and (F,,,;+£,—) for all buffers studie, 
were affected on the aver: age to less ‘ths an 0.1 mv by wide variations jp 
the rate of flow and in the ratios of the buffer and saturated potas. 
sium chloride solution forming the liquid junction, the differences 
(E;,—E),,)/k for the same pairs of buffers with and without chlorid 
are probably less than 0.002 pH unit. <A degree of turbulence in th 
formation of the junction sufficient to cause a mixing of only 4 parts 
of the saturated potassium chloride solution with 996 parts of th 
chloride-free buffer would raise the chloride content of the latter t 
0.02 m. Furthermore, at the concentrations used the difference 
the effective mobilities of the positive and the negative ions of 
chloride-free buffer is not altered appreciably when the solution is 
then made 0.02 m in sodium or potassium chloride. When the Hen 
derson equation [24] is used to calculate the differences in the poten- 
tials of the liquid junction between the reference electrolyte and th 
buffer solution with and without chloride, the major assumption 
made in the initial derivation of this equation are found to cane 
such a procedure should therefore be expected to give reasonably 
trustworthy values for these differences. The average differen 
calculated in this manner for the solutions listed in table 4, with and 
without chloride (hydrochloric acid and hydrochloric acid-sodium 
chloride excepted), is —0.003 pH unit.’ Since this average difference: 
is applied to each of the entries for £,.,+£; in table 2, no net chang 
results in the values for (;,—£;,,)/k in table 4 

Likewise, the temperature coefficients for the liquid-junctio 
potentials do not vary greatly. The calculated average liquid: 
junction potential for the buffers studied here is approximately 0.0! 
pH unit higher at 35° than at 25° C.* Again, the use of this averag 
difference causes no net change in the values in table 4. Until mor 
accurate data for the conductance of mixtures of ions become avai- 
able, it is recommended that the values in table 4 be used for chlori 
free buffers as well as for those that are 0.02 m in chloride and { 
temperatures up to or slightly beyond 35° C. The pH of these solv- 
tions varies with the concentration of the chloride ion and the temper- 
ature, and the appropriate pH values (references 1 to 5) must be used 
in addition to the corrections given in table 4. 


? The individual differences for phthalate, phosphate, phenolsulfonate, borax, and calcium hy 
buffers are, respectively, —0.004, —0.003, —0.003, —0.005, and —0.002 pH unit. The mobilitic 
and of Ps-- were estimated from data on other sulfonic acids. 

‘ The individual differences for the buffers listed in table 4 are, respectively, 0.002, 0.006, | 
0.009, 0.013, and 0.008 pH unit. 
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sSPECTROGRAPHIC DETERMINATION OF SODIUM, POTAS- 
SIUM, AND LITHIUM IN PORTLAND CEMENT WITH 
THE DIRECT-CURRENT CARBON ARC 


By Armin W. Helz* 


ABSTRACT 


Spectrographic methods, which are much less time-consuming than the usual 
chemical methods, are described for the determination of sodium, potassium, and 
im in portland cement. For the determination of sodium and potassium the 
unple is mixed with a mineral base containing silver for the internal standard. 
For lithium a standard mixture of graphite and strontium carbonate is added to 
the cement sample. In either case a fixed amount of the resultant mixture is 
placed in the crater of a graphite electrode, which is made the lower positive ele- 
t of a direct-current arc. The Na 8194.81 A, K 7698.98 A, and Li 3232.61 A 
s are used for the quantitative calculations. A discussion of the development. 
the method is given, followed by details of the adopted analytical procedure. 
spectrographie analyses of 41 cements for sodium and potassium are given 
rcomparison with the chemical values. 


CONTENTS 


Introduction__- 
Method of analysis for sodium and potassium 
1. Choice of method___ 
2. Choice of spectral region 
3. Development of the method___-- 
ll. Analytical procedure for sodium and potassium 
. Standard mixture- 
. Sample preparation and electrodes 
3. Making the spectrogram_ 
4. Photometry and analytical curves 
5. Procedure for high K,O content 
a Time required_--_- 
1\. Res ults of analysis for sodium and potassium 
\. Determination of lithium - 
Summe ary 
I. References 


I. INTRODUCTION 


"aaa analysis depends upon the fact that elements of a 
bstance can be made to radiate light, a spectrum of which has many 
ines characteristic of the elements in the source. The relative inten- 
ities of these lines for different elements bear a definite relation to 
he relative quantities, kind, and number of elements in the source, but 


“Research Associate at the National Bureau of Standards, representing the Portland Cement Associstion. 
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with the present state of knowledge of spectroscopy, this relation ; 
too complicated to predict. It has been found that by controlling th 

excitation used for the source and dividing analytical problems iny 
type cases, empirical data may be obtained that will give & rapii 
analytical ‘method. This method of analysis is ideally ‘Suited 
problems in which (1) a rapid determination is required, (2) the si 
of the sample is nec Beds small, (3) the rs srconn ol of th 
elements in question are small, or (4) the chemical method is difficy 
or uncertain. Many examples may be found in the literature on th 
subject where the accuracy of the results is as good as that of the bey 
chemical methods, and the time for the analysis i is a matter of minuts 
for the spectrographic method as compared to hours for the chemicg! 
method. 

The percentage error of spectrographic analysis tends to remaj 
constant regardless of the concentration, whereas, in chemical analysis, 
the absolute error tends to be constant. Thus, for high concentration 
chemical methods are preferable to spectrographic methods, whereas 
for low concentrations the reverse is true. In general, spectrographi 
analysis has proved advantageous for the dete rmination of ¢ oncentn- 
tions of elements below about 5 percent. 

The interest in spectrographic methods by the cement industry 
arises from the presence in portland cement of small quantities 0! 
metallic elements, the specific effects of which it is important to stud; 
Of particular significance at this time are the alkali elements, whic 
are under critical observation by various laboratories in connection 
with questions on the durability of concrete [1].!_ In the presen 
paper, the elements considered are potassium, sodium, and lithium 
The spectrographic method does not require the difficult separations 
that are involved in the chemical determinations of the alkali elements 
and it is much less time-consuming. This advantage of speed will 
even greater if the method is extended to include other minor elements 
such as magnesium, manganese, and titanium. 

No general method of spectrographic analysis has proved entire 
satisfactory. It is necessary to make a study of the type of s 
stance to be analyzed, because the intensities of the lines of the min 
elements depend not only on their concentrations but also in variow 
degrees upon the major elements and the physical and chemical stat 
of all the elements present. Portland cement, with its limited var 
ability of the major components, may be treated as a single type t! 
material for spectrographic analysis. Primary considerations i 
this application are the preparation of the sample and excitation 
the spectrum. Subsequent details of spectroscopy and photograph 
are readily decided according to established principles. The spect 
graphic procedure is outlined in detail for the benefit of those int! 
ceramic industry who are not familiar with this field. 


II. METHOD OF ANALYSIS FOR SODIUM AND 
POTASSIUM 
1. CHOICE OF METHOD 


In this study, the widely used internal-standard principle of Ger 
lach and the rotating-stepped-sector method of plate calibration ¥ | 
adopted. For discussions of these matters, papers by Duffendaci 


! Figures in brackets indicate the literature references at the end of the paper. 
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fe, and Smith [2], Twyman [3], Pierce and Nachtrieb [4], Serib- 
er [5], and Owens [6] may be consulted. But the problems involved 
n the treatment of the sample, the excitation of the spectrum, and 
he choice of the spectral region best suited to the present conditions 
ould be settled only by a planned program of experimental tests. 

Characteristics of portland cement, of immediate significance in its 
pectrographic analysis, are its low electrical conductivity, and its 
bility to sinter and melt at the temperature of the electric are. 
Procedures suitable for use with this type of material may be found 
n the literature relative to the analysis of minerals and inorganic 
‘hemicals [6 to 15], ceramic materials [16, 17], and ignited biological 
naterials [6, 7, 9, 18, 19, 24]. A common treatment is to place a 
mall quantity of the sample in the crater of a graphite electrode 
7, 16, 20). The sample is often mixed with graphite [11, 14, 17, 23] 
ya base material free from the elements sought [13, 19] and bearing 
vith it an internal standard ? [2, 4] and a spectroscopic buffer ? [24]. 

The sample mixture may be compressed into pellets and placed on a 
praphite electrode [7] or used as the electrode alone [21, 22]. A method 
esigned to use a solution may also be considered because the alkalies 
n portland cement are readily brought into solution by treatment with 
ilute hydrochloric acid. The solution may be injected into a flame 
9} or a small amount dried on a graphite electrode [8, 10, 12, 20] for 
we or spark excitation. 

The kind of excitation to use for best results is interrelated to a 
onsiderable extent with the condition of the sample. In the refer- 
ices given for the sample treatments discussed in the preceding para- 
raph, four kinds of excitation are noted: (1) flame, (2) direct-current 
ire, (3) alternating-current arc (stabilized by the use of a high-voltage, 
urrent-limited circuit), and (4) condensed spark. 

Several methods and a great many variations in procedure for the 
letermination of the alkalies in portland cement were explored, but 
or one reason or another were rejected in favor of the direct-current 
wre procedure. An alternating-current arc [8] with the sample mixed 
vith a powdered metal or graphite and compressed into pellets [21, 22 
dicated good possibilities except for requirements of rigidly controlled 
perating conditions. <A solution of the cement sample dried on a 
praphite electrode for use in an alternating-current are [8, 12] proved 
0 be impractical because of low intensity of the alkali lines. In the 
iethod finally adopted, a fixed quantity of the sample was placed in a 
upped graphite electrode and used in a direct-current arc. The 
preliminary survey indicated good possibilities and at the same time 
he simplest solution of the problem. 


2. CHOICE OF SPECTRAL REGION 


lt was desired to determine the alkali elements in concentrations 
etween about 0.05 percent and 1.5 percent, the latter value being 


inequal vaporization and excitation conditions from sample to sample and variations in 
properties from plate to plate, the intensity of a spectral line of the “unknown” element is 
itive to that of a line of an element either in the sample in a constant concentration or added 
This is the “internal standard,” which should possess the following properties: (1) The 
hould be in the same wavelength region as the analysis line, (2) the excitation characteristics 
hould be similar, (3) the two elements should have similar vaporization rates, (4) the inten- 
lines should be similar, and (5) a uniform dispersion of the reference material throughout the 
iid be obtained. For practical considerations the internal standard may not have all of the above 
variations in the results may be noted from sample to sample due to the presence of extra- 
sand variations in the excitation. These differences may be diminished by the addition of a 
ic buffer’’—an element (usually having a relatively high vapor pressure and a low ionization 
hich tends to maintain uniform excitation in the arc. 
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about the upper limit of concentration of any one alkali element , 
commercial cements. With the direct-current carbon arc the following 
lines have intensities convenient for quantitative work (expressed jy 
wavelength in angstroms):* Na, 3302.32, 8194.81 and K, 4044.}4 
7698.98. Lines in the 3,200 to 4,500 region have the advantage oye 
lines in the 7,300 to 8,600 region in that they may be photographed 
plates much more suitable for quantitative work. The K 4044.14 lin, 
however, occurs in a region of strong cyanogen bands, which interfer 
with measurements of intensity; but this background may be effe. 
tively reduced by a special exposure procedure or by the use of , 
nitrogen-free atmosphere [18, 19] for the arc. The 7,300 to 86% 
region has much less background interference, but here a difficulty 
is encountered with the K 7698.98 line, which exhibits self-reversal for 
the higher potassium concentrations.* (See fig. 1, B.) However 
certain refinements and improvements in the procedure made jt 
possible to use that line. Therefore the red region was chosen for the 
analytical procedure. This is illustrated in figure 1. Spectra D1 and 
D2 are in the 8,500 A to 7,600 A region. Each spectrum is divided 
into five strips of known relative intensity for photometric purpose 
explained later in the text. 


3. DEVELOPMENT OF THE METHOD 


The self-reversal of the potassium lines noted in the powder method 
was observed only for cements having the higher potassium concer- 
trations and was found to be subject to a partial control by an adjust- 
ment of the flow of air around the are and of the concentration ¢ 
the sample. Powdered graphite and alumina made the condition 


more acute, but silica improved the spectrum, forming, however, 
refractory melt that gave the are erratic burning characteristic 
A mixture of silica, alumina, and calcium carbonate, which wha 
mixed with the cement sample produces in the arc a fluid low-tempers- 
ture melt, effectively reduced the self-reversal of the potassium lines 
Under these conditions the potassium in the sample evaporated 
a steady rate for a period of about 1 minute. The addition of lithium 
carbonate to the mixture, for buffer action, further improved th 
behavior of the melt, but erratic results were still obtained, owug 
to the large volume of gas liberated in the arc. This sudden reless 
of carbon dioxide in some instances forced drops of the sample mu- 
ture to flow over the rim of the crater, but it was found that ths 
could be prevented by prefusion of the mixture or the use of a larger 
diameter electrode having a conical-drilled crater. Furthermore, tht 
use of this new electrode also decreased the rate of evaporation ( 
the alkali elements by requiring a longer time to reach the maximul 
temperature. 

Silver, introduced as the oxide, was selected for the internal standat 
because of the similarity in rate of evaporation of silver and t 
alkalies and the favorable supply of reference lines for sodium # 
potassium. This internal standard is not ideal but is a practic# 
compromise. For this reason, particular precision was required W 
such details in the procedure as the weight of the sample in the ele 
trode, the time of exposure, the arc current, the arc gap and its posit! 
in relation to the spectrograph. 

+ An angstrom is equal to 10-* cm. il 
4 “Self-reversal” is a decrease in the intensity of a spectral line caused by the absorption of radiation 


atoms in the central part of the source by atoms of the same kind in the outer part of the source 
appear very broad, blurred, and often completely transparent in the center. 
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III. ANALYTICAL PROCEDURE FOR SODIUM AND 
POTASSIUM 


1. STANDARD MIXTURE 


The standard mixture used in the preparation of the sample for 
analysis of potassium and sodium is made by thoroughly mixing 60 
percent (by weight) of SiO,, 14 percent of Al,O3, 20 percent of CaCO;, 
and 6 percent of Li,CO;. Highest purity powdered chemicals are 
used. The mixture is fused in a platinum boat by heating in an 
lectric furnace to 1,300° C, cooled and crushed to pass a No. 100 
sieve. An amount of finely divided silver oxide (alkali free) is added 
equal in weight to one-tenth of the weight of the fused and ground 
product just described. Thorough mixing by alternately grinding in 
an agate mortar and sifting through a No. 100 sieve completes the 
preparation of the standard mixture. Thirty grams of this mixture is 
a sufficient quantity for use in analyzing 100 cement samples. 


2. SAMPLE PREPARATION AND ELECTRODES 


For the determination of sodium and potassium, 0.300 g of the 
CaQ-Al,O3-SiO,-Li,0 mixture is mixed with 0.100 g of the portland 
cement sample. Small glass vials are convenient for weighing, mixing, 
and storing the samples. Time is saved if a number of vials are filled 
in advance with the required amount of the standard mixture. After 
he addition of the cement sample, mixing is accomplished with the use 
fan agate mortar and pestle and by sifting through a small No. 100 
ieve. Thirty milligrams of the mixture is placed in each of two 
raphite electrodes prepared for duplicate results. 
The electrodes are cut from second grade, spectroscopic carbons, 
hich are suitable for the determination of the alkalies in concentra- 
ions above 0.1 percent. The upper electrode consists of a 1-inch 
ength of \-inch-diameter stock with a tapered “pencil” point. The 
ower electrode similar to the design of Scribner and Corliss [25] con- 
ists of two parts, (1) a base and pedestal cut from \%- or \-inch-di- 
tmeter stock, which is clamped to the are stand and acts as a support 
or (2) the cap which carries the sample. With this design the sample- 
wearing area of the electrode obtains a more uniform and higher 
emperature than with the usual one-piece design. The base may be 
sed several times, but the cap and cathode are good for only a single 
xposure. The cap consists of a 15-mm length of %-inch-diameter- 
tock graphite with a hole in each end, one to fit the pedestal, which 
hay be cylindrical, but preferably the shape of a truncated cone, 
nd the other to receive the sample. The sample hole has conical sides 
larting at the top with a diameter of 8 mm and tapering down to 4-mm 
ameter ata depth of 4mm. For a depth of 3 mm beyond this point, 
ie hole continues with a cylindrical shape 4 mm in diameter. Special 
the tools make the cutting of these caps a simple and rapid operation. 


3. MAKING THE SPECTROGRAM 


A grating spectrograph having a dispersion of about 5 A per 
ulimeter was used. This high dispersion is incidental, however, and 
' hecessary for this procedure. Uniformity of illumination at the 
‘is aided by placing at the slit a convex lens [5], which produces an 
age of the arc on the grating of the spectrograph. Because the arc 
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in this work is an extended source with a large flame, and also becays; 
the vertical distribution of intensity of the internal standard lines jg 
different from that of the alkali lines, additional specification of th, 
illumination is necessary. Good illumination is obtained with the ap 
in such a position that the upper limit of its inverted image, or th 
rim of the lower electrode, coincides with the upper edge of tly 
grating. The image of the upper electrode then falls well within th 
grating area. This setting and the readjustment required during th 
burning of the are to keep the gap at a specified value may be facili. 
tated by projecting an image of the arc on a screen on the side of th 
arc opposite the spectrograph. A rotating stepped sector [4, 5] for 
producing a relative-intensity pattern, and a filter, complete th, 
optics outside the spectrograph. Light from the are passes in order 
filter, sector, lens, and slit. A Corning, Noviol B, G38 filter is needed 
for the sodium-potassium determination if a grating spectrograph is 
used. This filter cuts out the light of wavelength shorter than 4500 4 

The electrode holder, or are stand [20, 25], should be designed t 
provide rapid vertical adjustment of both electrodes, adjustment o! 
the upper electrode during the exposure to keep the gap constant 
and controlled ventilation. 

A 220-volt direct-current power supply with the current controll 
at 10 amperes is used. The lower sample-bearing electrode is mad 
positive [20] and the gap setat4 mm. An exposure time of 90 seconds 
is used, starting with the striking of the arc. The arc is convenient); 
started by the use of an auxiliary circuit, which momentarily induces 
a high-frequency, high-voltage, low-current spark across the gap [2 
However, the results are the same as when the simpler starting 
method is used of shorting with a carbon rod. Eastman IIR spec: 
trographic plates are used. The plates in this study were process 
in automatically rocked, temperature-controlled trays. 


4. PHOTOMETRY AND ANALYTICAL CURVES 


The determination of the relative intensities of the spectral line 
and their correlation with the concentration of the elements follov 
standard procedures, for which many references may be found (4, 5 
Measurements of the blackening of the lines are made with « » 
recording microphotometer. The lines measured are K 7698.98 A, Ag 
7687.78 A, Na 8194. 81 A, and Ag 8273.52 A. A graph of blackenm 
versus logarithm of the relative intensity is made for each region 
the spectrum and for each plate. The most suitable steps of the # 
kali and internal-standard lines are measured and calculated to rel 
tive intensity with the use of the plate-calibration curve. 

The alkali concentration is determined by reference to an anal 
tical curve giving the relation between concentration and relative ! 
tensity. This is determined by carrying out the entire spectrograpl 
procedure on cements of known alkali content. 


5. PROCEDURE FOR HIGH K,O CONTENT 


It has been stated that the direction of the arc flame, towards © 
away from the spectrograph, affects the self-reversal of the potassiu@ 
lines. This has been reduced by the choice of a standard mist 





Spe ctrographic Determination of Alkalies 35 


hut was not completely eliminated for high K,O concentrations. 
This difficulty may be surmounted in either of two ways. 

1) The ventilation of the are chamber may be controlled to keep 
the arc flame on the side of the upper electrode nearest the spectro- 
oraph. The draft required is very gentle and should not be strong 
enough to make the are burn erratically. If the are chamber is pro- 
nerlv ventilated, the self-reversal is not significant for K,O concentra- 
tions below 1.0 percent. If no draft is used, the limit of applicability 
of the procedure outlined in the preceding sections is lowered to about 
0.6 or 0.7 percent of KO. 

2) The procedure to be used here and also for K,O concentrations 
above 1.0 percent where controlled ventilation is used is followed 
exactly as outlined to the point of measuring the intensity of the lines. 
\ visual inspection is made to exclude from measurement those potas- 
sium lines obviously self-reversed. An analysis of the rejected samples 
for potassium is repeated with a more dilute mixture. Instead of 
taking 0.3 g of the standard mixture and 0.1 g of the cement sample 

sto 1), 0.450 g of the standard mixture and 0.050 g of the cement 

sample are taken (9 to 1). All other details of the procedure are the 
same, except of course the requirement of another analytical curve 
made with the more dilute mixture. 


6. TIME REQUIRED 


The procedure is rapid, requiring about 4 hours for the determina- 
tion of sodium and potassium in six samples as compared with 30 
hours for the chemical procedure. Increasing the number of samples 
and the extension of the spectrographic method to include other minor 
lements would increase the time advantage of the spectrographic 
method over chemical methods. The 4 hours for the spectrographic 
method does not include the time required to cut graphite electrodes 
and to prepare the standard mixture. These may be prepared in 
large quantities and kept on hand. 


IV. RESULTS OF ANALYSIS FOR SODIUM AND 
POTASSIUM 


The analytical, or working curve, when plotted in terms of the 
logarithm of the relative intensity (alkali line relative to the internal- 
standard line) and the logarithm of the concentration of the alkali, 
should be a straight line [4]. This was the case for KO, but the 
NO curve deviated from linearity owing to some sodium contami- 
nation either in the electrodes or in the standard mixture. In all the 
calculations the concentrations of the oxides of the alkali elements are 
used rather than those of the elements, to conform with the customary 
procedure for reporting analysis of cements.’ 

The reproducibility of the method is shown in table 1. Results 

given for 18 individual determinations of Na,O and K.O in the 

me sample. The deviations of values from the average illustrates 

tue errors due to fluctuations in the arc, inhomogeneity of the sample, 
LT 


i es for determining the analytical curves for Na;O and K,O were analyzed by J. J. Tregoning of 
ay and Silicate Products Division of this Bureau 


the ( 
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and inconsistencies in photometric measurements. The probable 
errors shown in the table (4.0 percent for duplicate determinations of 
Na,O and 5.8 percent for K,O) represent the precision of the method 
but the accuracy is demonstrated best by a comparison of the analy. 
ses of a large number of portland cements as is done in table 2, in which 
the average of the spectrographic results on each cement is compared 
with chemical results. The average percentage difference betwee, 
the chemical and spectrographic results is 0.06 for Na,O and 0.07 for 
K,O. This shows a much greater percentage error relative to the 
average concentration of Na,O and K,O than is indicated by the prob. 
able errors reported in table 1. This larger error cannot be ascribed 
alone to inaccuracies in the spectrographic method as it also includes 
inaccuracies in the chemical method. 


TABLE 1.—Reproducibility of the method for Na,O and K,O as indicated by 13 
determinations on the same sample of cement 


[Sample No, 1] 





Na;0 * | K,0» | 


% 

0.33 0. 
37 | : 

41 

. 38 


25 


. 37 
. 34 
4 
4 


4 
. 33 
. 35 
- 35 


. 36 
. 34 
40 
. 35 


. 38 
4 





Average... 0.36 














* Probable error for Na;O: 1 determination, 5.7% of the concentration; 2 determinations, 4.0. 
» Probable error for K3O: 1 determination, 8.2; 2 determinations, 5.8. 


Tables 3 and 4 give more detailed data for the first 21 cements 
listed in table 2. Two spectrographic “results” are given, a ‘‘result’ 
being the average of two determinations. The chemical values are 
the average of two or more determinations (with the exclusion of 
widely divergent values) and were obtained by five different labore- 
tories. The roman numerals at the head of columns of chemical 
results designate the laboratories, and the letters A and B designate 
the method (A, modified Berk-Roller [27] and B, the ASTM (23) 
method). The average difference between the maximum and minimum 
chemical results of different laboratories is 0.06 percent, which is the 
same as the difference between chemical and spectrographic results 


* Obtained through the courtesy of the Portland Cement Association; Berk-Roller method used. 
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jlso the average difference between the chemical results I A and 
| Bis greater than the average difference between — determina- 
ions for either I A or 1 B (not shown in the table). ‘These observa- 
‘ions emphasize that reproducibility alone is not a measure of the 


Taste 2.—Comparison of the results for Na;O and K2O (in percent) obtained by the 
spectrographic and the chemical (Berk-Roller) methods 


K;0 








Cement number | — 
Spectro- | Spectro- | « , 
graphic Chemical | Difference 
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curacy of either the spectrographic or the chemical method. The 
st estimate of the suitability of the spectrographic method for any 
wticular application is to be gained from a study of the data in 
able 2 comparing the spectrographic results with chemical results 
ota large number and large variety of cements. ; 
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T ABLI | Con parison o} ré sults for Na,O (perce nt obtained by the 8 pectrograr 
method and by different chemical lahoratories 


TABLI | ilis fo kK. percent obtained by, the 
ib j different che micat laborator ics 


0. OBS 


1 A,Modified Berk Roller met! 
B,ASTM method [® 
[he potassium line with tl! 1 sto 1) mixti for nents showed self-reversal 


were determined with a ¥ 
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V. DETERMINATION OF LITHIUM 


An extension of the spectrographic method to include all the alkali 
lements, that is, lithium, cesium, and rubidium was considered. 
‘ye determination of these elements by chemical methods is so diffi- 
it that the development of a spectrographic method would be useful 
n spite of the possibility of these elements occurring in cement only 
rely. It appeared possible that lithium might be found in appreci- 
‘ble amounts in portland cement, and therefore a spectrographic 
procedure for its determination was developed. The method was 
esigned to determine Li,O in amounts above 0.05 percent. 
In the selection of the lithium line to be employed for the determin- 
tion, it was noted that Li 3232.61 A had some advantages over 
possible lines in being of the desired intensity and located in a 
pectral region that would permit the simultaneous determination of 
ithium with a number of the other minor elements in cement. ‘These 
hdvantages seemed to outweigh the objection that, in its use, the 
multaneous determination of lithium with sodium and potassium 
sould not be possible. The Li 3232.61 A line was therefore selected 
r this study. 
Synthesized standards containing lithium were prepared in two 
ways. (1) Lithium carbonate was added to a CaQ-Al,O -SiQ, mixture 
give a product containing 2.89 percent of LiO. (2) Lithium 
rhonate was added to a low-total-alkali cement to give a product 
ontaining 2.00 percent of Li,O. These mixtures were calcined at 
1,300° C, then cooled and ground and a chemical analysis made to 
heck the lithium content. The above mixtures were then diluted 
th additional amounts of the base material to give a series of 
amples containing Li,O between 0.10 and 1.0 percent. 
The procedure used for the determination was similar to the one 
or sodium and potassium. A standard mixture consisting of 70 
reent of graphite and 30 percent of SrCO, (internal standard) was 
sed. The graphite may be obtained from the powdered product of 
ing and shaping electrodes. It is sieved to eliminate the size 
raction above No. 100 sieve. ‘Two-tenths gram of the graphite 
rCO; mixture is used with 0.2 g of the cement sample. Thirty 
ligrams of the resultant mixture is placed in a graphite electrode 
the are. A 40-second exposure on an Eastman spectrum analysis 
. | plate is used with an are current of 10 amp and are gap of 4 mm. 
ground- glass screen in the optical path near the are was used to 
nish the intensity of the spectrum in preference to the use of a 
nig — time. The Li 3232.61 A and Sr 3301.74 A 
) lines are measured for the analytical calculations. 
No lithium | was detected in 30 portland cements analyzed, indicating 
Li,O was not present in concentrations exceeding the limit of 
‘nsitivity for which the procedure was designed, about 0.05 percent. 
0 detect lower concentrations was considered of no interest at this 


There remains the possibility that lithium is present in some form 
&s easily detectable than that in the synthetic samples used for the 
tualytical curve, although the procedure employed (heating of the 
amples after the addition of lithium carbonate) diminishes the 
tobability of such an error. The reliability of the method was 





140 Journal of Research of the National Bureau of Standards 


checked by the use of materials of known lithium content, (1) a clay 
containing 0.21 percent of Li,O and (2) a spodumene containing 6,33 
percent of Li,O. The clay as obtained was analyzed spectrographically 
whereas the spodumene was added to several portland cements ¢ 
widely varying K,0 and Na,O contents in such amount as to produce 
mixtures containing 0.15 percent of Li,O. The alkali contents an 
spectrographic results for Li,O are given in table 5. 


TaBLe 5.—Comparison of results for Li,O (in percent) 


Spectro- 


Chemical graphic 


Sample * 


Clay 

Cement No. 3 
Cement No. 4 
Cement No. 5 
Cement No. 21 
Cement No. 43 





*® The chemical values for lithium were obtained as follows: Clay, chemical analyses; cements 3, 4 
and 21, the addition of spodumine; 43, the addition of lithium carbonate. 


It will be noted that in general the agreement between the spectr- 
graphic and chemical values for Li,O is good, the chief exceptions 
being noted with cements 5 and 21. This illustrates a condition 
sometimes experienced, of the effect of one element on the spectral 
intensity of another element. In cements 5 and 21, where the K,0 
content in the first case and the Na,O content in the second case an 
high, the Li,O result is too low. The chief reason for this was found 
to be an enhancement of the intensity of the strontium line (internal 
standard) by the extraneous alkali element, with only a minor effec! 
on the intensity of the lithium line. Problems of this character could 
be solved satisfactorily by further experimentation on technique and 
buffer agents, but the absence of quantities of Li,O as high as 0.0; 
percent in the cements examined made further perfection of tl: 
procedure seem unwarranted. 


VI. SUMMARY 


Procedures are described for the spectrographic determination 0! 
sodium, potassium, and lithium in portland cement by means of the 
direct-current carbon arc. Sodium and potassium are determined 
simultaneously by mixing the sample with a mineral base, by the use 
of a fixed amount of the mixture in an arc and recording the red regio! 
of the spectrum. For lithium, a mixture of the sample with a graphite 
base is used and the violet region of the spectrum recorded. Internal 
standards (silver and strontium) and spectrographic buffers are intro- 
duced with the standard mixtures. The Na 8194.81 A, K 7698.98 4, 
and Li 3232.61 A lines are used for the analytical calculations. Th 
procedures are rapid, requiring about 4 hours for the determinatil 
of sodium and potassium in six samples and about the same time 10! 
lithium. | 

The procedure was designed to minimize and control self-revers# 
vsamadl rate of evaporation of the alkali and reference elements 
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unequal distribution in the arc of intensity of the spectrum of the 
alkali and reference elements, and the effect of extraneous elements. 

An estimate of the reproducibility of the sodium and potassium 
determinations in terms of the probable error for the average of dupli- 
cate determinations is 4.0 percent of the Na,O concentration and 5.8 
percent of the K,0 concentration. The accuracy of the method can- 
not be expressed in such simple figures but is illustrated by a com- 
parison of chemical and spectrographic analyses of 41 portland cements 
where the spectrographic results were found to differ from the 
chemical results by an average value of 0.06 percent (of the sample). 
Chemical results obtained in five different laboratories are given in 
this comparison. 

The spectrographic examinations of 30 portland cements for lithium 
indicated that in no case was Li,O present in amounts in excess of 0.05 
percent. The method was checked by the analysis of several cements 
to which spodumene of known Li,O content had been added. 


The writer is indebted to the Optics Division of this Bureau for the 
facilities of the spectroscopic laboratory in which the experimental 
work of this research was performed, and especially to B. F. Scribner 
for the results of his preliminary study of this problem and for his 
continued interest and advice throughout the study reported in this 
paper. 
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HEATS, FREE ENERGIES, AND EQUILIBRIUM CONSTANTS 
OF SOME REACTIONS INVOLVING O,, H,, H,O, C, CO, 
CO,, AND CH,? 


By Donald D. Wagman, John E. Kilpatrick,’ William J. Taylor,? Kenneth 
S. Pitzer,® and Frederick D. Rossini 


ABSTRACT 


Values are presented for the following thermodynamic properties: The heat- 
content funetion, (#7 H,°)/T, the free-energy function, (F° — H,°)/T,the entropy, 
° the heat content, H° — H)°, and the heat capacity, C,°, for O» (gas) to 5,000 °K, 

gas) to 5,000 °K, H,O (gas) to 3,000 °K, N, (gas) to 5,000 °K, C (solid, graphite) 

» 1,500 °K, C (solid, diamond) to 1,200 °K, CO (gas) to 5,000°K, CQO, (gas) to 
3500 °K, and CH, (gas) to 1,500 °K; the standard entropy, S°, for H»O (liq) at 
15°C; the heat of formation, AHf°, the free energy of formation, AFf°, and the 

librium constant of formation, Af, from the elements, for H,O (liq) at 25°C, 
and for HO (gas), CO (gas), CO, (gas), and CH, (gas), to 1,500 °K; the increment 
n heat content, AH°, and the increment in free energy, AF®°, for the transition of 
raphite into diamond, to 1,200°K and to 20,000 atmospheres. 

From the foregoing, values were calculated for the increment in heat content, 
AH®°, the increment in free energy, AF°, and the equilibrium constant, K, for the 
following reactions, most of which are important in connection with the production 
of liquid hydrocarbon fuels from natural gas or coal and hydrogen: 


‘(solid, graphite) + CO,(gas) =2CO (gas). 
(solid, graphite) + H,O (gas) = CO(gas) + H,(gas). 
‘O(gas) + 1/20>(gas CO,(gas). 
)+ HO (gas) = CO,(gas) + H,(gas). 
1/20.(gas) = CO(gas) + 2H,(gas). 
+ CO.(gas) = 2CO(gas) + 2H2(gas). 
H.O(gas) = CO(gas) + 3H.2(gas). 
2H.O(gas) = CO,(gas) + 4H,.(gas). 


table of the fundamental constants used in the calculations is given. 
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volving O., H2, H,O, C, CO, CO:, and CH, 
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I. INTRODUCTION 


In connection with the work of the Thermochemical Laboratory and 
the American Petroleum Institute Research Project 44 at the National 
Bureau of Standards, values for thermodynamic properties of the 
simple substances oxygen, hydrogen, water, nitrogen, carbon, carbon 
monoxide, carbon dioxide, and methane have been recalculated. These 
data have been used to calculate the heats, free energies, and equili- 
brium constants of a number of reactions involved in the production of 
liquid hydrocarbon fuels from natural gas or coal and hydrogen. The 
results of these calculations are presented in this report. 


II. UNIT OF ENERGY, MOLECULAR WEIGHTS, CONSTANTS, 
ETC. 


The unit of energy used in this paper is the conventional thermo- 
chemical calorie, defined as follows [1]:* 


1 calorie=4.1833 international joules (NBS). (1 


The unit of pressure is the conventional standard atmosphere, 
defined as follows [2, 7]: 

1 atmosphere= 1,013,250 dynes/cm?’. (2) 

The constants in eq 1 and 2 are conventional constants, and are not 
subject to experimental error. 

The molecular weights were calculated from the values given in 
the 1941 table of International Atomic Weights [3]. 

The value of the PV product for an ideal gas (a real gas at zero 
pressure) at 0°C, (PV)3z8, which is equal to the product of the 
gas constant per mole and the absolute temperature of the ice point, 
(RT) wc, is taken from the value calculated by Cragoe [4]: 


(PV )er8= (RT) o¢= 2271.16 absolute joules/mole. 


The uncertainty in the foregoing relation is estimated to be +0.002 
percent. 

The relation between the absolute watt and the interrational watt 
(NBS) is taken from the value recommended by Crittenden and 
Curtis [5]: 


1 international watt (NBS)=1.00017 absolute watts. (4 


The uncertainty in the foregoing relation is estimated to be +0.005 
percent. 


‘ Figures in brackets indicate the literature references at the end of the paper. 
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The additional constants used in the present calculations are (1) 
the absolute temperature of the ice point, Tec; (2) the velocity of 
light, ec; (3) the Avogadro number (chemical seale) N; and (4) the 
Planck constant, A. The relation between the liter and the cubic 
centimeter has been utilized to convert several quantities to pressure- 
volume units, since they are often useful. The values recommended 
by Birge [6] for these constants have been used. 

The foregoing values of the constants lead to the consistent set of 
fundame ntal and derived constants that are given in table 1. The 
uncertainties of the constants have been reestimated in some cases by 
the present authors. 

The values of the constants in table 1 have been used to calculate 
the constants in the equations for the contributions of translation 
and of rotation of a rigid molecule to the several thermodynamic 
functions. The contribution of translation is given by the Sackur- 
Tetrode equation, and is exact for all molecules. The contribution of 
a rigid rotator is a useful approximation to the rotational contribution 
in actual molecules. See references [8, 9]. 


TABLE 1. rhe alues ade the constants 


Relation From 


1 liter = 1000.028 +0.004 cm! [6] 
1 atm=1,013,250 dyne/cm? Definition, [2] 
1 int. joule (N BS) = 1.00017 +0.00005 abs. joule 5] 
1 cal = 4.1833 int. joule (NBS) Definition, [1] 
= 4.18401 +0.00021 abs. joule d,c 
= 41.2930 +0.0021 cm? atm__.. | d,e, b 
=0.0412918 +0.0000021 liter atm___. d, c, b, a 
(PVE > Q=( RT hoc=2271.16 +0.04 abs. joule/mole..._. (4) 
= 22414.6 +0.4 cm! atm/mole h, b 
= 22.4140 + 0.0004 liter atm/mole h, b, ¢ 
c= 3 — +0.010 °K 
, ~~ Tec =8.31439 +0.00034 abs. joule/deg mole 
=8.31298 +0.00054 int. joule/deg mole 
= 1.98718 +0.00013 cal/deg mole -_ - 
= 82.0567 +0.0034 cm! atm/deg mole 
= ().0820544 +0.0000034 liter atm/deg mole 
N = (6.02283 +0.0022) «10"/mole 
k=(R/N) = (1.38048 +0.00050) x 10-"* erg/deg 
h = (6.6242 +0.0044) X 10-7" erg sec__- 
c= (2.99776 +0.00008)  10!° cm/sec_. 
Z= Nhc= 11.9600 +0.0036 abs. joule em/mole_. 
= 11.9580 +0.0036 int. joule cm/mole 
= 2.8585 +0.0009 cal cm/mole 
(27/R) = (he/k) = 1.43847 +0.00045 cm deg 


The equations in table 2 yield the values of the thermodynamic 
functions in units of calories, degrees, and moles. The logarithms 
are common logarithms, to the base 10. M is the molecular weight 
in grams, and TJ is the absolute temperature in degrees Kelvin. 
! is the value of the two equal moments of inertia, in g-cm’, of a 
linear molecule about axes perpendicular to the axis of the molecule. 
I,, I;, and J;, are the three principal moments of inertia, in g-cm ’, 
not necessarily unequal in value, of a nonlinear molecule. Separate 
equations are required for the rotational contributions of a diatomic 
or linear polyatomic molecule, and those of a nonlinear polyatomic 
molecule. The symmetry number, ¢, is the number of ways the 
molecule may be superimposed upon itself by rotation of the entire 
molecule. For linear molecules the value of ¢ is either 1 or 2; the 
equations for the rotational free energy and entropy of a linear mole- 
cule are given explicitly for these two cases. 
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TABLE 2.—Summary of equations, with constants, for calculating the values of the 
thermodynamic functions for translation (of all molecules) and rotation (of rigid 
j 
molecules) 


(H°—Hi)/T, C°, (F°—Hi)/T, and S° in cal/deg mole; Tin °K; M in g/mole; J, 1;, J:, and J; in g-en 
¢, a dimensionless integer; logarithms to the base 10.] 


TRANSLATION (OF ALL MOLECULES) 


H})/ T=C3=4.9680 
~H)/ T= —6.8635 log 11+-7.2820—11.4391 log T 
S° =6.8635 log !M—2.3141+11.4391 log T 
ROTATION (OF RIGID MOLECULES) 


I. Diatomic or linear polyatomic molecules 
H° — Hi)/ T= C® =. 1.9872 
1) o (Symmetry number) =1 
(} T = — 1.5757 log (7X 108") +2.7676—4.5757 log T 
symmetry number) =2 
"°—. 15) /T= 5757 log (7 108") +-4. 1450— 4.5757 log T 
S° = 4.5757 log (7 10%") — 2.1578+-4.5757 log 7’ 
II. Nonlinear polyatomic molecules. 
H°— Hi)/ T= C° =2.9808 
f° — }1§)/ T= — 2.2878 log (1; 1, Is 10"'7) +-4.5757 log #+3.0140—6.8635 log T’ 
S 2.2878 log (1; I, 13 10"'") — 4.5757 log ¢—0.0332+-6.8635 log T 


III. HEAT-CONTENT FUNCTION, FREE-ENERGY FUNC. 
TION, ENTROPY, HEAT CONTENT, AND HEAT CAPACITY 


1. METHOD OF CALCULATION 


Values of the thermodynamic properties of the simple molecules, 
O,, H,, H,O, N., CO, and COs, in the ideal gaseous state and to high 
temperatures, were calculated statistically from spectroscopic and 
other molecular data by Johnston, et. al. [10, 11, 12, 13], Giauque, 
et. al. [14, 15, 16}, Gordon [17, 18), and Kassel [19]. Simce such new 
spectroscopic data as have become available would not lead to values 
of the thermodynamic functions significantly different than thos 
obtained in the foregoing calculations, it was not considered necessary 
or justified to repeat these calculations in detail. The values required 
for the present investigation were obtained by converting the existing 
values to new ones based on the new values of the fundamental 
constants and the unit of energy. 

The important corrections required to be made arise from the 
following: (a) Change in the value of the gas constant R from 1.9869 
to 1.98718 cal/deg mole, the ratio being 1.00014; (b) change in the 
value of the constant of the Sackur-Tetrode equation for entropy of 
translation from —2.300 to —2.3141 cal/deg mole; (c) change in the 
value he/k from 1.43242 to 1.43847 cm deg, the ratio being 1.004224; 
(d) changes in the atomic weights of hydrogen and carbon from 1.0078 
to 1.0080 and from 12.00 to 12.010, respectively; (e) correction for 
rotational stretching in CO, and H,O. 

The foregoing corrections were handled in the following manner: 

(a) Each thermodynamic function was increased in value by tl 
factor 1.00014, the ratio of the new to the old value of the gas constant 

(b) A constant amount of +0.014 was subtracted from the value 
of the entropy for traaslation and added to the value of the free- 
energy function for translation, for each molecule. 

(c) The values previously calculated for the rotational and vibra- 
tional contributions to a given thermodynamic property were taket 
to be correct for a new, slightly higher temperature 7*, such that 
7T*/T=1.004224, the ratio of the new to the old value of he/k. This 
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procedure is possible because the previous calculations involved sum- 
mations with the partition function, in which hc/k appears as a 
multiplier of 1/7’, as in the expression 

Q= Sp e-te / (5) 

i 

The tabulated functions were interpolated from the temperatures 7* 
to even temperatures, after subtraction of the temperature-dependent 
term of the Sackur-Tetrode equation for the‘translational contribution 
to the entropy and the free-energy function. On a plot of any 
function f(7') against 7’, the temperature for which the calculation 
was originally made, the corrected value of f(7) for the temperature 
T is the value of the function at (1—1/1.004224) T=0.004217' degrees 
below 7. In practice, the slope of the function at 7’ was calculated 
by averaging Af(7)/AT for the intervals above and below 7. The 
correction to f(T) due to the change in he/k is then —0.00421T 
df(T)/dT. The correction varies slowly with the temperature. 

d) The entropy and the free-energy function were increased and 
decreased, respectively, by 3/2RIn(1+4M/M), where M is the mole- 
cular weight and A.M the change in the molecular weight for the given 
substance. Since AM/M is very small, this expression reduces to 
3/2RAM/M or 2.984M/M. 

(e) A small additive correction for rotational stretching, as given by 
Wilson [20], was made to each of the thermodynamic functions of H,O 
and CQ. 

The complete equations used for the conversion are as follows, with 
the old values primed: 


H° — H>3)/T=1.00014{ |(H° —H>)/TY 
0.00421 7Td|(H°—H>)/T\'/dT}- (6) 


Hs) /T=1.00014{[(F° — H>)/T)’ —0.00421[Td| (F°— H>)/TY /dT 
+-5/2R’|4+-0.014}—2.98AM/M. (7) 


S°=1.00014{ (S°)’—0.00421[Td(S°)’/dT—5/2R’|—0.014} 
+2.98AM/M. (8) 


OC: =1.00014{ (C)’ — 0.00421 Td(C3)’ /dT}- (9) 

In the case of methane, the statistical calculations made by one of 
the authors [21] were extended to include other temperatures. 

In the case of carbon (graphite) and carbon (diamond), the existing 
calorimetric data were subjected to graphical analysis to yield selected 
“best”? values for the heat-content function, free-energy function, 
entropy, heat content, and heat capacity, to 1,500° K for graphite and 
to 1,200° K for diamond. 


2. DATA USED IN THE CALCULATIONS 


The data used in the present calculations, according to the pro- 
cedure described in the preceding section, were from the following 
sources: 

Q,(gas).—Statistical calculations from Johnston and Walker [10, 11], 
who gave values for the free-energy function, the entropy, the heat 
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content, and the heat capacity, to 5,000° K. See also Gordon and 
Barnes [22]. 

H, (gas).—Statistical calculations from Davis and Johnston [12), 
who gave values for the free-energy function, entropy, heat content, 
and heat capacity to 5,000° K, and from Giauque [14], who gave values 
for the free-energy function. See also Gordon and Barnes [23] and 
Libby [24]. 

H.O (gas).—Statistical calculations from Gordon [17, 18], who 
gave values for the free-energy function, entropy, and heat capacity, 
to 3,000° K. See also Libby [24]. 

N. (gas).—Statistical calculations from Johnston and Davis [13], 
who gave values for the entropy, heat content, and heat capacity to 
5,000° K, and from Giauque and Clayton |15], who gave values for 
the free-energy function to 5,000° K, 

CO (gas).—Statistical calculations from Johnston and Davis [13], 
who gave values for the entropy, heat content, and heat capacity, to 
5,000° K, and from Clayton and Giauque |16], who gave values for the 
free-energy function to 5,000° K. See also Gordon and Barnes [22], 
and Kassel [25]. 

CO, (gas).—Statistical calculations from Kassel [19], who gave 
values for the free-energy function, heat content, and heat capacity, 
to 3,500° K. See also Gordon [26] and Gordon and Barnes [23]. 

CH, (gas).—Statistical calculations from Pitzer [21], who gave 
values for the free-energy function, heat content, entropy, and heat 
capacity to 1,500° K. 

C (solid, graphite).—Experimental data on heat capacity and heat 
content from Nernst [27] from 30° to 85° K, from Worthing [28] (in- 
direct measurements) from 1,170° to 2,400° K, from Magnus [29] 
from 370° to 1,200° K, from Schlapfer and Debrunner [30] from 
370° to 1,400° K, and from Jacobs and Parks [31] from 87° to 300° K. 
See also Magnus [32]. 

C (solid, diamond).—Experimental data on heat capacity and heat 
content from Weber [33, 34] from 283° to 600° K, from Magnus and 
Hodler [35] from 670° to 1,180° K, and from Pitzer [36] from 70° 
to 288° K. See also Bettendorf and Wiillner [37], Nernst and Linde- 
mann [38], and Robertson, Fox, and Martin [39]. 


3. RESULTS FOR 0O,, H;, H,0, N:, C, CO, CO, AND CH, 


The resulting values of the thermodynamic properties for O, (gas), 
H, (gas), H,O (gas), N, (gas), C (solid, graphite), C (solid, diamond), 
CO (gas), CO, (gas), and CH, (gas) are presented in tables 3, 4. 5, 6, 
and 7, which give, respectively, values of the heat-content function, 
(H°-—H,)/T, the free-energy function, (F°-H;)/7T; the entropy, S°, the 
heat content, H°—H,, and the heat capacity, C,°. 

A value for the entropy of liquid water at 25° C may be obtained 
from the value for gaseous water at 25° C and the standard entropy 
of vaporization of water, as given by eq 20: 


H,O (liq.) ; Soos-16= 16.716 +0.016 cal/deg mole. (10) 
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IV. HEAT OF FORMATION, FREE ENERGY OF FORMATION, 
AND THE EQUILIBRIUM CONSTANT OF FORMATION 


1. METHOD OF CALCULATION 


Values of AH, the increment in heat content for the given reaction, 
with all the reactants and products in their appropriate standard 
states at 0° K, were calculated by combining values of AH. 16, the 
increment in heat content for the given reaction at 25° C, with values 
Of Hoys..—H, for the reactants and products, by means of the 
relation 


AH, AH x5 16 — A(z w— 1), 1] 


where the last term on the right is the sum for the products less the 
sum for the reactants of Hox .4—H>, the heat content at 298.16° K 
(25° C) less that at 0° K. 

Values of AH° for the given reaction at other temperatures were 
then calculated by means of the relation 


AH°= AH, + A(H°— A). (12 


Values of AF°, the increment in free energy for the given reaction, 
with all the reactants and products in their appropriate standard 
states at the given temperature, were calculated by combining values 
of AH; for the given reaction with appropriate values of (F°—H))/T, 
the free-energy function, by means of the relation: 


AF°/T=AH,/T+ Al((F°—H,)/T), (13) 


where the last term on the right is the sum for the products less the 
sum for the reactants of (7°—H))/T, the free-energy function at the 
given temperature. 

Values of the equilibrium constant, and its logarithm, for the given 
reaction, were calculated from the relation [40] 


AF°=—RT In K. (14) 
2. DATA USED IN THE CALCULATIONS 


The following data were used in calculating the heat and free 
energy of formation of H,O(liq) at 25° C and of the formation of 
H,O (gas), CO (gas), CO, (gas), and CH, (gas) at various tempera- 


tures to 1,500° K: 
For the heat of formation of liquid water, 


H, (gas) + 1/2 O,(gas) = H,O(liq), (15 


the value previously reported from this laboratory [41] was corrected 
to the present molecular weight of water, yielding for the formation 
of liquid water according to eq 15: 





- 
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AHs.16= —68,317.4 +9.6 cal/mole. (16) 

On the basis of the present molecular weight, the values previously 

reported for the standard heat, free energy, and entropy of vapori- 

zation of water at 25° C [41], including a small change resulting 

from an improved value for the vapor pressure of water at 25° C [42], 
become 

H,0 (liq) = H,O (gas). (17) 

AH 98.16=10,519.5 +3.1 cal/mole. (18) 

AF 98.16 = 2,054.8 +1.0 cal/mole. (19) 


ASoos-16= 28.390 +0.012 cal/deg mole. (20) 


A new value for the heat of formation of carbon dioxide was 
recently reported [43]: 


C(solid, graphite) + O,(gas) =CO,(gas). (21) 
AH oys.16= —94,051.8 +10.8 cal/mole. (22) 


The value previously reported [41] for the heat of combustion of 
carbon monoxide to carbon dioxide remains unchanged: 


CO(gas) + 1/2 O2(gas) =CO,(gas). 
AH ox8.16= — 67,636.1 +28.7 cal/mole. 


Combination of eq 21, 22, 23, and 24 yields a new value for the 
heat of formation of carbon monoxide: 


C(solid, graphite) + 1/2 O,.(gas) =CO(gas). (25) 
AH %.16= —26,415.7 + 30.7 cal/mole. (26) 


The value for the heat of formation of methane is taken from 
reference [44]: 


gas) =CH, (gas). (27 


C(solid, graphite) --2H.(5 
AH 98.16 — 17,889 +75 cal/mole. (28) 

3. RESULTS FOR H,0, CO,CO;,AND CH, 
The resulting values of the thermodynamic properties for the 


formation of H,O(gas), CO(gas), CO.(gas), and CH,(gas), from the 
elements, O,(gas), H.(gas), and C(solid, graphite), are presented in 
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table 8, which gives values of the heat of formation, AHf°, the free 
energy of formation, AFf°, the equilibrium constant of formation, Kj 
and the logarithm of the equilibrium constant of formation, log,,Af 
all to 1,500° K. 

For the formation of liquid water from its elements at 25° C, the 
values of the standard heat, free energy, and entropy of formation 


H, (gas)+ 1/2 O, (gas)=H,0 (liq). (29 
AH oes. 16= —68,317.4 +9.6 cal/mole. 
AF o08.1= —56,689.9 +11.2 cal/mole. (31) 
ASSos 16= —38.997 +£0.019 cal/mole. (32 


V. HEAT AND FREE ENERGY OF TRANSITION OF 
GRAPHITE INTO DIAMOND 


The value for the heat of transition of graphite into diamond at 
25° C is taken from reference [43]: 


C (solid, graphite)=C (solid, diamond) 
AA os.16= 453.2 +20.3 cal/mole. 
From the foregoing value and the values for diamond and graphite 
given in tables 3, 4, 5, and 6, values were calculated for the change 
in heat content and free energy for the transition of graphite into 


diamond at 1 atmosphere, as presented in table 9. 


TABLE 9.— Values of the change in heat content and free energy for the transition of 
graphite into diamond, at a pressure of 1 atmosphere 











Temperature in degrees Kelvin 
C (solid, graphite) =C = 


298. 16| 300 | 400 | 500 | «00 | 700 | s00 | 00 | 1, 0001, 100 1, © 





(solid, diamond) = 
0 


- | srl 


Heat of reaction, AH®, in | | | | 
cal/mole.... | 576.6) 453.2) 452.1) 399.0) 358 = 310} 299) 291) 
| 


| } 
| | | | | | . 
in cal/mole -. 576. 6) 685.0) 686.7) 774.0) 873) 978)1, <= 1, 201) 1, 314) 1,428) 1, 543) 1, 663 


Free energy of reaction, AF°, 


The expressions for the thermal expansion and compressibility of 
graphite and diamond have been reported previously [45]. Assuming 
that (1/V) (0V/0P),r for graphite and diamond are constants inde- 
pendent of temperature and pressure, these values may be combined 
with the thermodynamic relation (OAF/0P);=AV to calculate the 
change in free energy for the transition of graphite to diamond at 
temperatures up to 1,200° K and pressures up to 20,000 atmospheres: 


AF?=AF°+A(F?— F°). (35 





the equilibrium consta 


TABLE 8.— Values of the heat of formation, a free energy of formation, AFf°, logarithm o a, the equilibrium constant of formation, logioKf, « 
in nl H,0 (gas), CO — te e Gon) and CH, (gas), from the elements in their standard states 














‘lemperature in degrees Kelvin 








Reaction 800 900 1,000 | 








Heat of formation, AHf°, in kilocalories per mole 








entiakiomeaa - —_———- 
—57.7979 | —57.8022| —58.042| -58.276| —58. 499 | —58.709 | 58.902} 59.080 | 


) 
solid, graphite)+1/203 (gas)=CO (gas)_. 
solid, graphite) +O3 (gas) =COs (gas) _. 


Fiy(cas)+1/203 (gas) = H30 (gas 

( 

C _ 
C (solid, graphite) +2H; (gas) =CH, (gas)... Rte ee 


Hy gas)-+1/203 (gas)= HO (gas) 

C (solid, graphite)+1/203 (gas) =CO (gas) 
solid, graphite) +O (gas) = COs (gas) 

C (solid, graphite)+2H; (gas) = C H, (gas) 


H,(gas)+1/203 (gas) =H 0 (gas) 

C (solid, graphite )+1/202 (gas)=CO (gas) -. 
C (solid, graphite)+O, (gas) = C Ox (gas) 

C (solid, graphite)+2Hs (gas) =C Hy, (gas) -. 


— 26. 4157 
—94. J518 


— 54. 6351 
— 32. 8079 
—94. 2598 
—12. 140 








40. 04695) 
24. 04790] 
69. 09145 
8. 8985 | 


— 26. 4131 
—¥4. 0520 
—17. 903 


— 54. 6152 
— 32. 8464 
—94. 2603 
—12. 104 


39. 78683 
23. Y2845 
68. 66801) 
8. 8177 


— 26. 317 
—94. 069 
— 18. 629 


— 53. 516 
—35. 007 
— 94. 325 
— 10. 048 


29. 23972 
19. 12672 
51. 53648 
5. 4899 


— 26. 295 
— 94. 091 
—19. 302 


— 407 
— 04. 167 
— 2. 401 


— 28. 511 
—04. 215 
— 20. 823 





— 04. 268 
—21. 166 


— 26. 635 | 


—26. 768 
—94. 318 
| —21.43 | 


Free energy of formation, ary’, in kilocalories per mole 





— 52. 358 
— 37. 184 
—94. 392 

—7. 840 








~48. 643 | 
—43. 677 | 
—94. 539 | 


— 49.912 
— 41. 526 
— 04. 497 
—3. 05 


—47. 349 | 
— 45. 816 
—94. 578 


—0. 55 +2. 01 








| 
— 46. 036 | 
—47. 942 
—94. 610 
4.61 


Logarithm of the equilibrium constant of formation, logwK/ 





22. 88551 


16. 25283 | 


41. 25820 
3. 4268 


| 
18. 63228 
14. 33621 
34. 40107 
2.0001 
! 


13. 22846 | 
11. 93193 | 
25. 82664 | 
0. 1494 | 


Equilibrium constant of formation, Kf 


11. 49776} 
11. 12559) 
22. 96647) 
—0. 4881 


15. 58315 
12. 96479 
29. 50309 
0. 9526 | 


10. oni04| 
10. 47772) 
20. 67675 
—1. 0075 





|, 


8. 88300) 
9. 94448) 
18. 80256) 
1, 4345 


| 


7. 89864 
9. 49826 
17. 23998) 
— 1. 7936 | 
| 


.114X10 @ 6. 12110 % |1. 73710 7.918X107 |1.158x10 
3.150K10* (1.31110 
(1. 73810 '7 | 8. 2! 

}1. 608X107? (7. 9: 


63810 8 
R00 X10 
347 X10 8 
677 X 1077 


. 683X102 4. 28810 18 (3.83010 (1.943X1u 8 \5. 14610 "! [1.151 X10 1 \7. 

.117X10 * |8. 48110 % |1.339X1u 1 11. 790X10 1 |2. 16910 4 |9. 22110 2 |8.549X10 " 11. 33510! 13. 004X110 10 18. 

- 234X10 © |4. 65610 % [3.43910 8! |1.812K10 4 |2. 51810 # a mexte® 6. 709 X10 % |9. 25710 33 14. 75110 * /6. 

: eee § 16.572X10*% |3.090X10 5 - 672X103 |100 8. 966 ' 4ll 0. 3250 9. 829X10-3 /3. 
| 


His(ens +1/204 (gas) = H;0 (gas) .... _— — 
( id, graphite) +1/203 (gas)=CO (gas) _. : 
soli i, graphite)+O;3 (gas) =COs (gas) 
solid, graphite)+2H, (gas) =CH, (gas) 


52™ 1 
32x 10 





TasBLE 11.—Values of the logarithm of the equilibrium constant, logwK, and the equilibrium constant, K, for some reactions involving Or, Ha, | H,0O, Clgraphis ), ¢ 














id, graphite) +CO2(gas) =2CO (gas). . 
lid, graphite)+H3O (gas)=CO (gas) +He (gas) 
) gas) +1/203 (gas) = COs (gas) 
) (gas) +H2O (gas) =CO, (gas) +Hs (gas) 


H, (gas) +1/202 (gas)=CO (gas)+2Hs; (gas). ------- 


Hy, (gas) +C Os (gas)=2CO (gas)+2Hy (gas) 
Hy, (gas)+H,O (gas)=CO (gas)+3H:2 (gas) 
CH, (gas) +2H30 (gas) =C Os (gas) +4Hs (gas) - -.-- 


id, graphite) +CO; (gas)=2CO (gas 

i, graphite)+H 30 (gas)=CO (ea) sis (gas) 
O (gas)+1/203 (gas) =C Ox; (gas) 
© (gas) +H»O (gas) =COy (gas) +H: (gas) 
Hy (gas)-+1/20% (gas) =CO (gas)+2He (gas) _- 
Hy, (gas)+ + (gas)=2CO (gas)+2Hs; (gas) .. 
Hy, £28) +s (gas)=CO (gas)+3H: (gas) -_..---- 
ar gas) +2H30 (gas) = COs (gas)-+4Hs (gas)----- 





Temperature in degrees Kelvin 








Logarithm of the equilibrium constant, logw:K 





—20. 99575 
— 15. 99896 
+45. 04367 
+ 4.99679 
+15. 14943 
— 29. 89351 
— 24. 89745) 


—19. 90067 


—20. s10s6| 


—15, 85786 
44. 73973 
4. 95303 

15. 11042 
—29. 62931 | 
—24. 67628 
—19. 72325) 


— 13. 28281 
—10. 11277 
32. 40960 
3. 17004 
13. 62866 
— 18. 78203 
—15. 61144) 
—12. — 


—8. 75242 
—6. 65217 
25. 00541 
2. 10025 
12. 82615 
— 12. 17882 
— 10. 05889 
—7. 93898 


— 5, 72851 
— 4. 29593 
20. 06491 
1. 43258 
12. 33631 
—7. 72859 
—6. 29601 
— 4. 86343 


—3. 57358 
—2. 61852 
16. 53817 
0. 95506 
12. 01203 
—4. 52613 
—3. 57077 
— 2. 61571) 


—1. 96284 
—1. 35664 
13. 89475 
0. 60620 
11. 78248 
— 2. 11227 
—1. 50580 
—0. 89960 


—0. 71538 
—. 37226 
11. 84070 

0. 34312 
11. 61365 

—0. 22705 
+. 11607 
+. 40919 


40, 27 
+. 


0. 137 
11. 
+1.% 





Equili 


brium constant, K 


. T5658 
59959 
74189 
15699 
29168 

3. 54995 
39314) 
. 23615 





1.012 10-2! 
1.002 10-16 
1.106 X 1045 
9.926 104 
1.41110" 
1.278 10-* 
1.266 X 10-35 
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1.546 X 10-2! 
1.387 X 10-1¢ 
5.49210 
8.975 X 104 
1.290 10!5 
2.348 X 10-% 
2.107 X 10-25 
1.891 X10-% 





5.214 10-" 
7.713X10-" 
2.568 X 10% 
1.479108 
4.253 X10" 
1.652 10-'* 
2.447 X 10-6 
3.623 X10-" 





1.768 X10-* 
2.228 X10-7 
1.013 10%5 
1.260 X 102 

6.710 X10" 
6.625 X 10-8 
8.732X10-"! 
1.151 X10-* 





1.868 X10-* 
5.059 X 10-§ 
1.161 X 10% 
27.08 
2.169 X 1012 
1.868 X 10-® 
5.058 X 10-7 
1.369X10-$ 


2.669 X 10-4 
2.407 x10" 
3.453 X 1016 
9.017 

1.02810" 


1.098 X10-? 
4.399 X 10-2 
7.848 X 108 
4.038 

6.060 X10!! 





2.978 X10 


7.722X10-% 


2.687 X10- | 3.12010 


2.423 X 10-* 





0.1260 


0.1926 

. 4244 
6.930 X 10!" 
2.204 =| 
4.108 X10)! 
0.5929 
1.306 
2.879 





1.900 
2.609 
1.582 10" 
1.374 
3.056 X10! 
19.32 
26.56 











11.58 


7.210 X 108 
0.9444 


2.392X10"' 


3.316X1v? 
3.13310? 


2.959 X lu? | 





39.77 


5.519X 107 
0.6066 





ibrium constant of formation, logiKf, and the equilibrium constant of formation, 


m the elements in their standard states 





Temperature in degrees Kelvin 
























































700 800 : 90 6} «(1,000 «| ~—s, 100 | 1,200 | 1,300 | 1,400 | 1,500 
formation, AHP in kilocalories per renal 

> TED -— : _ oe : 

—58. 709 —58. 902 | —59.030 | —59.239 |) —59. 384 | —59. 511 —59. 623 | —59. 724 | —59. 811 
— 26. 407 — 26. 511 —26.635 | 26.768; —26.909 | —27.056 —27. 212 —27. 376 | —27. 545 
—O4. 167 | —94. 215 | —%4. 268 | —-04.318 | —94.364 —4. 410 —04. 456 —94. 505 —94. 555 
~meor} -men| 2108) —26 | -20 | -a7 | -2e | -200 | > —22. 06 

| | | 

: of formation, arr in kilocalories per mole 

—49. 912 —48. 643 — 47. 349 | — 46. 036 —44. 710 | —43. 370 | —42.017 —40. 661 —39. 296 
—41. 526 — 43. 677 —45. 816 | — 47.942 | — fw. 053 — 52. 153 — 54. 235 — 56. 308 — 58. 370 
—04. 497 —94. 539 —94. 578 | —4. 610 — 94. 637 —04. 661 | —94. 677 —94. 690 | — 94. 707 

—3. 05 —0. 55 | +2. 01 4.61 7. 22 9. 85 12.50 15.14 17. 80 
the inna constant of formation, logwKf 

15. 58315) 13. 22846 11. 49776) 10. 06104! 8. 88300 7. eoses| 7. ona 6. 34747 5. 72542 

12. 96479 11. 93193 11. 12559 10. 47773 9 94448 9. 49826) 9. 11762 8. 78999 8. 50449 

29. 50309, 25. 82604 22. 0647 20. 67675 18. 80256 17. 23998 15. 91654 14. 78159 13. 79863 

0. 9526 0. 1494 —0. 4881 —1. 0075 —1. 4345 | — 1. 7936 —2. 1006 — 2. 3638 —2. 5927 
uilibrium constant of formation, Kf 
| | a | | 

83010 8 |1.943X1u ™ |3. 14610 | /1.151X10 © (7.638108 (7.918107 |1.158xK107 |2.226xK10% | 5.314105 
22110 "2 (8. 549X110"! |1.335K10 " |3.004X10 '° (8. 800K10*® |3.150K10" |1.311xK10* \6. 166X10* | 3.195108 
185X<10 * 6.70910 35 |9. 257 X10 3 | 4.75110 ® |6. 34710 '§ |1. 73810 "? |8. 25210 '§ |6.040K10 * | 6.29010 4 
966 3. 677X10-? ! 608 < 10-2 - 932X10-9 |4.327K10- | 2.554x10-% 


[lau ( 3250 ie 829 10-2 


a 

















t, K, for some reactions involving 02, H;, H.O, C(graphite), CO, CO,, and CH, 








Temperature in degrees Kelvin 















































700 | 800 | 900 | 1,000 | 1,100 | 1,200 | 1,300 | 1,400 1,500 
| 
ithm of the equilibrium constant, logwK 
] 7 l 
—3. 57358: —1. 96284 —0. 71538! +0. 27865 1. 08638) 1. 75658 2. 31863 2. 79835 3. 21033 
—2. 61852 — 1. 35664 —. 37226) +. 41655 1. 06373 1. 59959: 2. 05385 2. 44243 2. 77906 
16. 53817 13. 89475 11. 84070 10. 19920 &. 85796 7. 74189 6. 79887 5. 99165 5. 29440 
0. 95506 0. 60620 0. 34312 0. 13790} —0.02484) —0.15699] —0. 26478} —0.35592/} —0. 43127 
12. 01203 11. 78248 11. 61365 11. 48514 11. 37881 11. 29165 11. 21808 11. 15376 11. 09715 
—4. 52613 —2. 11227) —0. 22705 +1. 28504 2. 52066 3. 54995 4. 41921 5. 16211 5. 80289 
—3. 57077 —1. 50580) +. 11607 1. 42428 2. 49602 3. 39314 4. 15460) 4. 80634 5. 37177 
—2. 61571 —0. — +. 40919 1. 56218 2. 47118 3. 23615 3. — 4. 45042 4. 94050 
Equilibrium constant, K 
2.669 X 10-4 | 1.098X10-? | 0.1926 1.900 12.20 57.09 2.083102 | 6.28610? 1.623108 
2.407 X10-% | 4.399 10-2 . 4244 2.609 11.58 39.77 1.135 x 10? 2.770 10? 6.013 X 108 
3.453 X10" | 7.84810" | 6.93010"! 1.58210" 7. — 5.519107 | 6.293108 | 9.810105 1.970 X 105 
O17 4.038 2.204 1.374 0.9444 0.6966 0.5435 0.4406 0.3704 
1.02810"? | 6.06010" | 4.10810" | 3.05610"'| 2.39210"| 1.95710"| 1.65210" | 1.42510" | 1.25110" 
2.978 X 10-5 | 7.722K10-' | 0.5929 19.32 3.316X lv? 3.548 X10" | 2.626 10+ 452X105 6.352K105 
2.687 X 10- | 3.120107? | 1.306 26.465 3.133 X 102 2.473 X10 | 1.428104 6.402K 104 2.354 X 105 
2.423 X 10-* | 0.1260 | 2.879 36.49 2.959 X lu? 1.723X10° | 7.750K 10° 2.821 X10* | 8.720104 
' 
627528—45 (Face p. 156) 


1? 
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The first term on the right is given in table 9, and the second term on 
the right is the increment in free energy from 1 atmosphere to P 
atmospheres for diamond less that of graphite, at a given temperature. 
Values of (F?—F°) and A(¥?— F°) for various pressures and temper- 
atures are given by the relations: 


For C (solid, graphite)— 


(Fe— F°) =[0.12843 +-7.617 X 10-77'+- 4.848 X 10-77} (36) 
((P—1)—1.52 X 10-*(P*— 1) + 1.54 X 10-"(P*— 1)] cal/mole, 


For C (solid, diamond)— 
(F?— F°) =[0.082777 — 1.505 X 10-7 7'+-5.631 X 10-”T? 


43.6224 X 10-®73]{(P—1) —8.0 10-8(P?—1) 
+4.0X 10-(P®—1)] cal/mole, 





Cleotid, graphite) = C (solid, diamond) 


g 


i 
: 
; 


@ 
So 
i=] 


CAL/MOLE 
n 
°o 
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4 
> 
°o 
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Pe gait 8 074 
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i i i 
200 400 600 800 1000 1200 








TEMPERATURE, °K 
Ficure 1.—Plot of the free energy change for the transition of graphite into diamond. 


The scale of ordinates gives the increment in free energy, AF°, in cal/mole, for the transition, C (solid, 
graphite)=C (solid, diamond). ‘The scale of abscissas gives the absolute temperature in degrees Kelvin. 
As indicated, the three curves are for pressures of 1, 10,000, and 20,000 atmospheres. Values below zero 
indicate, within the limits of uncertainty, the region in which diamond is thermodynamically more stable 
than graphite. See text. 
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For the transition 


C (solid, graphite)=C (solid, diamond) 
A(F?— F°) [0.045653 +- 9.122 « 10-77 — 7.83 X10-"'T? 
3.6224 « 10-8 7*}[ P— 1]+ [1.886 & 107-°+ 1.17 X10-"T 
+- 6.92 * 107-7? —2.9 « 10-”7*)][ P?— 1]—[1.974 « 10-" (39 
1.174 10°"°7+-7.44  10-"7?— 1.45 x 10-*7")[ P® — 1]cal/mole. 
In eq 36, 37, and 39, P is in atmospheres, and 7’ is in degrees Kelvin, 
Values of the change in free energy for the transition of graphite 


into diamond at various temperatures and pressures are plotted in 
figure 1. 


VI. HEATS, FREE ENERGIES, AND EQUILIBRIUM CON. 
STANTS OF SOME REACTIONS INVOLVING O,, H,, H,0, C, 
Co, CO,, AND CH, 


From the data given in the preceding sections, values were calculated 
for the increment in heat content, the increment in free energy, the 
logarithm of the equilibrium constant, and the equilibrium constant, 
to 1,500° K, for the following reactions: 

C(solid, graphite +CO,(gas) =2CO(gas). (40 
C(solid, graphite) +H,O(gas) = CO(gas) + H2(gas). 
CO(gas) + 1/20,(gas) =CO,(gas). 
CO(gas) + H,O(gas) = CO,(gas) + H2(gas). 
CH, (gas) + 1/20,(gas) =CO(gas) + 2H,(gas). 
CH, (gas) +-CO,(gas) = 2CO(gas) +-2H,(gas). 
CH, (gas) + H,O(gas) = CO(gas) + 3H,(gas). 
CH, (gas) + 2H,O(gas) =CO,(gas) +4H,(gas). 

Values are presented for the above reactions, to 1,500° K, in table 10 

for the increment in heat content, AH®°, and the increment in free 


energy, AF°, and in table 11 (facing p. 156) for the logarithm of the 
equilibrium constant, log, and for the equilibrium constant, A. 


me inwvoale:na (> 
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In figure 2 are plotted the values of the logarithm of the equilibriyy 
constant for these reactions and for reaction (27). 





cO- HO . co, - H, 





CH,+2H,0 . co, *4Ho—~ 


CH,- CO, = 2CcO *2h, 





a | i. 





1 L di I a. i 1 


1 
200 400 600 800 1000 1200 1400 
TEMPERATURE, °K 





Figure 2.—Plot of the logarithm@of the equilibrium constant for eight react 


involving O,, Hz, H,O, C (graphite), CO, CO., and CH,. 


The scale of ordinates gives the logarithm (to the base 10) of the equilibrium constant, log,X, for t! 
reaction. The scale of abscissas gives the temperature in degrees Kelvin. The curves apply t 
following reactions (see text) 


(a) C(solid, graphite) +CO>2(gas) = 2C O(gas). 

(b) C(solid, graphite)+H»O(gas) = CO(gas)+H2(gas 
(c) CO(gas)+1/202(gas) = C O2(gas). 

(d) CO(gas)+H20(gas) = CO2(gas) + H2(gas). 

(e) CH«(gas)+CO2(gas) = 2C O(gas) +2H2(gas). 

(f) CH«(gas)+H20 (gas) = CO(gas)+3H2(gas). 

(g) CH«(gas)+2H20 (gas) = C O2(gas) +4H>(gas) 

(h) C(solid, graphite) +2H2(gas) = C Hi(gas) 
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HEATS OF ISOMERIZATION OF THE 18 OCTANES 
By Edward J. Prosen and Frederick D. Rossini 


ABSTRACT 


The heats of isomerization of all of the 18 octanes were determined by measure- 
ment of the ratios of their heats of combustion, in the liquid state for 17 and 
the solid state for 1, using the procedure previously described for the hexanes 
and the heptanes. 

Values of the heats of isomerization are reported for the condensed (liquid or 


solid) state at 25° C and for the gaseous state at 25° C and 0° K. 
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I. INTRODUCTION 


In continuation of the program of determining the heats of forma- 
tion and of isomerization of the paraffin hydrocarbons [1, 2, 3, 4]', 
calorimetric measurements have been made which yield values for 
the differences in the heats of combustion of all of the 18 octanes at 
25° C. These values for the differences in the heats of combustion, 
which, for the same state, are also the heats of isomerization at 25° C, 
have been combined with heats of vaporization (or sublimation for 
the one solid octane), measured by Osborne and Ginnings [5], to give 
values for the heats of isomerization in the gaseous state at 25° C. 
Combination of these with values for the decrease in heat content 
f the various isomers in the gaseous state from 25° C to 0° K, caleu- 
ated by Pitzer [6], gives values for the energies of isomerization for 
he ideal gaseous state at 0° K. 


*igures in brackets indicate the literature references at the end of this paper. References to the earlier 
I 1 in [1] 
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II. UNIT OF ENERGY, MOLECULAR WEIGHTS, ETc. 


The unit of energy upon which the values reported in this pape 
are based is the international joule. For conversion to the conyep. 
tional thermochemical calorie, the following relation is used [7]: 


4.1833 international joules (NBS)=1 calorie. 


The molecular weight of carbon dioxide, the mass of which was usw 
to determine the amount of reaction, was taken as 44.010 from thy 
1941 table of atomic weights [8]. 

The uncertainties assigned to the various quantities dealt with iy 
this paper were derived, where possible, by a method previously 
described [9]. 

Definitions of all symbols used are given in a previous paper (2) 


III. METHOD AND APPARATUS 


The same method and apparatus were used in the present investi 
gation as in the work on the hexanes [2] and heptanes [3]. 


IV. EXPERIMENTAL PROCEDURES 
1. CHEMICAL PROCEDURE 


The description of the octanes investigated is summarized in table! 
which gives the name of the compound, the normal state (liquid 
solid) at 25° C, the source of the material, and, where known, tly 
amount of impurity. For 13 of the octanes, determination of th 
amount of impurity was made by A. R. Glasgow, Jr., of the American 
Petroleum Institute Research Project 6 at this Bureau, from time 
temperature freezing and melting curves [19]. For the others, est- 
mates of the amount of impurity were made by the authors, consider- 
ing the method of purification and the purity of related isomers. Thi 
octanes used were so prepared and purified that the residual impurities 
must have consisted largely of isomeric octanes. Further, the impur- 
ties would be largely only those isomers having boiling points neat 
that of the given isomer. The values given in table 2 indicate th 
tolerances with regard to such impurities, and also the actual est: 
mated amounts of impurity. From these data, it appears that th 
amounts of impurity in the 18 octanes were not sufficient to affe 
their heats of combustion significantly, with the possible exception 0 
2,3-dimethylhexane. . 

The procedures for sealing the liquid hydrocarbons in glass ampoules 
purifying the oxygen used for combustion, igniting the hydrocarbot 
examining the products of combustion, and determining the amout! 
of reaction have been described previously [2]. 
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. 
- 


TC TABLE 1.—WSource and purity of the octanes investigated 
> Paper Amount of 
‘OnVen- impurity 


|: 


Mole fraction 


n-Octane (1 Liquid 0.0080-+0.0010 
ane (II do “ d .0035+0.0005 
9-Methylheptane do | <.03 
hylheptank do J <.03 
as Use 4-Methylheptane. - . do | y. <. 015 
Ol ho -Ethylhexane. . do <. 02) 
n td 9-Dimethylhexane do i <.015 
n ethylhexane do A, <. 057) 
-Dimethylhexane do .) 02 
with t '; Dimethylhexane do | A, Or 
ViOUSIT 
Isty methylhexane do , <. 005 
imethylhexane do my ® b(<. 02) 
. ) 2-Methyl]-3-ethylpentane do J <. 01 
pe r i“ Methyl-3-ethylpentane do | A, <. 006 


lrimethylpentane do | : <. 005 


imethylpentane do | <—. 002 
imethylpentane . do A, <.01 
lrimethylpentane do A <. 005 
Tetramethylbutane --. Solid <. 02) 





Che letters in this column have the following significance: 

Prepared by the American Petroleum Institute Hydrocarbon Research Project at the Ohio State 
University. These materials were made available by Cecil E. Bourd, Supervisor of the Project, and 
George Calingaert, Chairman of the Subcommittee on Physical Properties. 

Further purified by the American Petroleum Institute Research Project 6 at the National Bureau of 

andards by means of azeotropic distillation in columns of high efficiency (approximately 100 theo- 
al plates) at a high reflux ratio (100 to 1); 

purified from relatively pure 2,2,4-trimethylpentane prepared commercially (Rohm and Hass Company, 

Philadelphia, Pa.) by azeotropic distillation (as under B), by the American Petroleum Institute 


table | Research Project 6; 

, D, purified by crystallization by A. R. Glasgow, Jr., of the American Petroleum Institute Research Project 
{Ula Or 6, from material prepared by the Chemical Research Laboratory of the Ethyl Corporation, Detroit, 
vn. th Mich. This material was made available by George Calingaert, Director of Chemical Research of 

oe the Ethyl Corporation. 

ol th Estimated; see text. 

erica . ; 

» time For the combustion of the solid 2,2,3,3-tetramethylbutane, glass 
3. esti- ampoules could not be used, and, since the sample was required to be 
naides inclosed because of its volatility, gelatin capsules were used for the 
. Th purpose. All the capsules were from one lot of No. 000 gelatin cap- 
vurities sules.2. These were dried by heating in an oven at 100° C, cooled in 
mpur- a desiccator, and weighed in small closed bottles out of contact with 
S neal the air of the room. A heating time of 6 to 8 hours usually served 
te th to bring the capsules to constant weight. For each combustion experi- 
al esti- ment involving a gelatin capsule as the inclosure for the main com- 
at th bustible material, the mass (“‘dry’’) of the gelatin capsule was first 
affect determined. Then the solid sample was packed into the capsule and 
tion 0! the cover put on. The cover of the capsule was sealed to the body 


by dipping the entire capsule halfway under water for several seconds. 


po lit’s A 
‘ar bon “FP Sharp and Dohme, Baltimore, Md. 
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TABLE 2.—Tolerances of the amounts of close-boiling isomeric impurities in ji, 
octanes 


Amount of 
each such 
impurity 
required to 
Compound Close-boiling isomers change the 
heat of com- 
bustion by 
five parts 
in 100,000 


Mole 
fraction 
n-Octane (I) 
n-Octane (II) 
3-Methylheptane 0. 10 
2-Methylheptane.__. -Methylheptane OS 
3-Ethylhexane " . 06 

-Methylheptane 10 
3-Methylheptane -Methylheptane . 37 

| .14 
-Methylheptane . 08 
4-Methylheptane __. ? Sr eenne | . 37 
3-Ethylhexane 22 
‘ -Methy lheptane | . 06 
a. Meh 3-Methylheptane 
3-Ethylhexane --1)4-Methylheptane 
3-Methyl-3-ethylpentane 
2,4- Dimethylhexane 
2,2- Dimethylhexane 2,5-Dimethylhexane 
2,2,3,3-Tetramethylbutane 
3,4-Dimethylhexane 
2,3-Dimethylhexane 2-Methyl-3-ethylpentane 
2,3,3-Trimethylpentane 
2,4- Dimethylhexane ot Seer SS me 
2,5-Dimethylhexane -Dimethylhexane 
3-Trimethylpentane 
4-Trimethylpentane 
fethylheptane 
fethylheptane 
fethylheptane 


3,3-Dimethylhexane f 
ethylhexane 


3,4-Dimethylhexane 


Dimethylhexane 
Methyl-3-ethylpentane 
,»3-Trimethylpentane 
‘thylhexane 
-Dimethylhexane 


3 
5 
N 
3 


2-Methyl-3-ethylpentane 


‘ 
N 
N 
F 
N 
5 
k 
N 


3-Methy]-3-etbylpentane 4 
Methy]-3-ethylpentane 
.4-Dimethylhexane 


9 - 
2 ,5- Dimethylhexane 


3-Trimethylpentane 
i-Trimethylpentane 
3- 


Trimethylpentane 


2,2, 
2,2, 
23 .3,4-Trimethylpentane 
— ,3- Dimethylhexane 
,3-Dimethylhexane 
3,3-Trimethylpentane 





2,3,4-Trimethylpentane 
2,2,3,3-Tetramethylbutane ‘im »2- 
) 9 


Dimethylhexane 
4-Trimethylpentane 


-Estimated; see text. 


For the solid octane, two types of experiments were performed. | 
the first type were measured the temperature rise of the calorimeter 
and the mass of carbon dioxide formed in the combustion of a mass 
(not accurately known) of 2,2,3,3-tetramethylbutane together with a 
known mass of gelatin capsule. In the second type of experiment were 
measured in one series, I, the temperature rise of the calorimeter and 
the mass of carbon dioxide formed in the combustion of a known mass 
of gelatin capsule, and in another series, II, the temperature rise of 
the calorimeter and the mass of carbon dioxide formed in the com- 
bustion of a known mass of gelatin capsule together with a known 
mass of pure benzoic acid of known heat of combustion. The second 
type of experiment served to determine the amount of heat and of 
carbon dioxide evolved in the combustion of unit mass of gelatin, the 
I series involving the determination directly on 13 capsules (1 capsule 
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in experiment No. 1, and 4 capsules each in experiments 2, 3, and 4 
of table 4), and a II series involving the determination indirectly, 

but under conditions simulating more closely the combustion experi- 
ments, with two experiments (Nos. 5 and 6 of table 4) involving one 
capsule and a known prein A amount of pure benzoic acid (NBS 
Standard Sample 39 e). 

The average values obtained for the amount of heat evolved and 
the mass of carbon dioxide formed in the combustion of unit mass of 
gelatin capsule were combined with the data of the first type of experi- 
ment to obtain the te mperature rise of the calorimeter and the mass 
of carbon dioxide formed from the combustion of the 2,2,3,3-tetra- 
methylbutane alone. 


2. CALORIMETRIC PROCEDURE 


The calorimetric procedure was the same as that in the work on 
the hexanes [2] and heptanes [3]. 


V. RESULTS 


The experimental results of the present investigation are presented 
in tables 3, 4, and 5, which give the data yielding values of B in ohms 
per gram of carbon dioxide, as defined by eq 4 of the paper on the 
hexanes [2], for 17 of the octanes in the liquid phase, and for the eight- 
eenth octane (2,2,3,3-tetramethylbutane) in the solid phase, at 28° C 
The symbols have the same significance as in previous re “ee (2, 12]. 

In the experiments on 2-methylheptane, there was a loss of carbon 
dioxide during the analysis of the products of combustion of the calo- 
rimetric combustion experiments. In order to save repeating the 
latter experiments, a series of four noncalorimetric combustion expe- 
riments was performed to determine the ratio of the mass of carbon 
dioxide formed to that calculated stoichiometrically from the mass 
of sample. These four experiments yielded the following value for 
thisratio: Mean, 0.99970; standard deviation of the mean, + 0.00006. 
The values of the mass of carbon dioxide formed as given in table 3 
are calculated from the foregoing ratio. The uncertainty assigned 
to the final value for 2-methylheptane was increased by the uncer- 
tainty in the stoichiometric ratio. 

In table 6 are given, for the 18 octanes, values of the following: 
The constant B in ohms per gram of carbon dioxide, as given in 
tables 3 and 5; the ratio of the heat evolved, per mole of hydrocarbon, 
in the bomb process at 28° C for each isomer to that for n-octane; 
the difference, between n-octane and each isomer, in the heat of 
combustion in the bomb process at 28° C; and the heat of isomeriza- 
tion of n-octane into each isomer, at 25° C and 1 atm, for the liquid 
state for the 17 normally liquid octanes and for the solid state for 
2, 2, 3, 3-tetramethylbutane. The value of the heat of isomerization 
was obtained by means of the relation [2] 


(- AU sz)», : AU ps); (—AUsz),(1—B,/B,). 


For this calculation, the value of (—Al’,), at 28° C was taken as 
458.8 int. kj/mole [4]. 
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The values in column 6 differ from those in column 5 because of 
the difference in the heat capacities of the isomers over the range 25° 
to 28° C, the Washburn conversion to zero pressure [20), and the cop. 
version of AL to A(U+-PV)=AH. Account was taken of the effec 
of the gelatin and of benzoic acid on the Washburn correction of th 
results of the experiments given in tables 4 and 5. 

In table 7 are given values for the following: The heat of i isomeris Z 
tion, for the condensed phase, of n-octane into each isomer, at 2 
and 1 atm, as given in table 6; the difference, between mentee and 
each isomer, in the heat of vaporization (sublimation for 2, 2, 3, 3. 
tetramethylbutane), at 25° C, as taken from the work of Osborne 
and Ginnings [5], reduced to the ideal state [10]; the heat of isomeriza- 
tion, for the ideal gaseous state, of n-octgne into each isomer at 25° ( 
the difference, between n-octane and each isomer, in the value of 
Hi3o3.13—H3, the heat content at 298.16° K (25° C) referred to the 
absolute zero, for the ideal gaseous state [6]; and, finally, the energy of 
isomerization, for the ideal gaseous state, of n-octane into each isomer 
at 0° K. 


TABLE 7.— Heats of isomerization of the octanes in the gaseous state * 








| Difference 
|} inthe Heat of iso- 
| Heat of isomerization of | heats of | merization | /4%s.10— 
the condensed phase ° | yaporiza- | of the gas H5) a— 
Substance | H? (liq or solid)— fon e | H? (g)— 
| 


| (Hy; - 
H§ (liq) at 25°C — | (AHe®)a—| (ey | Hos 
| Corey | at asec | 

= 





int. kj/mole | kcal/mole | kcal/mole | kcal/mole kcal/mole | kcal/mole 
n-Octane.. : 0. 00 | 0.00 0. 00 0. 00 0. 00 0.00 
2- Methylheptane ~---.------| —5. 340. 92)—1. 2540. 22) . 4320. 01/—1. 6820. 22) —. 5040. 20)—1. 004.0.% 
3-Methylheptane. - . “ns 2. 56-0. 78|—0. 6140.19} .39-+0. 01/—1.00+0. 19) —. 694-0. 25|—0. 314 
4-Methylheptane. -. | - 8240.79) —.44+0.19) .43+0.01)—0. 8720. 19] —. 772-0. 25) —. 1040.31 
3-Ethylhexane. --_. eS . 58+0. a —. 14+0. a .44+0.01) —. 58240. 17)—1. 2040. 40) +. 624/ 
2-Dimethylhexane._____-- . 10+0. 65) —2. 89-0. 16) 1. 00+0. 01/—3. 89-0. 16| —0. 99-0. 35) —2. 904.0. 8 
a Dimethylhexane ¢.____.-. 2. 80-1. 17|—0. 6740. 28 64+0. 01/—1. 310. 28) —. 8840. 30)—0 430 
4-Dimethylhexane.._-__.-..-- . 22+0. 76)—1. 7340.18) . 89-0. 01|—2. 6240. 18}—1. 210. 40)/—1. 412 
2, 5-Dimethylhexane.________- . 60-1. 17|—2. 53-40. 28]. 86-0. 01)—3. 39-0. 28]—1. 114-0. 35|—2. 2 0.8 
Dimethylhexane 7. 76-0. 70|—1. 850. " - 940. 01) —2. 790. 17) —1, 220. 40)—1. 5720.8 


_o 

4-Dimethylhexane.._....... . 0441. 23)—0. 490. - 60+0. 01|—1. 090. 29]—0. 90-0. 30) —0. 1940.22 
- Methyl-3-ethylpentane..__ -- 23-0. 94) +-. 0540.22) .71+0. 01/—0. 6640. 22/—1. 0340.35) +. 3740.41 
- Methyl-3-ethylpentane..___. 3. 0440.88) —.7340.21| .81+0. 01/—1. 5640. 21/—1. 2140. 40) —. 3540.4 
2, 3-Trimethylpentane..__.-. 7. 10+1. 16|—1. 70+0. . 09-40. 01]—2. 79-40, 28/—1. 35-40. 45/—1. 440.8 
2, 4-Trimethylpentane . 38-+1. 03) —2. 24+0. 25) 1. 5140. 01!/—3. 75-0. 25) —1. 354.0. 45) —2. 4040.82 


3, 3-Trimethylpentane.... 3. 7241. 07|—0. 89--0. 26) 1. 0240. 01/—1. 9140. 26|—1. 23-40. 40|—0. 6820.8 

, 3, 4-Trimethylpentane..-_.- . 2541. 34) —1. 254-0. 32} 0. 90+0. 01) —2. 1540. 32 . 36-40. 45) —. 7940.5 

2, 3, 3-Tetramethylbutane._|—18. 841. 57|—4. 50+0. 38|}—. 334.0. 01) —4. 170. ; . 41+0. 45 —27 6+0. 8 
] 





® The uncertainty to be assigned to the value of the heat of isome rization of any one of the isom ers into any 
other one may be taken as + 0.25 kcal/mole for the condensed phase at 25° C, 5 kcal/mole for the gas ® 
25° C, and +0.50 kcal/mole for the gas at 0° K, except for 2,2,3,3-tetrame thylbutane for whieh the uncertainty 
is somewhat larger, and except for 2,3-dimethylhexane, for which see footnote ‘‘d”’ below. 

> All are liquid except 2, 2, 3, 3-tetramethylbutane, which is solid. 

¢ For the solid 2, 2, 3, 3-tetramethylbutane, the heat of sublimation is involved. 

4 The values assigned to 2,3-dimethylhexane may be uncertain by several times the amount indica! 
because of the uncertainty as to purity (see tables 1 and 2) 





Heats of Isomerization of Octanes 


VI. DATA OF PREVIOUS INVESTIGATIONS 


The only published papers describing the determination of the heat 
of combustion of one or more octanes are the following: 


iid Vita i | 
Investigators Year Reference Compounds 


1898 {14} n-Octane 


n-Octane 

2-Methylheptane 

1910 {15] 3-Ethylhexane 

2,5-Dimethylhexane 

3,4-Dimethylhexarie 
fn: -Octane 

[16] | 2,2,4-Trimethylpentane 

(202 3,3-Tetramethyibutane 
Banse and Parks. . 3c fi n-Octane 
37 18 | n-Octane 


Jessup 


Of the foregoing data, only those of Richards and Jesse [15] and of 
Jessup [16] yield values for the heats of isomerization of any of the 
octanes. 

The data of Richards and Jesse [15] yield the following values for 
the ratio of the heat of combustion of a given isomer to the heat of 
combustion of n-octane: n-octane, 1.0000; 2-methylheptane, 1.0010; 
3-ethylhexane, 0.9983; 2,5-dimethylhexane, 0.9990; 3,4-dimethylhex- 
ane, 0.9993. The corresponding values of heats of isomerization for 
he four isoparaffins become, for n-octane into a given isomer, in the 
liquid state, AZ=-+1.3, —2.2, —1.3, —0.9 kcal/mole, respectively. 
The corresponding values from the present investigation are —1.25, 
—().14, —2.53, —0.49 keal/mole, respectively. Since the values from 
Richards and Jesse may be uncertain by as much as 2 kcal/mole or 
more, it is seen that their data are substantially in accord with those 
of the present investigation within the respective limits of uncertainty. 

The data of Jessup [16] were obtained in an investigation in which 
circumstances made it possible to perform only two “complete com- 
bustion experiments on each of three isomers. His data yield, for 
ne heat of isomerization of liquid n-octane into liquid 2,2,4-tri- 

iethylpentane and into solid 2,2,3,3-tetramethylbutane, AH=—1.8 
ad —3.2 kcal/mole, respectively. The corresponding values from 
he present investigition are —2.24 and —4.50 kcal/mole. Because 
of the then unavoidable limitations placed upon these experiments 
of Jessup, the values of heats of isomerization deduced from them are 
arbitrarily assigned an uncertainty of +1.0 kcal/mole. It is seen 
that Jessup’s values are in accord with those of the present investiga- 
tion within the respective limits of uncertainty. 


VII. DISCUSSION 


\s has been previously pointed out, values of the heats of isomeri- 
zation of the octanes, along with those of the other paraffin hydro- 
arbons, are needed in calculating, in combination with free-energy 
functions, values of the free energies of isomerization and in studying 
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the influence of structure upon the energy content of the molecule. 

With regard to the correlation of the energy content with structure. 
a separate report will be made. It may be pointed out here, how- 
ever, that the data on the energies of isomerization of the octanes 
are in accord with the generalization previously made [1, 3] concerning 
the change with structure of the energy content of hydrocarbons, 

The use of these data on energies of isomerization for calculating 
equilibria and free energies of isomerization of the 18 octanes is dis. 
cussed in another report [13]. 
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REGION OF USABLE IMAGERY IN AIRPLANE- 
CAMERA LENSES 


By Francis E. Washer 


ABSTRACT 


The proper positioning of a lens with respect to the focal plane in a fixed-focus 
camera, such as an airplane-mapping camera, is governed by several factors that 
relate to the optical qualities of the lens. These factors are quality of imagery, 
depth of focus at a given stop opening, and curvature of field. This investigation 
shows that there is reasonably good agreement between observed depth of focus 
at a given stop opening and that predicted on the basis of geometric optics. Ob- 
served values of the maximum resolving power at various angular separations from 
the axis are generally lower at the larger stop openings than values predicted on 
the basis of physical optics. This lowering is doubtless a consequence of residual 
aberrations, inherent in an actual lens, which are more noticeable at large aperture 
ratios. A method of presenting the resolving power characteristics throughout 
the range of useful imagery in the form of sets of master curves has been developed. 
These curves show at a glance the variation of resolving power with distance from 
the plane of best axial imagery, the depth of focus for any observed value of the 
resolving power, the effect of field curvature on imagery in any given image plane, 
and the differing performance for tangential and radial imagery at various angular 
separations from the axis. Successive groups of these master curves show how 
variations in stop opening affect the performance of a lens. 


CONTENTS 


. Introduction 
Method 
Theory = een _— saaeee 

1. Maximum resolving power of an ideal lens 

2. Depth of focus for an ideal lens- 
Results of measurement_- - 
l. Variation in depth of focus on the axis with stop opening - 
2. Variation in resolving power on the axis with stop opening 
3. Effect of field curvature 
t. Change in depth of focus with angular se eparation from the axis- 
». Change in maximum resolving power with angular separation 


from the axis_ - anne : . 
6. Variation of resolving power in a given image plane with stop 
opening ; 


I. INTRODUCTION 


In the use of aerial photographs for mapping, the quality of defini- 
on is of considerable importance. Since quality of definition in a 
hotograph is directly related to the resolving power of the lens used 
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in the aerial camera, the factors affecting resolving power throughout 
the entire area of the image plane are of special interest. These fae. 
tors are depth of focus, stop opening, and curvature of field. Al] of 
these factors are considered in the proper positioning of the lens jy 
the camera. For a lens free from aberration, the image plane of best 
definition is unique and its position can readily be determined theo. 
retically, but for an actual lens it is necessary to resort to experiment 
and to compromise upon a position of the lens that yields best average 
definition throughout the image plane. This is usually done at the 
maximum stop opening of the lens, and this setting is then held, jy 
the fixed-focus airplane camera, for all smaller stop openings. 

Since the position of best average definition for a lens is generally 
determined with the lens not mounted in a camera, and since the lens 
may then be mounted in such manner that the distance from rea; 
vertex of lens to image plane (back focal length) is not that prescribed 
for best average definition, it is of interest to know how great a depar- 
ture from the prescribed position is allowable before the deterioratio: 
of imagery becomes sufficient to impair the usefulness of the camers 
These tolerances are best determined from information on the dept! 
of focus and curvature of field of lenses used in practice. The present 
paper therefore shows the results of measurement on resolving powe: 
for several typical airplane camera lenses. This information is pre- 
sented in such manner that the effect of depth of focus, stop opening 
and curvature of field is readily apparent. 


II. METHOD 


The precision lens testing camera! is well suited to the study of 
the resolving power characteristics of a lens throughout the region of 
usable imagery. The only modification in procedure from that pre- 
viously reported ? is the increasing of the range through which th 
camera back moves in making a resolving power negative. Fifty- 
seven steps instead of nineteen are made, and the negative shows the 
resolving power at definite intervals from 6.5 mm nearer to the lens 
to 6.5 mm farther from the lens than the plane of best axial focus 
The negatives are examined under a microscope and the resolving 
powers recorded for each image. A pattern is not considered resolved 
unless all coarser patterns are also resolved and the number of lines 
in the resolved patterns agree with that of the corresponding patterns 
in the test chart. In this manner, the values of the resolving powe' 
are determined for each of the 57 image planes at every 5° interval 
from the axis from 0° to the limiting angle set by the lens. By making 
such a negative for each stop opening, the region of usable imagery 's 
explored throughout the range of aperture selected. 

Although a large number of lenses have been studied in this manne! 
the results of only four are here presented, and of these four, ‘wo 
are selected for a more complete description. The lenses selected 
are representative of a wide range of focal lengths, angular fields, and 
function in mapping. Lens A is a wide angle airplane-mapping lens. 
nominal focal length is 150 mm, maximum aperture //6, and a useill 
field of view 45°. Lens B is a lens usable in copying cameras or [0 
general photography; it has a nominal focal length of 195 mm, max- 


mum aperture //5.6, and a normally useful field angle between 2 


1J. Research N BS 18, 449 (1937) RP984. 
*J. Research NBS 22, 729 (1939) RP1216. 





Region of Usable Imagery 177 


and 25°. Lens C is usable in copying cameras or for general photog- 
raphy; it has a nominal focal length of 210 mm, maximum aperture 
‘6.8 and a normally useful field angle of 25°. Lens D is frequently 
used in aerial cameras and also in copying cameras; it has a nominal 
focal length of 125 mm, maximum aperture f/8, and a normally useful 
field angle of 35°. Lenses A and & have been selected for the fuller 
description of their resolving power characteristics whereas C and D 
are limited to a partial description. 


III. THEORY 


A brief discussion of the resolving power characteristics as pre- 
dicted from theoretical considerations is necessary before proceeding 
to the results of measurements. This theory has long been estab- 
lished and predictions from it will be used in this article only for the 
purpose of providing a yardstick against which to compare the per- 
formance of the lenses studied. Since the predicted values of resolv- 
ing power are for an ideal lens, it is not surprising that some apparent 
contradictions to the theory will be found in the experimental results. 
These deviations result from aberrations inherent in any actual lens. 
These departures do not however signify that the theory is without 
value in the prediction of lens performance. On the contrary, the 
uumber of times that the observed data are of the same order of 
magnitude as the predicted values is sufficiently high that one may 
feel more confident than ever that the predicted values are sufficiently 
close to the true values as to render them a valuable aid. This is 
particularly true when dealing with resolving powers of the same 
order of magnitude as those of the more commonly used photographic 
emulsions. 


1. MAXIMUM RESOLVING POWER OF AN IDEAL LENS 


The smallest separation of two star imeges theoretically obtainable 
on the basis of diffraction of light is 
1.22F 


where X is the separation, F' the focal length of the lens, d the dia- 
meter of the lens and \ the wavelength of the incident light. For the 
purpose of the present work, it will be assumed that the same equa- 
tion is valid where the images are composed of parallel lines instead 
of points. It is also more convenient to express the maximum re- 
solving power, A,, in lines per millimeter. It can be shown that the 
equation, 


a] 


A,= lines per millimeter, 


where 5 is the f-number, is equivalent to eq 1. This relation is com- 
puted for a wavelength of 575 millimicrons, which is the effective 
wavelength of the light used in making resolving power negatives with 
the precision lens-testing camera. This equation is valid only on the 
axis of the lens. For extra axial points, it is necessary to use the two 
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additional approximate relations, one for maximum resolving 
power for tangential lines, ; 


1426 
Ae=—j—cos 


and another for maximum resolving power for radial lines, 


26 
Ase= ee “cos 6, 


where @ is the angular separation from the axis of the point consider 
Table 1 gives the values computed for the resolving power in both 
orientations at 5° intervals from 0° to 45° for a series of f-numbers 


ranging from f/4 to {/64 in a form convenient for use. 


TaBLeE 1.—Mazimum resolving power for an ideal lens as a function of stop opening 
in lines per millimeter 


Resolving power for tangential lines. 
Angular distance from axis. 
Relative aperture & s 
or f-number 


265 231 

190 | 166 

132 116 

117 | 84 
80 | 66 58 
58 | < 48 42 
41 | 37 33 | 29 
29 26 24 21 
20 18 16 14 


Resolving power for radial lines. 
Angular distance from axis. 


j 
20 25 30° 


2. DEPTH OF FOCUS FOR AN IDEAL LENS 


The term “depth of focus’’ as used throughout this paper refers to 
the range measured along the axis in the image space that may be 
traversed without the resolving power measured in the image plan 
falling below an arbitrarily selected value. It is not to be confused 
with “depth of field,” which refers to the range in the object spac 
for which all point objects are imaged with a Tadius of the circle of 
confusion remaining less than or equal to some arbitrarily selected 
value. 
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RE 1.—Azial resolving power versus position of the image plane for an ideal 
lens at several stop openings. 


here the maximum value of the resolving power theoretically possible according to physical optics falls 

n the scale range of the graphs, a straight-line intercept parallel to the axis of abscissas is shown to 

te that it is not significant to extend higher the two curves, asymptotic to the focal plane, that are 

i from geometric optics. Where this maximum falls outside the scale range of the graphs the asymp- 

urves are left open at the top but are to be assumed to extend beyond the range of the graph to the 

tical maximum value of the resolving power, which is indicated in the upper right-hand corner of the 

me foreach stop opening. ‘The zero of abscissas marks the position ef best axial focus, and positive values 
issas indicate positions farther from the lens. 
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It can be shown that the magnitude of the depth of focus d jy 
millimeters can be approximately computed from the relation 


(6 


where 6 is the f/-number and a is the resolving power, in lines per 
millimeter, for which the depth of focus is sought. Eq 6 is based oy 
geometric optics and is restricted to values of a lying between zery 
and A,,, where A, is the maximum resolving power as computed by 
eq 2, which is based on physical optics. According to the theory of 
geometric optics the value of d is independent of angular separation 
from the axis for a given 6, and consequently no additional equations 
for off-axis conditions are necessary. 

With the aid of eq 6 supplemented by eq 2, it is possible to predict 
the manner in which the resolving power varies on the axis as the 
image plane assumes various positions along the axis in the vicinity 
of the focal plane. This has been done, and the results are presented 
in figure 1 for a series of apertures. It may be noted that values of 
the resolving power that are considered acceptable for lenses used in 
aerial photography, such as 20 on the axis and 7 at the wider angles 
of separation, are obtained throughout a considerable range of settings 
even at the largest aperture. 


IV. RESULTS OF MEASUREMENT 


In order to present a comprehensive picture of the performance of 
actual lenses throughout the entire region of usable imagery, the 
results of measurements on the two lenses A and B are shown in figures 
2 to 8. For lens A, figures 2, 3, 4, and 5 show the resolving power 
curves for both tangential and radial imagery for all marked apertur 
ratios from {/6 to {/45, inclusive, and at each 5° interval of angular 
separation from the axis from 0° to 45°, inclusive. In a similar manner 
figures 6, 7, and 8 show the resolving-power curves for lens B for all 
marked aperture ratios from {/5.6 to {/32 and at each 5° interval of 
angular separation from the axis from 0° to 30°, inclusive. For any 
one of these master curves, the ordinates give the resolving power in 
lines per millimeter clearly resolved; the corresponding abscissas give 
the axial distance from the plane of best axial focus. Negative values 
of abscissas signify that the points are nearer to the lens than the 
plane of best axial focus, which reference plane marks the zero of 
abscissas for all curves shown. The lens identification by symbol, 
focal length, and f-number is given in the upper-right-hand corner o! 
the uppermost frame for each group. The type of imagery and the 
angular separation from the axis are indicated in the upper left-hand 
corner of each frame. Tangential imagery is signified by the letter 7 
and radial imagery by the letter R. 





Region of Usable Imagery 














2 s “~ 7 


POS iN OF THE IMAGE PLANE 


Resolving power versus position of the image plane for lens A at apertures 
f/6 and f/8. 
lving power throughout the region of useful imagery is shown for tangential (T and radial (R) 


intervals from 0° to 45°. The zero of abscissa marks the position of best visual focus on the axis at 
sitive values of abscissas indicate positions farther from the lens 
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FiGuRE 3.—Resolving power versus position of the image plane for lens A at aper 
j/11 and f/16. 
The resolving power throughout the region of useful imagery is shown for tangential (7') and radial F 
lines at 5° intervals from 0° to 45°. The zero of abscissas marks the position of best visual focus on the ais # 
f/6, and positive values of abscissas indicate positions farther from the lens. 
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POSITION OF THE IMAGE PL ANE 


ving pr versus position of the image plane f 


at apertures 
f/22 and f/32. 


ing power throughout the region of useful imagery is shown for tangential (7') and radial (R) 
ntervals from 0° to 45° I'he zero of abscissas marks the position of best visual focus on the axis at 
tive values of abscissas indicate positions farther from the lens. 
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Figure 5.—Resolving power versus position of the image plane for lens A at ape’ 
145 
f/45. 


The resolving power throughout the region of useful imagery is shown for tangential (7°) and radia! 
lines at 5° intervals from 0° to 45°. The zero of abscissas marks the position of best visual focus on the 
f/6, and positive values of abscissas indicate positions farther from the lens. 
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Figure 6.—Resolving power versus position of the image plane for lens Bat apertures 
f/5.6 and f/8. 
The resolving power throughout the region of useful imagery is shown for tangential] (7') and radial (2) 


ues at 5° intervals from 0° to 30°. The zero of abscissas marks the position of best visual focus on the axis at 
and positive values of abscissas indicate positions farther from the lens. 
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Figure 7.—Resolving power versus position of the image plane for lens B at apert 
f/11 and f/16. 


lhe resolving power throughout the region of useful imagery is shown for tangential (7°) and radial 
lines at 5° intervals for 0° to 30°. The zero of abscissas marks the position of best focus on the axis at / 
and positive values of abscissas indicate positions farther from the lens. 


Thoughtful examination of the curves in figures 2 to 8 yields 3 
wealth of qualitative information on particular phases of lens perform- 
ance and on the manner in which these separate phases form ou , 
coherent whole. However, in view of the complexity of the con- locu: 
figuration when viewed as a whole, it is desirable to discuss in separat may 
sections the various aspects of lens performance that are here shows pow 
in the integrated form. Chief among these aspects are Cepth of focus lor a 
effect of stop opening, variation with angular separation from th of a: 
axis, and field curvature. In these sections data derived for lenses ( the 
and D from sets of curves similar to those in figures 2 to 8 (but noi ment 
shown to conserve space) are presented and discussed in conjunctiol dept! 
with the derived data for lenses A and B. In these discussions, tht is pr 
experimental results are also compared with those predicted from taine 


theory. rt 
catio 
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Ficure 8.—Resolving power versus position of the image plane for lens B at apertures 
D2 2 
f/22 and {/32. 


The resolving power throughout the region of useful imagery is shown for tangential (7) and radial (R) 
I ntervals from 0° to 30°. The zero of abscissas marks the position of best focus on the axis at //5.6, 
j p sitive values of abscissas indicate positions farther from the lens. 


1. VARIATION IN DEPTH OF FOCUS ON THE AXIS WITH STOP 
OPENING 


When an infinitely distant object is imaged by a lens, the depth of 
focus is a measure of the range of positions along the optic axis that 
may be assumed by the image plane without reduc tion of the resolving 
power below a specified value. In figures 2 to 8 the depths of focus 
for any observed resolving power is simply the length of the segment 
of a straight line parallel to the axis of abscissas cut off by the curve at 
the selected value of resolving power. Thus the length of this seg- 
ment for the maximum observed resolving power gives directly the 
depth of focus in millimeters. However, for comparison purposes, it 
is preferable to select some value of the resolving power that is at- 
tained at each stop opening for the lens consider ed. 

The least axial resolving power allowable under most federal specifi- 
cations dealing with lens performance i is 20 lines per millimeter. It is 
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therefore pertinent to study the range throughout which this value js 
found and also to note the effect of stop opening. In figure 9, the 
variation of depth of focus with stop opening is shown for the four 
lenses A, B, C, and D for the resolving power of 20 lines per millimeter. 
or as near thereto as feasible. The actual value of the resolving 
power for which the measurements are made is indicated on curve | j; 
each frame. 

The values used for plotting curve 1 are computed with the aid 
of eq 6 for the indicated resolving power. Curve 2 is plotted from 
the experimental data, each circle on the curve showing the observed 
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Figure 9.-—Azial depth of focus versus f-number for a fixed resolving power 


Curve 1 shows the computed relation for an ideal lens, and curve 2 shows the observed relatior 
actuallens. The value of the fixed resolving power in lines per millimeter is indicated on curve | for each k 


value of d for the given f-number. It is noteworthy that the depth 
of focus for an actual lens does not always increase with diminishing 
aperture. In fact in some instances it may decrease, as is seen 
the curve for lens A. Here the depth of focus decreases from /' 
to f/11, and then begins to increase and continues to increase iD § 
linear fashion but at a slope less than that of curve 1. Thus # 
{/6 and f/8 the depth of focus is greater than that predicted, whereas 
at all smaller f-numbers it is substantially less than predicted. Th 
predicted values are not exceeded for lenses B, C, and D. However 
the increase with diminishing aperture appears to become more ul 
form at f/11 for B and C, and at f/16 for D. 

It is clear from these four examples that the increase of depth 
focus with diminishing aperture does not proceed in regular fashiot 
until an aperture appreciably smaller than the maximum is reached 
The value of the aperture ratio that marks the beginning of regu! 
increase is usually f/11, although in some instances it begins as eat! 
as {/8 and in other instances as late as f/16. However, on the bass 
of many other examples not reported here, f/11 is the transition poll! 
in the majority of cases studied. In this region of uniform increas 
the depth of focus, while increasing steadily with diminishing ape 
ture, is nonetheless lower than that predicted by eq 6 and in mo 
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cases 15 percent or more lower. The anomalous behavoir at large 
aperture ratios, such as is shown by lens A at f/6 and 7/8, can only 
he accounted for by the presence of lens aberratiors, which in some 
nstances may so broaden the base of the resolving power curve that 
values of the depth of focus greater than predicted are found. If the 
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POSITION OF THE IMAGE PLANE 
iURE 10,—Azial resolving power versus position of the image plane for a lens 
possessing multiple foci at several stop openings. 


f abscissas marks the position of best visual focus on the axis at //6.8, and positive values of 
ndicate positions farther from the lens. 


axial imagery of such a lens as A be examined either visually with 
# microscope or photographically with finer test charts, then it is 
probable that two or more foci will be seen at full aperture. The 
axial separation of these foci may vary from a few tenths of a milli- 
Meter to more than a millimeter. 

lhis phenomenon has been observed on other lenses and in figure 10, 
resolving-power curves for lens E at f/6.8, f/8, f/11, and f/16 show the 
fect ia a striking manner. This lens shows three foci on the axis at 
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4/6.8, drops to two foci at /8 and to a single focus at f/11 and f/15 
It is evident that the outer zones are focused at a point nearer to the 
lens than the inner zones, the separation of foci in this instane 
amounting to nearly a millimeter. It is clear from the curves tha: 
if the depth of focus is measured, say at 17 lines per millimeter, th 
depth of focus will decrease with diminishing aperture so long as cop. 
ditions contributing to a multiple focus effect are present. Althoug) 
only this one example is given, it is by no means unusual to find ty 
or more foci on the axis with either visual or photographic means oj 
observation. In the present instance it is evident that the princip, 
focus is the one farthest from the lens, for it remains fixed with respec 
to the zero of abscissas for all aperture ratios shown. This effect may 
also account for the lowering of central definition with diminishing 
aperture which is sometimes observed. If the focal plane is locate 
by the maximum nearest to the lens, the resolving power will be 
markedly lower at f/11 and /16 than at //6.8 and #/8. 
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Figure 11.—Azial depth of focus versus line width (or reciprocal of resolving power 
in lines per millimeter) for aperture f/11. 


Curve 1 shows the computed relation for an ideal lens, and curve 2 shows the observed relat 
actual lens. 


From curves 2 to 8 it is clear that the depth of focus diminishes 
with increasing resolving power. To provide a more clear cut pictul 
of this effect, the curves in figure 11 show the manner in which th 
depth of focus varies with resolving power for a given f-number, 1! 
this instance {/11. The reciprocals of resolving power, or line spacings 
rather than values of resolving power are plotted as abscissas in orde! 
that the theoretical curve computed from eq 6 appear as a straigil 
line instead of a rectangular hyperbola. Curve 1 shows the variation 
predicted and curve 2 shows that actually found. Aside from the 
observed values of depth of focus being lower than the predicte 
there are no important deviations from theory. The curves do sho¥ 
the greater freedom in selection of a focal plane that exists for the 
lower resolving powers. 
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2. VARIATION IN RESOLVING POWER ON THE AXIS WITH STOP 
OPENING 


The upper limit of resolving power, expressed in lines per milli- 
meter, is determined by aperture ratio and wavelength of the image 
forming light. The capacity of a lens to yield this maximum resolving 
power 1s in practice impaired by lens aberrations. In the present 
instance, the patterns of the test chart are not sufficiently fine to 
show variations in the maximum resolving power at the larger stop 
openings. This is evident when it is noted that the maximum 
resolving power on the test chart for lens A is 78 lines per millimeter, 
which is higher than can be expected for {/22 and lower than predicted 
for {/16. For lens B, the test-chart maximum is 56, which is greater 
than the predicted value at {/32 and less than predicted at //22. 
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F-NUMBER 
FicurE 12.—Mazimum axial resolving power versus slop opening. 


irve 1 shows the computed relation for an ideal lens for a wavelength of 575 millimicrons, Curve 
shows the observed relation for an actual lens. 


? 


_ However, despite these limitation in the test-chart scale, some 
information on the change in resolving power with stop opening can 
be obtained from the curves in figures 2 to 8. To facilitate study, 
the maximum observed resolving power is plotted against f-number 
as curve 2 in figure 12 for lens A and B. Curve 1 in each frame 
shows the maximum theoretical resolving power. 

_ For lens A, it is evident that the maximum fixed by the test chart 
is obtained for f/6, {/8, and f/11. However, the resolving power 
drops to 56 at {/16, which drop can be attributed to lens aberrations. 
The value of 56 is still maintained at f/22, and since the theoretical 
is only 65, it can be stated that the lens aberrations definitely affect 
the resolving power less at f/22 than at f/16. The two curves are 
nearly in coincidence at f/32 and #/45, thus clearly showing the drop 
in resolving power with diminishing aperture as well as demonstrating 
that the drop in resolving power at small apertures results from 

627528457 
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theoretical limitations based on aperture rather than on the effe: 
of lens aberrations. 

In the case of lens B, the limitations imposed by the test chgy 
preclude any effect becoming noticeable for most of the range; hoy. 
ever, it is clear that lens aberrations are effective at f/22 and /') 
in preventing the lens from performing as predicted by theory. — 


3. EFFECT OF FIELD CURVATURE 


It is apparent on examination of a given series of curves at a spec. 
fied stop opening, for either tangential or radial imagery, that th 
position of maximum resolving power occurs at different distances 
from the position of best axial definition for the various anguly 
separations from the axis. This effect is a result of curvature of t] 
image field. The curves of primary and secondary curvature cx 
readily be visualized if, while examining figures 2 to 8, one holds thy 
page at an angle of 20° to 30° with the line of sight and looks along tly 
page from the bottom to the top. The central regions of the tangential 
curves then describe the curve of primary curvature, and similarly 
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Figure 13.—Field curvature. 


Curve P shows the primary curvature and is the locus of the positions of maximum resolving power 
tangential lines at the indicated angular separations from the axis. Curve S shows the secondary curvs'® 
and is the locus of the positions of maximum resolving power for radial lines at the indicated angular 
rations from the axis. The zero of abscissas marks the position of best axial focus and positive values 0! 8° 
sas indicate positions farther from the lens. 
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the central regions of the radial curves show the curve of secondary 
curvature. 

The primary and secondary curvatures for lenses A, B, C, and D 
are shown in conventional form in figure 13. The curves, as drawn, 
represent an average value of the field curvature for several apertures 
and consequently do not portray precisely the observed field curvature 
at any single aperture. It was noted in the course of this investiga- 
tion that the shape of the field curves changes appreciably when 
soing from the largest stop opening to the smallest. However, it was 
not considered worthwhile to show this change in detail. 

It is this effect produced by field curvature that makes it pointless 
to search out an image plane that yields say 200 lines per millimeter on 
the axis, because so placing the image plane automatically fixes the 
resolving power obtainable for the 30° region at 14 lines per milli- 
meter or less. It is clear from the curves in figures 2 to 8, and 13 that 
no single image plane can include the position of maximum resolving 
power for all angular separations from the axis for either type of 
imagery. Consequently, selection of an image plane that will yield 
values of the resolving power sufficient for good usable photography 
in all parts of the image plane for both types of imagery is preferable to 
selection of the plane of best axial imagery. It may also be noted 
that this shift of the maxima of resolving power for the various angular 
separations from the axis justified setting the upper limit of the test- 
chart patterns to exceed that of the average emulsion but not neces- 
sarily to exceed the theoretical resolving power of the lens at the 
larger stop Opening. 


4. CHANGE IN DEPTH OF FOCUS WITH ANGULAR SEPARATION 
FROM THE AXIS 


Examination of the curves in figures 2 to 8 shows that for a given 
‘number and specified resolving power the depth of focus is not 
generally constant. This behavior is markedly different from that 
which would be manifested by an ideal lens for which the depth of 
focus is constant under these conditions. Accordingly, curves illus- 
trating this lack of constancy in depth of focus for the lenses studied 
are shown in figure 14. The resolving power in each case is 20 lines 
per millimeter or as near thereto as feasible, the actual value being 
indicated. The variation of depth of focus with angular separation 
from the axis is shown for a series of f-numbers for each of the lenses 
A,B, C,and D. The theoretical depth of focus on the axis is marked 
by a short line projecting from the axis of ordinates. It must be 
borne in mind that the ordinates for any value of the abscissas repre- 
sent only the range throughout which the particular resolving power 
for that type of imagery is observed, and the curves do not show 
the space relationship of these ranges with respect to one another. 
lt is possible to draw curves, showing the depth of focus, that pre- 
serve the space relationships but since this would involve the effect 
of field curvature, the change in depth of focus with angle would be 
somewhat obscured. 

lt is noteworthy that of these four examples only lens B preserves 
‘hear constancy of depth of focus for tangential and radial lines with 
angular separation from the axis, and that even here there is an 
abrupt change between 20° and 25°. Since this change is more 
Pronounced at the larger stop opening, it is probable that the increase 
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in depth of focus is a consequence of vignetting by the lens boundaries 
which produces a sharp reduction in effective aperture between 20° 
ond 25°. Too, it may be noted that in the region of 0° to 20° the 
‘heoretical maximum is not exceeded. Wide variations in the depth 
of focus with angle are shown by lenses A, C, and D, and the varia- 
tions for tangential lines make the greater changes. Since these 
changes in many instances carry the depth of focus above the theo- 
retical maximum it is reasonable to suppose that this is another 
manifestation of the multiple focus effect discussed in section 1. 

The fact that the depth of focus for tangential lines is usually 
sreater than that for radial lines is probably a consequence of the 
asymmetric use of aperture for tangential imagery, thus permitting 
the aberrations to increase the depth of focus at the expense of a 
lowered maximum observable resolving power. 


5 CHANGE IN MAXIMUM RESOLVING POWER WITH ANGULAR 
SEPARATION FROM THE AXIS 


Examination of the curves shown in figures 2 to 8 shows that the 
maximum resolving power decreases with increasing angular separa- 
tion from the axis. This effect is more pronounced for tangential 
than for radial lines, although definitely present in the latter. This 
effect is compounded from lens aberrations and limitations set by 
aperture ratio. For the larger stop openings the decrease evident 
from the curves arises mainly from aberrations because the upper 
limit set by the stop opening is well above the upper limit set by the 
test charts. The variation in the maximum value of resolving power 

.b angular separation from the axis is shown graphically in figure 

The solid line curves represent the variation in maximum 
rsolving power with angle as computed with the aid of eqs 4 and 5. 
The broken curves are the results of observation. 

It is interesting to note that, save for one exception, the resolving 
power decreases monotonically from a maximum on the axis to the 
last observed value for both tangential and radial lines. However, 
the resolving power for tangential lines decreases more rapidly than 
does the resolving power for radial lines. This behavior is of course 
predicted by theory, but it is interesting to note that this behavior 
is present for those stop openings where the least resolving power 
predicted for any angle is nonetheless higher than the upper limit 
of the test chart, as can be seen from the curves for f/6 and f/8. At 
fll, the effect of diminishing aperture is evident, particularly for 
tangential values, although even here it is clear that aberrations are 
chiefly responsible for the failure of the lens to yield full theoretical 
resolving power. At {/16, the limits imposed by stop opening become 
more evident, the theoretical and observed curves are closer together 
and their similarity becomes more pronounced. At //22, the effects 
of aberrations are practically nil and the two curves nearly coincide. 
This close coincidence of the two curves continues through f/32 and 
1/45, 

From the study of these curves it is evident that lens aberrations 
are of primary importance in limiting the maximum resolving power 
of lens A for the aperture ratios from f/6 to f/11. From #/16 on, 
the maximum values of the observed resolving power conform more 
and more closely to the values predicted by theory; thus these 
observations tend to validate the accuracy of eqs 4 and 5. 
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Ficure 15.—Mazimum resolving power versus angular separation from the am 
for a series of stop openings. 


The solid curve shows the computed relation for an ideal lens, and the broken curve shows the observed 
relation for an actual lens. T signifies tangential imagery and RF radial imagery. 
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6. VARIATION OF RESOLVING POWER IN A GIVEN IMAGE PLANE 
WITH STOP OPENING 


A general improvement in definition throughout the image plane 
is frequently encountered when the imagery obtained with a lens 
set at a small stop opening is compared with that obtained with the 
same lens set at a larger stop opening. At first thought, this behavior 
appears to be at variance with the fact that the resolving power 
of a lens for any given part of the image plane is a maximum at the 
r stop opening. However, because of curvature of the image 
field (discussed in section 3) the image plane for an actual lens can- 
not usually be so placed as to contain all the resolving power maxima 
for all angular separations from the axis. On the contrary, a given 
mage plane can include only a few of the resolving-power maxima 

1 lens set at maximum stop opening; and, as is evident from 

examination of the master curves shown for lens A at 7/6 in figure 2, 
the small depth of focus present for maximum resolving power at the 
various angular separations from the axis further serves to limit the 
number of resolving-power maxima obtainable in the image plane. 

If, as is the case for airplane cameras, the image plane is main- 
tained at a fixed distance from the lens, which distance is established 
as yielding best average definition at full stop opening, it is clear that 
when a smaller stop opening is used the depth of focus for each maxi- 
mum increases. Consequently, there is greater likelihood of the given 
image plane including more resolving-power maxima than it did at the 
tull stop opening and, as a result, the definition considered throughout 
the image plane will show improvement. As a matter of fact, some 
mprovement is likely even if more maxima are not included because 
the increased depth of focus for any observed resolving power makes it 
probable that some regions will be benefited by an increase in resolving 
power even if the maximum is not reached. The manner in which this 
effect is operative can be seen from the curves shown in figures 2 to 8. 
For lens A, if one selects an image plane 0.20 mm nearer the lens than 
he plane of best axial focus, the resolving power for tangential lines at 
5° from the axis increases from 28 to 56 lines per millimeter when the 
aperture is reduced from f/6 to f/11. Similar gains can be noted in 
other parts of the field. 

This improvement does not continue indefinitely with diminishing 
stop opening. For lens A, reducing the f-number below //22 causes a 
general decrease in resolving power throughout the image plane. 
This arises from the fact that the limitations set by aperture radio 
described in section 2) becomes noticeable. For apertures between 
‘6 and f/22, the upper limit of the resolving power for the various 
regions of the image plane is set preponderantly by lens aberrations; 
but for f/32 and f/45, the upper limit is set by aperture ratio. 

It must be mentioned that there are additional reasons for the 
betterment of photographic quality with diminishing stop opening. 
Chief of these is the reduction of the large variation in effective ex- 
posure over the photographic plate which otherwise would be produced 
by vignetting. The difference in relative illumination is of course 
much less when the diaphragm stop is the only limiting aperture.® 

Wasnineton, November 14, 1944. 
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RELAXATION OF STRESSES IN ANNEALING GLASS 
By Arthur Q. Tool 





ABSTRACT 


An empirical equation representing relaxation of stresses in annealing glass 
jis derived. Although the derivation is based on Maxwell’s equation for viscous 
flow, consideration is also given to the changing viscosity as a glass anneals at 
a constant temperature. ‘This new equation has been applied to data that were 
obtained by Adams and Williamson on the relaxation of stresses in annealing 
glass at various temperatures. It is shown that this equation applies very 
satisfactorily to their data. Also, it appears to be much more suitable than 
the reciprocal relation which these authors proposed as a substitute for the 
wholly unsatisfactory exponential relation that is derived directly from Maxwell’s 
juation. This exponential relation is unsatisfactory simply because it does 
not take into account the changing viscosity of annealing glass. Finally, the 
new equation leads to a clearer understanding of the nature and behavior of 
glass because it does not neglect the change in viscosity. 
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I. INTRODUCTION 


In procuring data for a previous publication [1],' it was noted that 
the relaxation time ? of the strains in glass during an annealing at a 
onstant temperature was generally much less when the relative re- 
lardations exceeded a thousand millimicrons than when they were 
under fifty. This should not be the case if the viscosity remains 
prastant as long as the temperature is constant and if the stresses 

elax according to Maxwell’s law. At the time, it was thought possible 
iat the observed results might indicate that glass, even at relatively 
high annealing temperatures, is appreciably plastic. On further 
ousideration, the degree of plasticity indicated by the results seemed 
inre asonable, although glass is undoubtedly plastic to a certain degree 
at very low annealing temperatures. It also seemed improbable that 
he indicated increase in the relaxation time with decreasing load was 
be result of observational errors, although the method of measuring 


[Yaseen 
Figures in brackets indicate the literature references at the end of this paper. 


‘0 a case in which the stress, F, relaxes according to the relation, F= Foe~*/d, the constant, d, is the 
“aration time and F¢ is the stress when the time ¢=0. 
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the large relative retardations was comparatively inaccurate. More. 
over, a more extensive set of data showing practically the same cop. 
ditions was obtained by Adams and Williamson [2] at about the same 
time. 

These results show either that Maxwell’s law does not apply to the 
relaxation of stresses in annealing glass, possibly because the glass jg 
plastic, or that the viscosity of an annealing glass at a constant 
temperature changes and thus causes the factor representing the re- 
laxation time in Maxwell’s equation for inelastic flow to be a variable 
On finding that Maxwell’s equation was inadequate, Adams and 
Williamson used a reciprocal relation to express the decay of the 
relative retardation. That this subject is still a matter of some interest 
is shown by the fact that Weyl and Pincus [3] have recently mentioned 
the results of Adams and Williamson and their use of the reciprocal 
relation. For this reason a recently presented theory [4, 5] of th 
behavior of annealing glass has been applied to their data and the 
results are presented in this paper. 


II. CHANGING VISCOSITY OF ANNEALING GLASS 


In any thorough study of annealing glass, it will soon be found that 
all of the properties of the glass are changed by the annealing, even 
though the temperature is unchanged. Moreover, it will be found 
that the properties are a function of the annealing temperature as well 
as of the actual temperature at which the glass is tested. That is, if 
the annealing is continued until the equilibrium condition is reached 
at any annealing temperature, this temperature as well as the standard 
temperature at which the glass is tested should always be given con- 
sideration. Also during an annealing at a constant temperature and 
before equilibrium is reached, the condition .of the glass * at any time 
is always that which approximates an equilibrium condition at a 
higher or lower temperature which, as the annealing proceeds, con- 
tinually approaches that of the annealing [6). This equilibrium, or 
fictive, temperature, which continually changes while the glass is an- 
nealing at an actual temperature, 7’, will be designated by rt. Further- 
more, the time rate of change in 7 is always roughly proportional to 
T—r if the annealing is accomplished while 7 is constant. 

If the coefficient of a property of a glass, when 7 is constant, ca 
related reasonably well to the temperature by some expression, such 
as D=D,(1—a(T—T,)—. . .) in the case of density, for example, it 
will usually be found that a similar expression with different constants 
relates the coefficient equally well to the fictive temperature, r, whe! 
7 is constant. Thus, in the above relation for the change in density 
with 7, the variation in the initial density because of any change 11" 
is expressed equally well by the relation Dp=D’)(1—a(r—7)—. . 
As a result of his studies of annealing glass, Twyman [7] conclud 
that the inelastic deformability (reciprocal of viscosity) of the average 
glass doubled approximately with each increment of about 8° C in th 
actual temperature. Consequently, the relaxation time, #, which § 
closely related to the coefficient of viscosity, should vary approxi- 
mately according to the expression, )=0d,e~7" as the temperatur 
changes, if k=11.5° C. That is, although the modulus of elasticity 
E, also changes with temperature, its variation is doubtless insigni!- 
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3 Certain glasses behave in an unusual manner, but these are excluded from the present discussio0 
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cant when compared to the very rapid change in the coefficient of 
viscosity, n= Ee. If the changes in the relaxation time with 7 and r 
can be expressed reasonably well by similar relations, it follows that 
K=1/8=Kye"/%e" ~should express approximately the value of the 
reciprocal of the relaxation time for different values of JT andr. Ifa 
glass is cooled or heated so slowly that it never departs materially 
from the equilibrium condition, T=r, it follows that k=gh/(g+A), 
provided it is assumed that Twyman’s empirical relation was based 
on measurements on deformations that developed almost entirely 
while the glass was never far from equilibrium. A study of his experi- 
ments suggests that this assumption is applicable, since the deforma- 
tions which he observed almost certainly were developed mainly at 
comparatively high annealing temperatures at which 7 never departs 
materially from 7 unless the cooling rate is quite high. That is, the 
deformations were developed mainly before the glass departed from 
the equilibrium curve of its liquid form. 

If the expression for 1/8 is introduced into Maxwell’s equation 
for inelastic flow [1, 2}, a modification, dF/di=EdS/dt—K,Fe"’9 e*”, 
is obtained that should apply reasonably well to the development of 
inelastic deformation in stressed glass, provided the variations in EF 
with 7 and 7 are not too great. 

As it has been observed [8] that the manner in which temperature 
affects the velocity of approach by any coefficient (density, refrac- 
tivity, ete.) to its equilibrium value is roughly the same (provided the 
coefficient has a given displacement from its equilibrium value regard- 
less of the temperature) as the manner in which temperature aifects 
the mobility of the glass, it seems that the viscosity governs the rate 
of these approaches. Consequently, it appears that 


dr/dt=K,(T—r)e™ Ger/h (1) 


should represent approximately the velocity of the change in 7 as 
both T and + change. In part at least, this differential equation 
allows for the lack of strict proportionality between dr/dt and T—r 
in an annealing glass at a constant temperature. Its solution (in 
terms of exponential integrals) 


Eii(y) — Ei(yo) = — Koe™t (2) 


if T is constant and if y=(7—r)/h and yp is its value when t=0) 
applies rather satisfactorily to data on the changes in 7 as they are 
derived from experimental results on the varying coefficients of an 
annealing glass. 

lf a glass is cooled very rapidly from temperatures well above the 
annealing range, + is very high when atmospheric temperatures are 
reached, and if the cooling is particularly severe, it may be more than 
50° C above the temperatures at which the rapid expansion and endo- 
thermic effects become appreciable; whereas, if a glass is cooled very 
slowly so that it can follow its equilibrium curve as a liquid through- 
out the annealing range and possibly to even lower temperatures, 
Tmay be reduced as much as 70° to 90° C below the beginning points 
of the heat effects mentioned. Somewhere in the range 10° to 40° C 
below these points, the relaxation time is approximately 1 hour, and 
the point at which K=1 is, therefore, roughly the midpoint of the 
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range of possible values of r in glass at ordinary temperatures, Th, 
is, without exceptional treatments, the possible range of r-values ; 
the average glass is about 140° C, and when the glass is in equilibriyy 
(T=r) at a temperature near the midpoint of this range, the viscosity 
coefficient, 7—= E08, (8 expressed in seconds) should be approximate) 
2.210" if H=6X10" at such temperatures. This assumed yaly 
for E is only slightly less than the average value for a number , 
glasses [9] at aia temperatures. Consequently, the above esti. 
mate of the value of 9 is presumably rather high, since E undoubted) 
decreases appreciably as the temperature approaches the annealin: 
range. Although » has been determined experimentally [10, 11, 
temperatures that are so low that they probably were near the uppe 
limits of, or were well within, the annealing ranges of the glasses 
tested, the results cannot safely be compared with the above estima: 
because no auxiliary determinations on the locations of the rapid 
expansion ranges w ~e made. However, the 7-values for the lowes 
temperatures used in these experimental determinations ranged from 
10'* to 10. 

If a slowly cooling glass follows its equilibrium curve throughout 
the range of temperatures corresponding to the range of 7-value 
that may exist even after ordinary temperatures are reached and als 
if k=10° C, the n-values for this interval of cooling will range, a- 
cording to the above estimate for the midpoint value, from 2X10" 
to 210" approximately. However, if the glass is cooled so rapidly 
through this temperature interval that + does not change, and i 
g=h=20° C, the initial 7-value on reaching the lower limit of th 
temperature interval will be about 210", and any stresses in th 
glass will begin to relax almost as rapidly as they would have relaxed 
if the glass were in equilibrium at the midpoint. However, the decre- 
ments in the relaxations as time progresses will be greatly different i 
the two cases. Moreover, if the glass after being brought to equ- 
librium at the lower limit of the temperature interval is heated » 
rapidly to the upper limit that + does not change, the initial -valu 
will again be that of the glass in equilibrium at the midpoint of the 
interval. (These values for g, h, and the resulting n, although they ar 
not improbable, were introduced merely as aids in conceiving th 
effects of a varying or unvarying r.) Effects of this nature, which ar 
the result of the previous history of the glass, were noted and studied 
by the authors of the reports cited. 


III. EQUATION FOR THE RELAXATION OF STRESSES IN 
UNANNEALED GLASS 


In a severely chilled piece of glass, the value of r may be very high 
and the residual stresses, which comprise a balanced system a¢ 
which are the result of the temperature gradients that existed as the 
glass passed through the annealing range, may approach the breaking 
limit of the glass. The high value of + and of these stresses are mall 
tained without appreciable decay at ordinary temperatures becaus 
the viscosity is so extremely high; but, when the glass is reheated, th 
viscosity is so reduced, as the annealing range is approached, that: 
and the average absolute magnitudes, F’, of the stresses immediately 
begin to decrease at an appreciable rate. 

Usually, stressed glasses are double refracting unless the stress 
are the same at all points and in all directions. Except i ceria 
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limited portions, this uniformity of the stresses is never met in glass 
strained by chilling. If the strained piece of glass has some simple 
geometrical form, such as a rectangular prism, it is usually possible to 
pass a beam of plane polarized light through it along a path which 
makes it possible to bisect, throughout most of the path at least, the 
angle between the directions of principal stresses with that of the 
plane of polarization. Because of the photoelastic action of most 
glasses, the component oscillation in the direction of one principal 
stress is retarded (or accelerated) in its velocity of propagation rela- 
tive to the component in the other direction. This relative retarda- 
tion, 6, (usually measured in millimicrons, my, per centimeter of path) 
is proportional to the average difference between the principal stresses, 
which usually vary both in sign arid magnitude along the path. That 
is, it does not measure the maximum difference nor does it give a cer- 
tain indication of the average magnitude of the stresses, particularly, 
if those seldom used paths along which one principal stress is always 
practically zero are excepted. Moreover, as the stresses decay in an 
annealing glass, 6 will not be a certain indicator of the rate of decay 
unless care is taken to preserve the pattern of the stress distribution 
throughout the time required for the test 

Under suitable conditions it may be assumed ‘ that 6 and F are 
proportional throughout an annealing test. That is, as in the case of 
the elastic modulus £, it will be assumed that the changes in the photo- 
elastic constant with 7 and 7 are insignificant compared to the changes 
in » with r and 7. Moreover, experience indicates that 6 may also 
be considered Jas} proportional to the average absolute value of r—T, 
which is usually large and positive in a severely chilled glass if 7’ is not 
too high. That is, z=6/6.=F/Fy)=(r— T)/(r— T), if the conditions at 


any time ¢ are compared to the conditions at the time from which f 
ismeasured. On writing A=(r—7),/h and substituting hAz for r—T 
in eq 1, the relation, d(Azx)/dt= —(Kge™) Azxe4*, is obtained for 
annealing at a constant temperature. The solution of this differential 
equation 1s 


Fi(— Az) — Ei(— A) = — Kee” “t= — Ket, (3) 


since r=1 when t=0. Clearly, Kr is the value of the reciprocal of the 
relaxation time that will be reached provided the annealing is con- 
tinued until r=T7Z. In certain cases this condition exists from the be- 
ginning and Maxwell’s equation, dF /dt=—K,F, applies. (In the case 
of a balanced system dS/dt=0, since the sum of the elastic and in- 
elastic strains remains constant.) That is, in view of the assumed 
proportionality of F and 4, 


In (6/59) = Art. (4) 


Moreover, jn exceptional cases, 7’ exceeds the average value of +r 
In such cases A is negative. 


No doubt this assumption merely approximates the actual relation between the relaxing stresses and the 
anging fictive temperature. In regard to this relation, it will be noted that the density changes that ac- 
ny the changes in r will react on the stresses if r-gradients exist in the glass. That is, appreciable 
nts doubtless are present in any severely strained glass, and the accompanying differential density 
hould as a rule hasten the relaxation. This follows because the highest values of r and the conse- 
hest rate in its decrease are presumably near the surface where the compressive stresses approach 
‘maximum. The comparatively large increase caused in the density by the decrease in r in the region of 
Hese Stresses should reduce them to a greater extent than the tractions in the interior are enhanced by; lesser 
nereases in density. 


‘ 
b 
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IV. APPLICATION OF THE EQUATION TO 
EXPERIMENTAL DATA 


To procure data that have the required accuracy for thoroughly 
testing eq 3 or any other proposed expression concerning the decay 
of the double refraction in annealing glass is not an easy matter 
In the first place, sizeable pieces of glass are usually considera) 
necessary for the tests. Consequently, temperature gradients of 
considerable magnitude are likely to be present at the time of testing 
unless exceptional care is taken. These gradients make the effectiy, 
temperature of the glass rather uncertain. Also, they cause th 
stresses to relax more quickly in some parts of the glass than in other 
Through this differential action, the double refraction pattern 
gradually changed, and, as a result, 6 may actually increase for , 
time if the observations are made along some possible paths through 
the glass. In the second place, K changes so rapidly with temperatur 
that the temperature range available for the tests is rather limited 
That is, the decay is so rapid at the upper limit of the range that 
much of the double refraction effect has disappeared before the glass 
can be heated to that point; while it is so slow at the lower limi 
that the tests are extremely tedious unless the strains are large 
As a result, the observations obtained in a test at high temperature 
are usually few in number and made when 6 is small, whereas thos 
obtained at low temperatures, although comparatively numerous 
are made when 6 is very much larger. In view of the dependen 
of K on 7, data of this sort will yield values for k that are too larg 
by varying degrees if the simple Maxwellian equation is used. More- 
over, because of the limited range available, the temperature intervals 
between tests designed for the determination of k are often scarcely 
more than double the possible errors in the temperature deter 
minations. 

Presumably, because of the effect of a changing fictive temper 
ture, r, Adams and Williamson found that eq 4 was inapplicable to 
their data and, therefore, proposed the reciprocal relation 


1/5—1/8)=A’t. 


Consequently, to test the relative applicability of eq 3, 4, and 5 to 
the decay of double refraction in annealing glass, they were applied 
to the data obtained by Adams and Williamson (page 622, reference? 
on a light flint while it was being annealed at 300° C. The exper 
mental data and resulting computed curves are shown in figure |. 
In fitting these data, eq 3 has the obvious advantage of an addi- 
tional constant. In arriving at tentative values for the three con- 
stants, 5, A, and Ky, required in using this equation, the method 
followed was essentially, (a) the drawing of a smooth curve that seem 
ingly approximated the trend of the experimental data, (b) the choos 
ing of three trial points on this curve in order to obtain two equations 
for the purpose of eliminating Kz, and finally (c) the search for a valu 
of A that would satisfy the equation resulting after the elimination. 
A few trials will yield such a value if an adequate table of exponent! 
integrals is used. In the case of the above-mentioned data for! 
light flint, the first and third trial points coincided, respectively, wit 
the first and eleventh experimental points, while the second was nea! 
the third experimental point, These trial points yielded the values 
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5.56 and 6.97X10-° for A and Ky, respectively, and the computed 
curve resulting from the use of these constants is that designated 
by Ain figure 1. Since this curve fits the experimental data reason- 
ably well, a betterment was not attempted. 

As Adams and Williamson pointed out, eq 5 does not yield a 
reasonable fit for these particular experimental data, and eq 4 is even 
less satisfactory as curves B and C, respectively, indicate. In com- 
puting these curves, only two trial points are required to determine 
tentative values for 4d) and A’ in one case and 6, and Ky, in the other. 
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Figure 1.—Decay of double refraction in light flint annealing at 300° C. 
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ve A: Circles, experimental data by Adams and Williamson; crosses, points computed by equation 3. 
Curve B: Points computed by equation 5. 
Curve C: Points computed by equation 4. 


The trial points used in these cases were the first and third of those 
used previously. 

In addition to this test at 300° C, Adams and Williamson deter- 
mined the decay of double refraction in a light flint at three higher 
temperatures that are within the upper half of the annealing range. 
Table 1 presents the residuals (R;) obtained through the use of eq 3 
ind also those (Rs) obtained by Adams and Williamson through the 
use of eq 5. In applying the first of these equations to their data, it 
usually was endeavored to balance approximately the positive and 
hegative residuals both in magnitude and distribution, also to make 
those residuals as few as possible which exceeded an assumed maxi- 
um probable error of +5my in the 6-determinations. In only a 
few cases was it necessary to depart markedly from these aims, pro- 
vided the results on ordinary crown are excepted. 
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For simplicity in applying eq 3, the observed initial relative retgy. 
dation was chosen as the value of 4. The ratios z=6/9 for th 
remaining observations were then multiplied by various trial valye 
for A and corresponding sets of values for Krt and Kr were computed 
The A-value that yielded a spread that seemingly approximated th 
minimum possible in the K,-values and the average of its set of 
K,-values were then used to compute values of 6 and the residuals 
(6,—6,) for the various times of observation. Usually it was neces. 
sary to adjust these A- and K7-values a few times in order to procur 
more satisfactory sets of differences between the experimental (j, 
and computed (6,) retardations. 


TABLE 1. Comparsien wh re siduals obtained 


LIGHT FLINT 
T=385°C, d=415, Kr=0.21, . 


0. 000 | 0. 183 0. 367 0. 567 
415 280 223 1S8 
—15 +16 +12 +4 | 

0 —4 —2/ —! 


T=400°C, d= 280, Kr=0.60, / 

ay 2 ws yy | 
0.000} 0.167| 0.433 0. 733 | 
280 195 | 135 97 | 
+9 —1 | —8 | —6 | 


0 +1 | —2 -2 





T=415°C, b= 293, Kr=1.74, A=1.5 
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In particular, the results obtained on ordinary crown at 450° C 
offered exceptional difficulties, since there is evidence of a serious 
temperature change after the eleventh hour of this test. As a matter 
of fact, the authors state that most of the measurements on this 
glass were of a preliminary nature and, from their remarks, it may be 
assumed that some of the results possibly are affected by considerable 
errors in the temperature determinations. Furthermore, most of the 
relative retardations were determined by an inferior method in the 
case of this glass. 

In addition to the light flint, Adams and Williamson also tested a 
medium flint not only at points restricted mainly to the upper half of 
the annealing range but also at a temperature, 306° C, which is far 
below that range as it is limited by normal practice. For this reason, 
the results obtained by applying eq 3 to the data on medium flint also 
are included in table 1. In the case of both of these glasses it will be 
noted that the averages of the residuals resulting from the use of 
eq 3 usually are less than one-half the averages of the residuals obtained 
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Figure 2.—Natural logarithms of Kr for various glasses at various temperatures. 
The values of Kr were computed by equation 3 from data obtained by Adams and Williamson for various 

flasses annealing at various temperatures. } 

_ Curve 1, extra-heavy flint; 2, light flint; 3, medium flint; 4, heavy flint; 5, barium flint; 6, ordinary crown; 
borosilicate crown; 8, light barium crown; 9, dense barium crown. 


by Adams and Williamson from the use of eq 5. Usually, this experi- 
ence was repeated when eq 3 was applied to data for the seven other 
glasses which these authors tested. Consequently, the residuals for 
these remaining glasses are not presented. However, the natural 
logarithms of the values for Ky are indicated in all cases in figure 2. 


V. REMARKS CONCERNING THE VALUES OBTAINED 
FOR A, Kr, AND k 


When 4) was large and the testing temperature was comparatively 
low, the values found for A ranged from 5.6 to 3, whereas, when 4, 
was relatively low and the testing temperature was high, they ranged 
fom 1.5 to —0.6. Otherwise the A-values were intermediate and 
came in the order anticipated except in a few cases, such as the values 
obtained from the tests on medium flint at 390° and 406° C (see table 
!). However, such reversals should not be uncommon, since a glass in 


627528—45——_8 
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equilibrium at any point in the upper half of the annealing range eq) 
be cooled sharply enough to be as severely strained as a glass that has 
been cooled less sharply from a higher temperature but still sharply 
enough to prevent a too great recession in r. , 

With regard to small values of A, two minor points observed whilp 
applying eq 3 to these data deserve mention. The first is concerned 
with the paucity of data on light barium crown in the test at 540° ¢ 
Two observations are sufficient to yield a value for Kr or A’ if eq 4 
or 5 is used, but are inadequate for the determination of A and K; of 
eq 3. Since 162 my is not a very large value for 4, (in view of the 
7.5-em path in this case), it seemed possible that eq 4 might be ade. 
quate, but its use yielded the result K7= 1.38, which is somewhat ou: 
of line with the values of Ky (0.066 and 0.264) that were obtained 
from the results of the tests at 500° and 516° C, respectively. That is 
these results indicate that a better value for Ky at 540° C is 2) 
although the introduction of this value yields the result A= —0.565 
when eq 3 is applied to the data obtained at 540° C. Since this tem. 
perature is not particularly high, a negative value for A, in this case 
seems somewhat improbable, but it calls attention to the fact that 
such values should be common whenever the stresses are relaxing at 
an annealing temperature that is higher than others at which tly 
glass just previously was annealing. For instance, in heating a 
strained glass at a moderate rate through the annealing range, r—T 
decreases not only because + is decreasing but also because T is 
rising, and it, consequently, often becomes negative while the stresses 
are still rather large. 

The second point concerns the results A=0.2 and Kr=0.4 obtained 
for the test on ordinary crown at 500° C and in which 6)=120 mg 
When these values were used, the residuals were 0, 0, —2, and 
in the order of decreasing 6. The low value of A indicates that eq 4 
should apply almost equally well, and that this is the case is shown by 
the residuals 0, +2, —1, and 0, which are obtained through its use if 
Kr is given the value 0.45. The corresponding residuals obtained by 
Adams and Williamson are +3, —6, —2, and +10. In this case also 
the number of observations is unsatisfactory and the testing tempera- 
ture, 500° C, is scarcely high enough to cause a close approach to the 
condition, (r— 7),.=0, ordinarily; but the result calls attention to the 
fact that eq 4 alone applies if A approximates zero. 

Although it would be quite difficult, if the testing temperatures 
were in the upper half of the annealing range, it is still possible to % 
manipulate the chilling treatments and the reheating to the annealing 
temperatures that the values of (r— 7), are known, at least approx 
mately, in tests such as these that are being discussed. In that case 
when A is determined, h becomes known. Unfortunately, it is impos 
sible to determine the A-values for any of the glasses tested by Adams 
and Williamson because there is nothing in their data to indicate the 
values of 7». ; 

When the natural logarithms of the values found for Kr are plotted 
as ordinates against the temperatures of the tests, as in figure 2, the 
result is much the same as that shown by Adams and Williamson 1 
a similar presentation of the relation that their A’-values bore to the 
same temperatures. Since the lines representing the trend of the 
plotted points in figure 2 were located graphically, judgment playeé 
a considerable role in their positions, curvature and slopes. ine 
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they are all approximately straight, these lines indicate that the 
j-values of the glasses do not change materially within the annealing 
ranges. Consequently, since the relaxation time, v=1/Kr, has been 
expressed in hours, k=(7'—T;)/In Kr approximately if the relaxation 
time is 1 hour at 7. In preparing these curves and computing the 
following values of k, only two determinations of In Ky were ignored. 
These were the results derived from the data on barium flint and on 
ordinary crown at the temperature 490° C in both cases. These 
results were ignored because it seemed possible that errors of 20° were 
made in recording the temperatures. 

The values obtained for & are as follows: 11.8° C, extra heavy flint; 
12.6°, light flint; 10.2°, heavy flint; 12.8°, medium flint; 11.6°, barium 
flint; 13.4°, ordinary crown; 10.5°, borosilicate crown; 10.6°, light 
barium crown; and 12.9°, dense barium crown. Although the prob- 
able errors in these determinations are undoubtedly large, the average 
of the k-values is 11.8° C, which agrees closely with Twyman’s experi- 
ence. From the change in the constant A’, Adams and Williamson 
obtained an average value for their constant, M4,=(logye)/k, that 
approximates 0.33. This corresponds to a k-value approximating 
13.2°. However, presumably for the purpose of simplifying the com- 
putations in preparing their sample annealing schedules, they gave 
p=(logio2)/M, the value 10°, in a relation which is essentially A’= 

(7-T)/? rather than 9.1°, which their average value for M, 
indicates or the value, 8°, which was found by Twyman. Compara- 
tively few if any of the optical glasses ordinarily produced have a 
k-value as high as 14.4°, which corresponds to p=10°. This rather 
high value for p suggests that the annealing times computed by Adams 
and Williamson are too short and particularly so for the comparatively 
low annealing temperatures. However, it must be considered that 
they were merely demonstrating a method of preparing annealing 
schedules and were not recommending schedules that could be used 
for all of the glasses that may be included in the types tested. This 
is obvious, since it is well known that the upper limits of the ranges 
of suitable annealing temperatures for the various glasses that may 
be included within a single type often cover a range of 50° C or more. 

Dimensionally, both AK and A’L, when L represents the length of 
the light path in glass, are time reciprocals and are directly related to 
the fluidity of the glass. The minute was the time unit used by Adams 
and Williamson in making their observations. Consequently, at the 
temperature 7 at which 1/A’Z becomes 1 minute, M,=M,T= 
T logiwe)/k. Such temperatures, 7, are far above the annealing range 
and are high enough to cause considerable deformation in a short time. 
Within the annealing range, the hour is a suitable unit for measuring 
the relaxation time, 8=1/K. In that case, the temperatures, 7, at 
which K=1 and In K=0, are suitable for annealing optical blanks 
that are not too large. These temperatures can be determined from 
igure 2 and are as follows: 390° C, extra-heavy flint; 407°, light 
lint; 405°, heavy flint; 413°, medium flint; 475°, barium flint; 513°, 
ordinary crown; 532°, borosilicate crown; 532°, light barium crown; 
and 585°, heavy barium crown. As the annealing temperatures of 
pieces of glass with a thickness of 2 cm, Adams and Williamson 
recommended the following corresponding temperatures: 382° C, 
408", 407°, 418°, 494°, 514°, 541°, 552°, and 592° C. 

Finally, the results obtained for K at very low temperatures, in the 
cases of the light and medium flints, indicate that * decreases appre- 
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ciably from its magnitude at ordinary annealing temperatures whep. 
ever a glass is cooled to temperatures near or below the lower limit 
of its practical annealing range. That is, the K-values obtained fo 
these glasses indicate that curves representing In K, as in figure 2, shoul 
become appreciably steeper at very low annealing temperatures 
Such a decrease is to be expected, since it seems probable that ther 
ig a continual decrease in the available flow planes as a glass coo); 
through these low temperatures and gradually changes from a highly 
viscous liquid condition to an elastic solid condition in which it shows 
no more viscous and plastic deformability than many crystallin. 
solids. It will be noted, however, that glasses that have been cooled 
very rapidly from high temperatures should show, because of their 
high r-values, a much greater inelastic deformability than well 
annealed glasses. In fact, it should be possible to chill fine threads 
of glasses or plastics, especially those with very low annealing ranges, 
so severely that they will anneal at atmospheric temperatures until 
r decreases sufficiently to raise the viscosity and plasticity beyond 
that possibility. 
VI. CONCLUSION 


When a glass is cooled very slowly from high temperatures, it con- 
tinues to behave as a liquid throughout the annealing range and ap- 
parently to even lower temperatures if the cooling rate is extremely 
low. Naturally, its properties, including viscosity, all change a- 
cordingly, that is, their temperature coefficients are generally much 
higher than when the glass or plastic behaves as a solid at tempers 
tures far below the annealing range. 

When the cooling rate is comparatively high, the glass departs from 
the equilibrium curve of its liquid condition at temperatures that are 
determined by the rate and nature of the glass. These temper 
tures of departure will be above the usual annealing range if the rate is 
exceptionally high. If the rapid cooling is continued after such de 
partures, the glass begins to behave as a solid even before the a- 
nealing range is left behind. That is, even within the annealig 
range, a glass, although greatly undercooled, will react temporarily a 
a solid to comparatively sudden temperature changes. 

Far below their annealing ranges, glasses and many of the s0- 
called plastics behave as solids; but, when they are reheated, they 
drift toward the equilibrium curves of their liquid conditions as soot 
as the annealing range is reached. When the undercooling is con- 
siderable and the rate of reheating is comparatively low, this drift » 
very marked and often begins 50° C or more below the lowest a- 
nealing temperatures. Drifts of this kind are accompanied by larg’ 
contraction and exothermic effects, which practically cease as soon ®& 
the equilibrium curve of the liquid condition is neared. 

If the heating is continued after the equilibrium curve is reachet 
the glass becomes superheated, and this superheating will be quite 
large if the heating rate is high enough and the prior cooling rate ws 
comparatively low. However, the rising temperature ultimately tf 
duces the viscosity so greatly that the glass or plastic returns rapid 
to the equilibrium curve, and this return is accompanied by lam 
“rapid expansion” and endothermic effects. 

When the rate of reheating is much higher than the prior rate 0! 
cooling, the drift while the glass is in an undercooled condition § 
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quite negligible, and the curve followed by a glass or plastic while in 
ts solid condition at comparatively low temperatures continues beyond 
the point at which it intersects the curve of the equilibrium condition. 
This point or temperature of intersection has been referred to in pre- 
rious papers by such terms as the “‘effective annealing,” “equilibrium,” 
or “fietive’’ temperature of the glass while it is in the particular solid 
condition corresponding to a particular intersection point. This tem- 
perature, Which can be varied at will over a wide range by reheating 
the glass, was termed the fictive temperature because it is, in most 
cases, at least as effective in determining the properties of a glass as 
the actual temperature. 

Although most of the properties of glass are affected considerably 
by any large change in the fictive temperature, any rise or fall in it 
causes an exceptionally marked decrease or increase, respectively, in 
the viscosity and, consequently, greatly affects the relaxation of any 
stresses that may be in the glass. This leads to an interesting situa- 
tion, since the fictive temperature is usually changing continually even 
when the stresses are relaxing at a constant real temperature. 

Based on experience, an empirical expression was proposed pre- 
viously for the relation between the time rate of change in this tem- 
perature and the magnitudes of it ‘and the real temperature. In this 
report, this expression has been used in connection with Maxwell’s 
equation for the development of inelastic deformation (that is, for the 
flow in stressed visco-elastic materials) in order to obtain an expression 
for the relaxation of stresses in an annealing glass. The resultant 
expression has been applied to the data obtained by Adams and 
Williamson on the stress relaxation at different annealing tempera- 
tures and for different kinds of glasses, and it has been found to fit 


experimental data better than previously suggested expressions fit 
them. 
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